key: cord-1042688-dmtb3icf
authors: Ahmed, E.; Elgazzar, A.S.
title: On fractional order differential equations model for nonlocal epidemics
date: 2007-06-15
journal: Physica A
DOI: 10.1016/j.physa.2007.01.010
sha: 9b4aa91c25a7dac9f3eb55b10b7ef5bfde8e6df0
doc_id: 1042688
cord_uid: dmtb3icf

A fractional order model for nonlocal epidemics is given. Stability of fractional order equations is studied. The results are expected to be relevant to foot-and-mouth disease, SARS and avian flu.

Recently three major epidemic diseases have occurred namely foot-and-mouth disease, severe acute respiratory syndrome (SARS) and avian (bird's) flu. Hopefully this will increase the awareness of modeling infectious diseases spreading that is an important topic in mathematical biology [1] . There are different approaches to this topic, e.g. ordinary differential equations, difference equations, partial differential equations and coupled map lattice. Here we use fractional order differential equations (FOD). The reason is that FOD are naturally related to systems with memory which exists in most biological systems. Also they are closely related to fractals which are abundant in biological systems. Consider the following evolution equation [2] :

kðt À t 0 Þf ðt 0 Þ dt 0 .

If the system has no memory then k(tÀt 0 ) ¼ d(tÀt 0 ), and one gets f(t) ¼ f 0 exp(Àl 2 t). If the system has an ideal memory, then 

where G(a) is the Gamma function. This system has the following solution:

where E a (z) is the Mittag-Leffler function given by

Following a similar procedure to study a random process with memory, one obtains the following fractional evolution equation:

where P(x, t) is a measure of the probability to find a particle at time t at position x. We expect that the above result will be relevant to many complex adaptive systems and to systems where fractal structures are relevant since it is argued that there is a strong relevance between fractals and fractional differentiation [3] . For the case of fractional diffusion equation the results are q aþ1 Pðx; tÞ qt aþ1 ¼ D q 2 Pðx; tÞ qx 2 ; Pðx; 0Þ ¼ dðxÞ;

ðÀ1Þ n z n n!GðÀbn þ 1 À bÞ .

For the case of no memory a ¼ 0 ) M(z; 1/2) ¼ exp(Àz 2 /4). The paper is organized as follows: in Section 2, we study the stability of FOD. The stability conditions are derived and several examples are given. The stability conditions for some fractional order differential coupled map lattices are concluded in Section 3. Applications to nonlocal epidemics is introduced in Section 4. The stability conditions for the disease-free state are discussed. Some conclusions are presented in Section 5.

2. Stability of fractional order differential equations [4] [5] [6] Consider the following system: 

where the fractional derivative in Eq. (1) is in the sense of Caputo. The equilibrium solutions are defined by f(x eq , y eq ) ¼ 0, g(x eq , y eq ) ¼ 0 and it is locally asymptotically stable if all the eigenvalues l of the Jacobian

" # evaluated at the equilibrium point satisfies the following condition [4, 5] :

(2)

The condition in Eq. (2) poses an interesting question namely: What are the conditions that all the roots of the polynomial equation

satisfy Eq. (2) where all the coefficients in Eq. (3) are real? For a ¼ 1, the solution is the Routh-Hurwitz conditions [7] a 1 40; a 1 1 a 3 a 2 40;

For aA[0, 1), these conditions are sufficient but not necessary. Since most biologically interesting systems are 1,2 and 3-dimensions, we will study the problem (3) for n ¼ 1, 2, 3.

is the determinant of the corresponding Sylvester (n+l)(n+l) matrix. The Sylvester matrix is formed by filling the matrix beginning with the upper left corner with the coefficients of f (x), then shifting down one row and one column to the right and filling in the coefficients starting there until they hit the right side. The process is then repeated for the coefficients of g (x).

Using the results of Ref. [3] , if D( f )40(o0) then there is an even (odd) number of pairs of complex roots for the equation f(x) ¼ 0. For n ¼ 3, this implies that D (f)40, and all the roots are real while D (f)o0 implies that there is only one real root and one complex root and its complex conjugate. For n ¼ 3, we have

Proposition 1.

(i) For n ¼ 1, the condition for (3) is a 1 40.

(ii) For n ¼ 2, the conditions for (3) are either Routh-Hurwitz conditions or (viii) For n41, the necessary and sufficient condition for (3) is

Proof. Case (i) is obvious. For (ii), the two roots are l AE ¼ Àa 1 AE ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi

If both roots are real or complex conjugates with negative real parts then condition (3) is equivalent to the Routh-Hurwitz case. If the two roots are complex conjugate with positive real parts, then the two roots become l AE ¼ Àa 1 AE i ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi 4a 2 À ða 1 Þ 2 q =2 and one gets Eq. (5).

To prove (iii) not that if n ¼ 3, D (p)40, then all the roots of p(l) ¼ 0 are real hence Routh-Hurwitz conditions are both necessary and sufficient for (3).

To prove (iv) not that if n ¼ 3, D (P)40, then the roots of P(l) ¼ 0 are one real and a complex conjugate pair thus

and a 1 40 ) b42b, a 2 40 ) b 2 sec 2 y42bb44b 2 ) y4p/3, where y ¼ |arg(l)|. The second part is proved similarly.

To prove (v) if n ¼ 3, D(P)40, then Eq. (7) is

The last equality is not valid if both a 1 40, a 2 40.

To prove (vi) use the fact that for general n

To prove (vii) use Eq. (8) and P(l) ¼ P(Àl)8l ) P(l) contains only even power of

And all the roots x should be negative.

To prove (viii) not that if P(z) has no roots in the region |arg(l)|oap/2, hence the function 1/P(z) will be analytic in this region. Using Cauchy theorem H C f ðzÞ dz ¼ 0 for all f(z) analytic within and on the curve C, and that P(z) is polynomial of degree 41 this completes the proof. & One can give explicit examples where the Routh-Hurwitz conditions are not valid yet Eq. (3) is satisfied for explicit a, e.g. a ¼ 1 2 . The first example is

For n ¼ 3, one can solve the polynomial equation explicitly and apply Eq. (3) to the solutions. The solution method is:

Define y ¼ lÀa 1 /2, then the polynomial becomes

choose the roots such that rs ¼ Àp/3. For all n43, if D 1 ; D 2 ; . . . ; D n are Routh-Hurwitz determinants

. . . , then the conditions

. . . ; n À 2; a n 40; D nÀ1 ¼ 0, (9) are sufficient conditions that Eq. (3) is valid for all aA[0, 1).

This conjecture can be proved for n ¼ 4. The case n ¼ 3 is proved in Proposition 1. Now we apply Proposition 1 to derive the value of a at which chaos or instability disappears at some models. Fractional order Lotka-Volterra (FLV) predator-rey model is given by

where b, g are positive constants. The equilibrium solutions of FLV are (0, 0) and (g, b) . The eigenvalues corresponding to (0, 0) are b, Àg respectively hence (0, 0) is unstable for all aA[0, 1]. The eigenvalues corresponding to (g, b) are AEi ffiffiffiffiffi bg p hence |arg(l)| ¼ p/2, then it is locally asymptotically stable for all aA[0, 1). This agrees with the numerical results of Ref. [5] .

Fractional order Chen (FOC) model is given by [7] D a x ¼ aðy À xÞ; D a y ¼ ðc À aÞx À xz þ cy; D a z ¼ xy À bz,

The equilibrium solutions for FOC are (0, 0, 0) and ðx n ; x n ; x n2 =bÞ, where x n ¼ AE ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi bð2c À aÞ p . The eigenvalues of (0, 0, 0) are l ¼ Àb, l 2 +(aÀc)lÀa(2cÀa) ¼ 0. Using Proposition 1 part (vi), then (0, 0, 0) is unstable for all aA[0, 1]. The eigenvalues of the internal equilibrium solution is l 3 +(aÀc+b)l 2 +bcl+2ab(2cÀa) ¼ 0. Since all the coefficients of this polynomial are positive and it is easy to check that D(P)o0, using Proposition 1 part (iv) then the internal solution of FOC will be stable for ao 2 3 which is in excellent agreement with the numerical result obtained by and Chen [8] .

A similar result is obtained for fractional order Lorenz system

where r, b, c are positive constants and r41. The fractional order Rossler system (FOR) is defined by [9] D a x ¼ Àðy þ zÞ; D a y ¼ x þ ay; D a z ¼ 0:2 þ xz À 10z; a ¼ 0:6 There are one equilibrium z ¼ Ày, x ¼ Àay, y ¼ ½À10 þ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 100 À 0:8a p =2a, and the characteristic polynomial is l 3 þ ðay À a þ 10Þl 2 þ ð1 À ya 2 À 10a À yÞl þ 10 þ 2ay ¼ 0.

It is direct to see that a 2 o0, a 1 40 hence by solving the characteristic equation directly one gets the following roots: l ¼ À9:986; l ¼ 0:3 À 0:954i; l ¼ 0:3 þ 0:954i.

This implies that the equilibrium solution is locally asymptotically stable if ao0.8 which is in excellent agreement with the result obtained numerically in Ref. [9] : aA(0.7, 0.8).

Considering the fractional order Chua model D a x ¼ aðy À xÞ À gðxÞ; D a y ¼ s½Àaðy À xÞ þ z, D a z ¼ Àcðy þ rzÞ,

½10ðjxj À 10Þ þ 9b À 1sgnðxÞ if jxj410; 

which implies that this equilibrium is locally asymptotically stable for aA[0,1). The eigenvalues for a ¼ 0.

188i which implies that this equilibrium is locally asymptotically stable for ao0.9.

3. Stability conditions for some fractional order differential coupled map lattices Spatial effects are important in many biological systems. Thus generalizing fractional order systems (FOS) to include them is important. The standard approach is fractional order partial differential equations. However since most biologically interesting systems are nonlinear [10] , one gets fractional order nonlinear partial differential equations whose existence and uniqueness has not been established yet. Therefore, we use coupled map lattices (CML) [11] as an alternative approach to include spatial effects in FOS. Consider the 1-system CML

where D is a positive constant. The homogeneous equilibrium solution of Eq. (10) satisfies f(u eq ) ¼ 0 and it is stable if all the eigenvalues of the circulant matrix [12] ð1 À DÞf 0 ðD=2Þf 0 0 . . . satisfies Eq. (2). Since the eigenvalues of the circulant matrices [12] are known, the local asymptotic stability conditions become arg ð1 À DÞf 0 ðu eq Þ þ Df 0 ðu eq Þ cos 2pk n 4 ap 2 , k ¼ 0; 1; . . . ; n À 1 ð11Þ

For n ¼ 3, the characteristic polynomial of A is

where a 1 ¼ À½l 1 À g 1 þ l 2 À g 2 þ l 3 À g 3 , a 2 ¼ ðl 1 À g 1 Þðl 2 À g 2 Þ À m 1 m 2 þ ðl 1 À g 1 Þðl 3 À g 3 Þ À m 1 m 3 þ ðl 3 À g 3 Þðl 2 À g 2 Þ À m 3 m 2 and a 3 ¼ À½ðl 1 À g 1 Þðl 2 À g 2 Þðl 3 À g 3 Þ þ 2m 1 m 2 m 3 À ðl 1 À g 1 Þm 2 m 3 À ðl 3 À g 3 Þm 2 m 1 À ðl 2 À g 2 Þm 1 m 3 .

Using part (vi) in Proposition 1, the disease free state is unstable if a 3 o0.

Finally vaccination may affect the parameters l i but not m i . Also it is important to know whether the vaccine works once administered or it takes time till it becomes effective.

A fractional order evolution equation describing a random process with memory is concluded. This will be relevant to many complex adaptive systems. The stability conditions for both fractional order differential equations and fractional order differential coupled map lattices are derived. A fractional order form of Lajmanovich-Yorke model is introduced. This model emphasis the nonlocal interaction beside the local interaction in the epidemic spreading. Both types of interaction are widely observed in the recent FMD, SARS and bird's flu outbreaks.

Memory effects and macroscopic manifestation of randomness

Fractional calculus and the evolution of fractal phenomena

Stability results for fractional differential equations with applications to control processing

Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models

Higher algebra

Chaos in the fractional order Chen model and its control

Chaos and hyperchaos in fractional order Rossler equations

Mathematical models in biology

Theory and Applications of Coupled Map Lattices

A deterministic model for gonorrhea in a non-homogeneous population

Stochastic and deterministic models for SIS epidemics among a population partitioned into households

An epidemic model on small-world networks and ring vaccination

One of us (E. Ahmed) wishes to thank Prof. H. M. Yehia for discussions and help.

Generalizing to 2-system CMLj ¼ 1, 2, y, n. The homogeneous equilibrium solutions are given by f(x eq , y eq ) ¼ 0, g(x eq , y eq ) ¼ 0 and is locally asymptotically stable if all the eigenvalues of the following matrix B satisfies Eq. (2):where the Jacobian matrix is evaluated at the equilibrium point.

Now we present a nonlocal interacting epidemic model. The nonlocal interactions in the epidemic spreading is widely observed in many outbreaks, especially the recent ones of FMD, SARS and bird's flu. The Lajmanovich-Yorke model [13, 14] is a special case of coupled map lattices. It is a globally interacting model that emphasizes the diffusion of a disease from other infected farms beside its spread within the same farm. We have generalized it to the inhomogeneous case [15] and introduced it as an approximation for the epidemic spreading on small-world networks [15] . Consider n patches, where each one contains a certain number of individuals (say animals). In general these patches are not identical. Infection spreads from infected animals within the patch and due to those diffusing from other patches. The fractional order form of the model is given by D a y i ¼ l i y i ð1 À y i Þ þ m i ð1 À y i Þ X jai y j À g i y i ; i ¼ 1; 2; . . . ; n,where y i is the number of infected individuals in the ith patch. The first term represents the infection within the patch with a rate l i . The second term represents the effect of other patches both nearby and far away at a rate m i . The recovery rate is represented by g i . Since the effect of other patches (second term in Eq. (14)) is significantly smaller than the first one, we expect that m i 5l i . The disease is eradicated if y i ¼ 0, 8 i ¼ 1, 2, y, n. So the stability of this solution is studied. Using the same procedures in the previous section, the stability of the zero solution (disease free state) is that all the eigenvalues x of the matrix A satisfy |arg(x)|4ap/2 where A ¼:: :: :: :: :::: :: :: :: :: m n m n m n . . . l n À g n ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðl 1 À g 1 þ l 2 À g 2 Þ 2 À 4½ðl 1 À g 1 Þðl 1 À g 1 Þ À m 1 m 2 q l 1 À g 1 þ l 2 À g 2 0 @ 1 A 4 ap 2 .