Next Article in Journal
Magnetic Trails: A Novel Artificial Pheromone for Swarm Robotics in Outdoor Environments
Previous Article in Journal
Pattern Recognition for Human Diseases Classification in Spectral Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Stability and Numerical Scheme of Fractional Differential Equations with Application to Biology

1
Equipe de Recherche en Modélisation et Enseignement des Mathématiques (ERMEM), Centre Régional des Métiers de l’Education et de la Formation (CRMEF), Derb Ghalef, Casablanca 20340, Morocco
2
Laboratory of Analysis, Modeling and Simulation (LAMS), Faculty of Sciences Ben M’Sick, Hassan II University of Casablanca, Sidi Othman, Casablanca P.O. Box 7955, Morocco
Submission received: 24 April 2022 / Revised: 2 June 2022 / Accepted: 10 June 2022 / Published: 15 June 2022
(This article belongs to the Section Computational Biology)

Abstract

:
The fractional differential equations involving different types of fractional derivatives are currently used in many fields of science and engineering. Therefore, the first purpose of this study is to investigate the qualitative properties including the stability, asymptotic stability, as well as Mittag–Leffler stability of solutions of fractional differential equations with the new generalized Hattaf fractional derivative, which encompasses the popular forms of fractional derivatives with non-singular kernels. These qualitative properties are obtained by constructing a suitable Lyapunov function. Furthermore, the second aim is to develop a new numerical method in order to approximate the solutions of such types of equations. The developed method recovers the classical Euler numerical scheme for ordinary differential equations. Finally, the obtained analytical and numerical results are applied to a biological nonlinear system arising from epidemiology.

1. Introduction

Fractional differential equations (FDEs) are the differential equations with non-integer powers of the differentiation order. In recent years, FDEs have gained importance in both theoretical and applied aspects of several fields of science and engineering such as biology [1,2], epidemiology [3,4,5], control theory [6], viscoelasticity [7], engineering [8], and bioengineering [9]. For other recent developments in the literature on theoretical and numerical studies of FDEs, see for example [10,11,12] for the inverse problem associated with FDEs, ref. [13] for Hermite–Hadamard-type inequalities, [14,15,16] for numerical methods of FDEs, and [17,18,19] for analytical and numerical solutions of some FDEs.
Currently, stability analysis of FDEs has been investigated by many authors. Li et al. [20] dealt with the stability of nonlinear dynamic systems describing by the Caputo fractional derivative with a singular kernel [21]. Delavari et al. [22] analyzed the stability of fractional-order nonlinear systems. They presented an extension of the Lyapunov direct method for Caputo-type fractional-order systems using Bihari’s and Bellman–Gronwall’s inequalities [23,24]. The stability of FDEs with Hadamard fractional derivative [25] was investigated by Wang et al. [26] utilizing a new fractional comparison principle. In [27], the authors focused on the Ulam stability of a generalized delayed differential equation of fractional order. A more recent study presented in [28] discussed the Mittag–Leffler stability of FDEs using the new generalized Hattaf fractional (GHF) derivative [29], which includes many fractional derivatives available in the literature such as the Caputo–Fabrizio fractional derivative [30], the Atangana–Baleanu fractional derivative [31], and the weighted Atangana–Baleanu fractional derivative [32]. The stability in the sense of Ulam–Hyers of FDEs with GHF derivative was studied in [33] using a new version of the Gronwall inequality. However, this paper extends the fractional comparison principle to the GHF derivative and the results related to Mittag–Leffler stability given in [28] using class K functions. Furthermore, the present paper establishes new interesting results concerning the stability, as well as the asymptotic stability of FDEs with the GHF derivative by means of the Lyapunov direct method and class K functions.
On the other hand, numerical methods have become indispensable tools to find the approximate solutions of both ordinary differential equations (ODEs) and FDEs. Oftentimes, it is impossible or complicated to find the exact solution for many nonlinear systems modeling real phenomena. Therefore, the second main objective of this research is to develop a new numerical method for solving FDEs with the GHF derivative.
The rest of the paper is outlined as follows. Section 2 describes the basic concepts and extends the fractional comparison principle to the GHF derivative. Section 3 analyzes the stability, the asymptotic stability, and the Mittag–Leffler stability of FDEs with the GHF derivative. Section 4 proposes a new numerical method to solve nonlinear FDEs. Section 5 presents an application of our analytical and numerical results to a biological system. Finally, the conclusion is given in Section 6.

2. Preliminaries

In this section, we present the necessary concepts and results related to the GHF derivative that are used throughout this paper.
Definition 1
([29]). Let α [ 0 , 1 ) , β , γ > 0 , and f H 1 ( a , b ) . The GHF derivative of order α in the Caputo sense of the function f ( t ) with respect to the weight function w ( t ) is defined as follows:
C D a , t , w α , β , γ f ( t ) = N ( α ) 1 α 1 w ( t ) a t E β [ μ α ( t τ ) γ ] d d τ ( w f ) ( τ ) d τ ,
where w C 1 ( a , b ) , w , w > 0 on [a,b], N ( α ) is a normalization function obeying N ( 0 ) = N ( 1 ) = 1 , μ α = α 1 α , and E β ( t ) = k = 0 + t k Γ ( β k + 1 ) is the Mittag–Leffler function of parameter β.
Some examples of normalization functions are as follows:
  • N ( α ) = 1 ,
  • N ( α ) = 1 α + α Γ ( α ) .
The main properties of the GHF derivative defined by (1) are given in detail in [29,34]. Furthermore, (1) is reduced to the Caputo–Fabrizio fractional derivative [30] when w ( t ) = 1 and β = γ = 1 , to the Atangana–Baleanu fractional derivative [31] when w ( t ) = 1 and β = γ = α , as well as to the weighted Atangana–Baleanu fractional derivative [32] when β = γ = α .
Now, denote C D a , t , w α , β , β by D a , w α , β . By [29], the generalized fractional integral associated with D a , w α , β is given by the following definition.
Definition 2
([29]).The generalized fractional integral operator associated with D a , w α , β is defined by
I a , w α , β f ( t ) = 1 α N ( α ) f ( t ) + α N ( α ) R L I a , w β f ( t ) ,
where R L I a , w β is the standard weighted Riemann–Liouville fractional integral of order β defined by
R L I a , w β f ( t ) = 1 Γ ( β ) 1 w ( t ) a t ( t τ ) β 1 w ( τ ) f ( τ ) d x .
Now, we recall an important theorem that we will need in the following. This theorem extends the Newton–Leibniz formula introduced in [35,36].
Theorem 1
([34]).Let α [ 0 , 1 ) , β > 0 and f H 1 ( a , b ) . Then, we have the following properties:
I a , w α , β D a , w α , β f ( t ) = f ( t ) w ( a ) f ( a ) w ( t ) ,
and
D a , w α , β I a , w α , β f ( t ) = f ( t ) w ( a ) f ( a ) w ( t ) .
In addition, we will need the following definition and lemma.
Definition 3.
A continuous function ψ : [ 0 , + ) [ 0 , + ) is said to belong to class K if it is strictly increasing and ψ ( 0 ) = 0 .
Lemma 1.
(Fractional comparison principle) Let x ( t ) and y ( t ) two functions defined on [ t 0 , + ) with D t 0 , w α , β x ( t ) D t 0 , w α , β y ( t ) and x ( t 0 ) y ( t 0 ) . Then, x ( t ) y ( t ) , for all t t 0 .
Proof. 
We have D t 0 , w α , β x ( t ) D t 0 , w α , β y ( t ) . By applying the fractional Hattaf integral to both sides of this inequality and using (4), we obtain
x ( t ) x ( t 0 ) w ( t 0 ) w ( t ) y ( t ) y ( t 0 ) w ( t 0 ) w ( t ) ,
which leads to
x ( t ) y ( t ) + w ( t 0 ) x ( t 0 ) y ( t 0 w ( t ) .
Since x ( t 0 ) y ( t 0 ) , we deduce that
x ( t ) y ( t ) , for all t t 0 .
This completes the proof. □
Remark 1.
Lemma 1 extends two results, one presented in Theorem 2.4 of [22] with the Riemann–Liouville fractional derivative and the other in Lemma 6.1 of [20] with the Caputo fractional derivative.

3. Stability of FDEs with the GHF Derivative

In this section, we study the stability of the following nonautonomous FDE with the GHF derivative expressed by
D 0 , w α , β x ( t ) = f ( t , x ( t ) ) ,
where x ( t ) I R n is the pseudo-state variable, f : [ 0 , + ) × Ω I R n is a continuous locally Lipschitz function satisfying in particular f ( t , 0 ) = 0 , and Ω is a domain of I R n that contains the origin x = 0 . When f ( t , x ) = f ( x ) , (6) becomes an autonomous FDE of the form
D 0 , w α , β x ( t ) = f ( x ( t ) ) .
Since α [ 0 , 1 ) , we chose x ( 0 ) = x 0 as the initial condition for (6) and (7).
First, we give some definitions of stabilities that will be used in the remainder of this paper. We begin with the definition of stability and asymptotic stability.
Definition 4.
Let x = 0 be an equilibrium point for the system (6):
(i) 
The equilibrium point x = 0 is said to be stable if, for any ϵ > 0 , there exists a η > 0 such that for each initial condition x ( t 0 ) = x 0 satisfying x 0 < η , the solution x ( t ) of (6) satisfies x ( t ) < ϵ for all t t 0 . Otherwise, we say that x = 0 is unstable.
(ii) 
The equilibrium point x = 0 is said to be asymptotically stable if it is stable and lim t + x ( t ) = 0 .
For the Mittag–Leffler stability, we have the following definition.
Definition 5.
The trivial solution of (6) is said to be Mittag–Leffler stable if
x ( t ) [ m x ( t 0 ) E β ( λ ( t t 0 ) β ) ] ν ,
where t 0 is the initial time, λ 0 , ν > 0 , m ( 0 ) = 0 , m ( x ) 0 , and m ( x ) is locally Lipschitz on x I R n with the Lipschitz constant m 0 .
It is very important to note that the Mittag–Leffler stability implies asymptotic stability.
Theorem 2.
Let x = 0 be an equilibrium point for the system (6). If there exist a continuously differentiable function V ( t , x ) : [ 0 , + ) × Ω I R and a class K function ψ satisfying the following conditions:
V ( t , x ) ψ ( x ) ,
D 0 , w α , β V ( t , x ) 0 ,
then x = 0 is locally stable. If (8) and (9) hold globally on I R n , then x = 0 is globally stable.
Proof. 
According to (9) and Theorem 1, we have
V ( t , x ( t ) ) V ( 0 , x ( 0 ) ) w ( 0 ) w ( t ) .
Since w ( 0 ) w ( t ) for all t 0 , we obtain:
V ( t , x ( t ) ) V ( 0 , x ( 0 ) ) , t 0 .
It follows from (8) that
x ( t ) ψ 1 V ( 0 , x ( 0 ) ) , t 0 .
Then, the equilibrium x = 0 is stable. □
Remark 2.
Theorem 2 extends the asymptotical stability result with the Caputo derivative given in Theorem 3.2 of [37] to the GHF derivative.
Theorem 3.
Let x = 0 be an equilibrium point for the system (6). If there exist a continuously differentiable function V ( t , x ) : [ 0 , + ) × Ω I R and class K functions ψ i (i = 1, 2, 3) satisfying
ψ 1 ( x ) V ( t , x ) ψ 2 ( x ) ,
D 0 , w α , β V ( t , x ) ψ 3 ( x ) ,
then x = 0 is asymptotically stable.
Proof. 
By (11) and (12), we have
D 0 , w α , β V ( t , x ) ψ 3 ψ 2 1 V ( t , x ) .
It follows from the fractional comparison principle presented in Lemma 1 that V ( t , x ) is bounded by the nonnegative solution of the following FDE:
D 0 , w α , β y ( t ) = ψ 3 ψ 2 1 y ( t ) , y ( 0 ) = V ( 0 , x ( 0 ) ) .
We have D 0 , w α , β y ( t ) 0 . Then,
y ( t ) y ( 0 ) , t 0 .
If y ( 0 ) = 0 , then y ( t ) = 0 for all t 0 .
If y ( 0 ) 0 , then lim t + y ( t ) = 0 . Assume the contrary, then there exists a ϵ > 0 such that y ( t ) ϵ for all t 0 . We have
D 0 , w α , β y ( t ) = ψ 3 ψ 2 1 y ( t ) ψ 3 ψ 2 1 ( ϵ ) = ψ 3 ψ 2 1 ( ϵ ) y ( 0 ) y ( 0 ) .
Then
D 0 , w α , β y ( t ) λ y ( t ) ,
where λ = ψ 3 ψ 2 1 ( ϵ ) y ( 0 ) > 0 . By applying Corollary 1 of [28], we obtain
y ( t ) y ( 0 ) E β α λ t β N ( α ) + ( 1 α ) , t 0 .
This is a contradiction. Therefore, it follows from the fractional comparison principle that
V ( t , x ( t ) ) y ( t ) .
By (11), we obtain
x ( t ) ψ 1 1 y ( t ) ,
which implies that
lim t + x ( t ) = 0 .
Then, x = 0 is asymptotically stable. □
Remark 3.
Theorem 3 extends the result of the asymptotic stability with the Caputo fractional derivative introduced in Theorem 6.2 of [20] to the GHF derivative with a nonsingular kernel.
Theorem 4.
Let x = 0 be an equilibrium point for the system (6). If there exist a continuously differentiable function Let V ( t , x ) : [ 0 , + ) × Ω I R such that is locally Lipschitz with respect to x and a class K function ψ satisfying:
k 1 x p V ( t , x ) k 2 ψ ( x ) ,
D 0 , w α , β V ( t , x ) k 3 ψ ( x ) ,
where k 1 , k 2 , k 3 , and p are arbitrary positive constants, then x = 0 is Mittag–Leffler stable. If (16) and (17) hold globally on I R n , then x = 0 is globally Mittag–Leffler stable.
Proof. 
It follows from (16) and (17) that
D 0 , w α , β V ( t , x ) k 3 k 2 V ( t , x ) .
According to Corollary 1 of [28], we obtain
V ( t , x ) V ( 0 , x ( 0 ) ) E β η t β ,
where η = k 3 α k 2 N ( α ) + k 3 ( 1 α ) . Using (16), we obtain
k 1 x ( t ) p V ( 0 , x ( 0 ) ) E β η t β ,
which leads to
x ( t ) m x ( 0 ) E β η t β 1 p ,
where m ( x ) = V ( 0 , x ) k 1 . Therefore, the equilibrium x = 0 is Mittag–Leffler stable. □
Remark 4.
Theorem 4 generalizes the result concerning the Mittag–Leffler stability presented in Theorem 3 of [28]. It suffices to take ψ ( z ) = z q , where z [ 0 , + ) and q is an arbitrary positive constant.
Theorem 5.
Let x = 0 be an equilibrium point for the autonomous system (7) and V ( x ) be a continuously differentiable function in a neighborhood U I R n of the origin satisfying the following conditions:
(i) 
V ( 0 ) = 0 and V ( x ) > 0 for all x U { 0 } ;
(ii) 
D 0 , w α , β V ( x ) 0 for all x U { 0 } .
Then, x = 0 is stable.
Proof. 
Let ϵ > 0 such that B ( 0 , ϵ ) ¯ U , where B ( 0 , ϵ ) ¯ denotes the closed ball with center 0 and radius ϵ be defined by B ( 0 , ϵ ) ¯ = { x I R n : x ϵ } . Furthermore, we define the open ball with center 0 and radius ϵ by B ( 0 , ϵ ) = { x I R n : x < ϵ } .
Since V is continuous on the compact subset S ( 0 , ϵ ) = { x I R n : x = ϵ } , we deduce there exists a x ^ S ( 0 , ϵ ) such that
V ( x ^ ) = inf x S ( 0 , ϵ ) V ( x ) = ϱ .
According to (i), we have ϱ > 0 . Consider the following subset of U:
U 1 = { x B ( 0 , ϵ ) : V ( x ) < ϱ } .
Then, U 1 is a neighborhood of the origin because it is an open ball containing 0.
Let ϕ ( t , x 0 ) be a solution of (7) with initial condition ϕ ( 0 , x 0 ) = x 0 U 1 . According to (ii), we have D 0 , w α , β V ϕ ( t , x 0 ) 0 . Then,
V ϕ ( t , x 0 ) V ϕ ( 0 , x 0 ) = V ( x 0 ) , for all t 0 .
Hence,
ϕ ( t , x 0 ) S ( 0 , ϵ ) , for all t 0 .
Indeed, suppose the contrary. Therefore, there exists a t 1 0 such that ϕ ( t 1 , x 0 ) S ( 0 , ϵ ) . This implies that V ϕ ( t 1 , x 0 ) ϱ . Then, V ( x 0 ) ϱ , which is contradicted with x 0 U 1 .
Now, we prove that
ϕ ( t , x 0 ) B ( 0 , ϵ ) ¯ , for all t 0 .
In fact, assume the contrary. Then, there exists a t 2 0 such that ϕ ( t 2 , x 0 ) B ( 0 , ϵ ) ¯ , which implies that ϕ ( t 2 , x 0 ) > ϵ . Let Λ = { t 0 : ϕ ( t , x 0 ) > ϵ } . Hence, there exists a sequence ( t n ) in { t 0 : ϕ ( t , x 0 ) > ϵ } such that lim n + t n = Λ . Thus, lim n + ϕ ( t n , x 0 ) ϵ . Therefore,
ϕ ( Λ , x 0 ) ϵ .
For t < Λ , we have ϕ ( t , x 0 ) ϵ . Then, lim t Λ ϕ ( t , x 0 ) ϵ , which implies that
ϕ ( Λ , x 0 ) ϵ .
From (23) and (24), we obtain ϕ ( Λ , x 0 ) = ϵ . This contradicts (21). We conclude that if any initial condition of System (7) satisfies x 0 U 1 , then the solution of (7) satisfies ϕ ( t , x 0 ) B ( 0 , ϵ ) ¯ , for all t 0 . This implies that the equilibrium x = 0 of System (7) is stable. □

4. Numerical Scheme

In this section, we introduce a numerical method to approximate the solution of the FDE with the GHF derivative given in (6).
From Theorem 1, Equation (6) can be converted to the following fractional integral equation:
x ( t ) x ( 0 ) w ( 0 ) w ( t ) = 1 α N ( α ) f t , x ( t ) + α N ( α ) Γ ( β ) 1 w ( t ) 0 t ( t τ ) β 1 w ( τ ) f τ , x ( τ ) d τ .
Let t n = n h , where n I N and h is the time step duration. We have
x ( t n + 1 ) = x 0 w ( 0 ) w ( t n ) + 1 α N ( α ) f t n , x ( t n ) + α N ( α ) Γ ( β ) w ( t n ) 0 t n + 1 ( t n + 1 τ ) β 1 w ( τ ) f τ , x ( τ ) d τ .
By applying the rectangular integration to the integral in the right-hand side of (26), we obtain
0 t n + 1 ( t n + 1 τ ) β 1 w ( τ ) f τ , x ( τ ) d τ = k = 0 n t k t k + 1 ( t n + 1 τ ) β 1 w ( τ ) f τ , x ( τ ) d τ k = 0 n w ( t k ) f t k , x ( t k ) t k t k + 1 ( t n + 1 τ ) β 1 d τ = h β β k = 0 n w ( t k ) f t k , x ( t k ) A n , k β ,
where
A n , k β = ( n k + 1 ) β ( n k ) β .
Therefore, we obtain the following numerical scheme:
x n + 1 = x 0 w ( 0 ) w ( t n ) + 1 α N ( α ) f t n , x n + α h β N ( α ) Γ ( β + 1 ) w ( t n ) k = 0 n w ( t k ) f t k , x k A n , k β .
Remark 5.
If α = β = 1 and w ( t ) = 1 , then (28) becomes
x n + 1 = x n + h f t n , x n .
Hence, the classical Euler numerical scheme for ODEs is recovered. Indeed,
x n + 1 = x 0 + h N ( 1 ) Γ ( 2 ) k = 0 n f t k , x k A n , k 1 = x 0 + h k = 0 n f t k , x k = x 0 + h k = 0 n 1 f t k , x k + h f t n , x n = x n + h f t n , x n .

5. Application to Biology

In this section, we apply our main analytical and numerical results to the biological system describing the dynamics of an epidemic disease:
D 0 , 1 α , β S ( t ) = A μ S ( t ) κ S ( t ) I ( t ) , D 0 , 1 α , β E ( t ) = κ S ( t ) I ( t ) ( μ + ε ) E ( t ) , D 0 , 1 α , β I ( t ) = ε E ( t ) ( μ + r ) I ( t ) , D 0 , 1 α , β R ( t ) = r I ( t ) μ R ( t ) ,
where S ( t ) , I ( t ) , and R ( t ) are the fractions of susceptible, exposed, infectious, and recovered individuals at time t, respectively. The biological meanings of the parameters are presented in Table 1.
Since the first three equations of (30) do not depend on the last one, System (30) can be reduced to the following model:
D 0 , 1 α , β S ( t ) = A μ S ( t ) κ S ( t ) I ( t ) , D 0 , 1 α , β E ( t ) = κ S ( t ) I ( t ) ( μ + ε ) E ( t ) , D 0 , 1 α , β I ( t ) = ε E ( t ) ( μ + r ) I ( t ) .
In fact, when the variable I ( t ) is determined by (31), then we easily obtain R ( t ) from the last equation of (30).
Clearly, Model (31) has a unique disease-free equilibrium Q 0 = ( S 0 , 0 , 0 ) , where S 0 = A μ . Furthermore, the basic reproduction number of (31) is given by
R 0 = ε κ S 0 ( μ + ε ) ( μ + r ) .
By a simple computation, Model (31) has a unique endemic Q * = ( S * , E * , I * ) , where S * = A μ R 0 , E * = μ ( μ + r ) ( R 0 1 ) ε κ and I * = μ ( R 0 1 ) κ .
Let Ω = { ( S , E , I ) I R + 3 : S S 0 } . For R 0 < 1 , construct a Lyapunov function as follows:
V ( S , E , I ) = ρ ( S 0 S ) + E + μ + ε 2 ε I ,
where ρ = 2 R 0 2 R 0 . We have
D 0 , 1 α , β V ( S , E , I ) = ρ D 0 , 1 α , β S + D 0 , 1 α , β E + μ + ε 2 ε D 0 , 1 α , β I ρ μ ( S 0 S ) μ + ε 2 E ( μ + ε ) ( μ + r ) 2 ε ( 1 R 0 ) I .
Thus,
D 0 , 1 α , β V ( S , E , I ) σ 1 V ( S , E , I ) ,
where σ 1 = min { μ , μ + ε 2 , ( μ + r ) ( 1 R 0 ) } .
Let X = ( S , E , I ) I R 3 with the norm X = | S | + | E | + | I | . Hence,
σ 2 X Q 0 V ( S , E , I ) σ 3 X Q 0 ,
where σ 2 = min { ρ , 1 , μ + ε 2 ε } and σ 3 = max { ρ , 1 , μ + ε 2 ε } . By applying Theorem 4, we deduce that the disease-free equilibrium Q 0 of (31) is Mittag–Leffler stable in Ω when R 0 < 1 .
For R 0 > 1 , consider the following Lyapunov function:
L ( S , E , I ) = S * Φ ( S S * ) + E * Φ ( E E * ) + μ + ε ε I * Φ ( I I * ) ,
where Φ ( z ) = z 1 ln z , for z > 0 . It is obvious that Φ ( z ) attains its global minimum at z = 1 and Φ ( 1 ) = 0 . Then, Φ ( z ) 0 for all z > 0 . Hence, L ( S , E , I ) 0 with L ( S * , E * , I * ) = 0 .
By applying Corollary 2 of [34], we obtain
D 0 , 1 α , β L ( S , E , I ) 1 S * S D 0 , 1 α , β S + 1 E * E D 0 , 1 α , β E + μ + ε ε 1 I * I D 0 , 1 α , β I .
Using A = μ S * + ( μ + ε ) E * , κ S * I * = ( μ + ε ) E * and μ + r ε = E * I * , we obtain
D 0 , 1 α , β L ( S , E , I ) μ S S S * 2 + ( μ + ε ) E * 3 S * S E I * E * I S E * I S * E I * = μ S S S * 2 ( μ + ε ) E * Φ S * S + Φ E I * E * I + Φ S E * I S * E I * .
Hence, D 0 , 1 α , β L ( S , E , I ) 0 when R 0 > 1 . It follows from Theorem 5 that the endemic equilibrium Q * of (31) is stable when R 0 > 1 .
In the absence of disease, System (31) reduces to the following linear system:
D 0 , 1 α , β S ( t ) = A μ S ( t ) .
From Lemma 2 of [28], the exact solution of (35) is given by
S ( t ) = A μ + N ( α ) a α S ( 0 ) A μ E β α μ a α t β ,
a α = N ( α ) + μ ( 1 α ) . Now, we apply the numerical scheme presented in (28) in order to approximate the solution of (35). For all numerical simulations, we chose A = 0.01 , μ = 0.01 , and the normalization function as follows:
N ( α ) = 1 α + α Γ ( α ) .
The comparison between the exact and numerical (approximate) solutions of (35) is displayed in Figure 1 for different values of α , β , and h.
In the presence of disease, System (31) cannot be solved analytically. Based on our numerical method, we approximate the solution of (31). Therefore, we chose κ = 0.4 , ε = 0.085 , r = 0.095 , and β = 0.8 . By computation, we have R 0 = 3.4085 > 1 . In this case, the solution of (31) converges to the endemic equilibrium Q * ( 0.2934 , 0.0744 , 0.0602 ) , which biologically means that the disease persists in the population. Figure 2 illustrates this observation for different values of α .

6. Conclusions

In this paper, we first investigated the qualitative properties of solutions of FDEs with the new generalized Hattaf fractional derivative, which includes several forms of fractional derivatives with non-singular kernels such as the Caputo–Fabrizio and Atangana–Baleanu fractional derivatives. In addition, we proposed a new numerical method to approximate the solutions of such types of FDEs. The obtained results extend and improve many results existing in the literature concerning the fractional comparison principle, stability, asymptotic stability, as well as Mittag–Leffler stability. Furthermore, the proposed numerical method includes the classical Euler numerical scheme, and it was applied to a nonlinear system describing the dynamics of an epidemic disease, such as COVID-19.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No data available.

Acknowledgments

The author would like to thank the anonymous reviewers for their constructive comments and suggestions, which have improved the quality of the manuscript.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Hattaf, K.; Yousfi, N. Global stability for fractional diffusion equations in biological systems. Complexity 2020, 2020, 5476842. [Google Scholar] [CrossRef]
  2. Magin, R.L. Fractional calculus models of complex dynamics in biological tissues. Comput. Math. Appl. 2010, 59, 1586–1593. [Google Scholar] [CrossRef] [Green Version]
  3. Cheneke, K.R.; Rao, K.P.; Edessa, G.K. Application of a new generalized fractional derivative and rank of control measures on Cholera transmission dynamics. Int. J. Math. Math. Sci. 2021, 2021, 2104051. [Google Scholar] [CrossRef]
  4. Qu, H.; Rahman, M.U.; Ahmad, S.; Riazd, M.B.; Ibrahim, M.; Saeed, T. Investigation of fractional order bacteria dependent disease with the effects of different contact rates. Chaos Solitons Fractals 2022, 159, 112169. [Google Scholar] [CrossRef]
  5. Zhang, L.; Rahman, M.U.; Ahmad, S.; Riaz, M.B.; Jarad, F. Dynamics of fractional order delay model of coronavirus disease. Aims Math. 2021, 7, 4211–4232. [Google Scholar] [CrossRef]
  6. Naji, F.A.; Al-Sharaa, I. Controllability of impulsive fractional nonlinear control system with Mittag–Leffler kernel in Banach space. Int. J. Nonlinear Anal. Appl. 2022, 13, 3257–3280. [Google Scholar]
  7. Meral, F.C.; Royston, T.J.; Magin, R.L. Fractional calculus in viscoelasticity: An experimental study. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 939–945. [Google Scholar] [CrossRef]
  8. Pinar, Z. On the explicit solutions of fractional Bagley-Torvik equation arises in engineering. Int. J. Optim. Control Theor. Appl. 2019, 9, 52–58. [Google Scholar] [CrossRef] [Green Version]
  9. Magin, R.L. Fractional Calculus in Bioengineering; Begell House: Danbury, CT, USA, 2006. [Google Scholar]
  10. Cao, X.; Lin, Y.-H.; Liu, H. Simultaneously recovering potentials and embedded obstacles for anisotropic fractional Schrödinger operators. Inverse Probl. Imaging 2019, 13, 197–210. [Google Scholar] [CrossRef] [Green Version]
  11. Cao, X.; Liu, H. Determining a fractional Helmholtz equation with unknown source and scattering potential. Commun. Math. Sci. 2019, 17, 1861–1876. [Google Scholar] [CrossRef]
  12. Lai, R.-Y.; Lin, Y.-H. Inverse problems for fractional semilinear elliptic equations. Nonlinear Anal. 2022, 216, 112699. [Google Scholar] [CrossRef]
  13. Srivastava, H.M.; Sahoo, S.K.; Mohammed, P.O.; Baleanu, D.; Kodamasingh, B. Hermite–Hadamard Type Inequalities for Interval-Valued Preinvex Functions via Fractional Integral Operators. Int. J. Comput. Intell. Syst. 2022, 15, 8. [Google Scholar] [CrossRef]
  14. Yang, S.; Liu, Y.; Liu, H.; Wang, C. Numerical methods for semilinear fractional diffusion equations with time delay. Adv. Appl. Math. Mech. 2022, 14, 56–78. [Google Scholar]
  15. Butt, A.I.K.; Ahmad, W.; Rafiq, M.; Baleanu, D. Numerical analysis of Atangana–Baleanu fractional model to understand the propagation of a novel corona virus pandemic. Alex. Eng. J. 2022, 61, 7007–7027. [Google Scholar] [CrossRef]
  16. Rahman, M.U.; Arfan, M.; Deebani, W.; Kumam, P.; Shah, Z. Analysis of time-fractional Kawahara equation under Mittag–Leffler Power Law. Fractals 2022, 30, 2240021. [Google Scholar] [CrossRef]
  17. Odabasi, M.; Pinar, Z.; Kocak, H. Analytical solutions of some nonlinear fractional-order differential equations by different methods. Math. Methods Appl. Sci. 2021, 44, 7526–7537. [Google Scholar] [CrossRef]
  18. Ashyralyev, A.; Dal, F.; Pinar, Z. A note on the fractional hyperbolic differential and difference equations. Appl. Math. Comput. 2011, 217, 4654–4664. [Google Scholar] [CrossRef]
  19. Ashyralyev, A.; Dal, F.; Pinar, Z. On the numerical solution of fractional hyperbolic partial differential equations. Math. Probl. Eng. 2009, 2009, 730465. [Google Scholar] [CrossRef]
  20. Li, Y.; Chen, Y.Q.; Podlubny, I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag–Leffler stability. Comput. Math. Appl. 2010, 59, 1810–1821. [Google Scholar] [CrossRef] [Green Version]
  21. Podlubny, I. Mathematics in Science and Engineering. In Fractional Differential Equations; Academic Press: Cambridge, MA, USA, 1999; Volume 198. [Google Scholar]
  22. Delavari, H.; Baleanu, D.; Sadati, J. Stability analysis of Caputo fractional-order nonlinear systems revisited. Nonlinear Dyn. 2012, 67, 2433–2439. [Google Scholar] [CrossRef]
  23. Rao, M.R. Ordinary Differential Equations; East-West Press: Minneapolis, MN, USA, 1980. [Google Scholar]
  24. Slotine, J.J.E.; Li, W. Applied Nonlinear Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1991. [Google Scholar]
  25. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  26. Wang, G.; Pei, K.; Chen, Y.Q. Stability analysis of nonlinear Hadamard fractional differential system. J. Frankl. Inst. 2019, 356, 6538–6546. [Google Scholar] [CrossRef]
  27. Brzdek, J.; Eghbali, N.; Kalvandi, V. On Ulam stability of a generalized delayed differential equation of fractional order. Results Math. 2022, 77, 26. [Google Scholar] [CrossRef]
  28. Hattaf, K. Stability of fractional differential equations with new generalized Hattaf fractional derivative. Math. Probl. Eng. 2021, 2021, 8608447. [Google Scholar] [CrossRef]
  29. Hattaf, K. A new generalized definition of fractional derivative with non-singular kernel. Computation 2020, 8, 49. [Google Scholar] [CrossRef]
  30. Caputo, A.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  31. Atangana, A.; Baleanu, D. New fractional derivatives with non-local and non-singular kernel: Theory and application to heat transfer model. Therm. Sci. 2016, 20, 763–769. [Google Scholar] [CrossRef] [Green Version]
  32. Al-Refai, M. On weighted Atangana–Baleanu fractional operators. Adv. Differ. Equ. 2020, 2020, 3. [Google Scholar] [CrossRef] [Green Version]
  33. Hattaf, K.; Mohsen, A.A.; Al-Husseiny, H.F. Gronwall inequality and existence of solutions for differential equations with generalized Hattaf fractional derivative. Math. Comput. Sci. 2022, 27, 18–27. [Google Scholar] [CrossRef]
  34. Hattaf, K. On some properties of the new generalized fractional derivative with non-singular kernel. Math. Probl. Eng. 2021, 2021, 1580396. [Google Scholar] [CrossRef]
  35. Djida, J.D.; Atangana, A.; Area, I. Numerical Computation of a Fractional Derivative with Non-Local and Non-Singular Kernel. Math. Model. Nat. Phenom. 2017, 12, 4–13. [Google Scholar] [CrossRef] [Green Version]
  36. Baleanu, D.; Fernandez, A. On some new properties of fractional derivatives with Mittag–Leffler kernel. Commun. Nonlinear Sci. Numer. Simul. 2018, 59, 444–462. [Google Scholar] [CrossRef] [Green Version]
  37. Zhang, F.; Li, C.; Chen, Y.Q. Asymptotical stability of nonlinear fractional differential system with Caputo derivative. Int. J. Differ. Equ. 2011, 2011, 635165. [Google Scholar] [CrossRef]
Figure 1. The exact and numerical solutions of (35) for different values of α , β , and h.
Figure 1. The exact and numerical solutions of (35) for different values of α , β , and h.
Computation 10 00097 g001
Figure 2. The numerical solution of (31) for different values of α .
Figure 2. The numerical solution of (31) for different values of α .
Computation 10 00097 g002
Table 1. Biological meanings of the parameters of model (30).
Table 1. Biological meanings of the parameters of model (30).
ParameterBiological Meaning
ANatality or recruitment rate
μ Natural death rate
κ Transmission rate of disease
ε Transfer rate from class E to class I
rRecovery rate of the infectious individuals
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hattaf, K. On the Stability and Numerical Scheme of Fractional Differential Equations with Application to Biology. Computation 2022, 10, 97. https://0-doi-org.brum.beds.ac.uk/10.3390/computation10060097

AMA Style

Hattaf K. On the Stability and Numerical Scheme of Fractional Differential Equations with Application to Biology. Computation. 2022; 10(6):97. https://0-doi-org.brum.beds.ac.uk/10.3390/computation10060097

Chicago/Turabian Style

Hattaf, Khalid. 2022. "On the Stability and Numerical Scheme of Fractional Differential Equations with Application to Biology" Computation 10, no. 6: 97. https://0-doi-org.brum.beds.ac.uk/10.3390/computation10060097

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop