Next Article in Journal
Primality, Fractality, and Image Analysis
Next Article in Special Issue
A Review of the Application of Information Theory to Clinical Diagnostic Testing
Previous Article in Journal
Category Theory for Autonomous and Networked Dynamical Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of Ebola Disease in the Framework of Different Fractional Derivatives

by
Khan Muhammad Altaf
1,† and
Abdon Atangana
2,*,†
1
Department of Mathematics, City University of Science and Information Technology, Peshawar 25000, Pakistan
2
Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein 9300, South Africa
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Entropy 2019, 21(3), 303; https://doi.org/10.3390/e21030303
Submission received: 18 February 2019 / Revised: 28 February 2019 / Accepted: 7 March 2019 / Published: 21 March 2019
(This article belongs to the Special Issue Applications of Information Theory to Epidemiology)

Abstract

:
In recent years the world has witnessed the arrival of deadly infectious diseases that have taken many lives across the globe. To fight back these diseases or control their spread, mankind relies on modeling and medicine to control, cure, and predict the behavior of such problems. In the case of Ebola, we observe spread that follows a fading memory process and also shows crossover behavior. Therefore, to capture this kind of spread one needs to use differential operators that posses crossover properties and fading memory. We analyze the Ebola disease model by considering three differential operators, that is the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu operators. We present brief detail and some mathematical analysis for each operator applied to the Ebola model. We present a numerical approach for the solution of each operator. Further, numerical results for each operator with various values of the fractional order parameter α are presented. A comparison of the suggested operators on the Ebola disease model in the form of graphics is presented. We show that by decreasing the value of the fractional order parameter α , the number of individuals infected by Ebola decreases efficiently and conclude that for disease elimination, the Atangana–Baleanu operator is more useful than the other two.

1. Introduction

Ebola caused many deaths in Western Africa, especially in the outbreak of 2014. It includes more than 16 thousand laboratory cases with 70% death cases, which is regarded the deadliest outbreak in history since 1976 with 20 Ebola threats. It is evident that in each outbreak, the first case of infection occurred due to contact with infected animals such as monkeys, fruit bats, etc., which shows the spread of the virus through indirect contact [1]. It is documented in [2] that some percentage of the Ebola-Zaire type survived after two weeks on glass at 4 °C and (10%) on plastic, and on surfaces (3%). Moreover, 0.1% to 1 % of the Ebola virus particle can remain up to 50 days at 4 °C [3]. The survival of the Ebola virus in the environment due to poor sanitary and hygienic conditions considerably become another source of Ebola infection in Africa. In Africa, regions were affected greatly by the Ebola virus outbreak due to their inhabitants being involved in hunting food, being close to the rain-forest, and harvesting forest fruits for food [4,5].
The Ebola disease outbreaks and their transmission have been documented in many articles (see [6,7,8,9,10] and the references therein) and the main focus was to study the human population and the direct transmission. Some models of the type SI, SIR, SEIR, and other types also considered the dynamics of the Ebola disease outbreaks [9,11,12]. Recently, in [13] studied an Ebola virus disease through a simple mathematical model of the type SIR with the inclusion of environment effect. Due to the fact that Ebola virus survives in the environment, this warrants that future epidemics can occur. Thus, the inclusion of the environment effect in Ebola disease spread should be studied more and some preventive and other measures should be used to protect people further from this deadly infection. Therefore, based on the model presented in [13], we aim to study the Ebola disease model in the framework of the fractional calculus. The reason for the use of the fractional calculus in Ebola disease is that it has many advantages. Some of them are the heredity and memory effects, the parameter estimations are better, the crossover behavior of the model, and effective strategies for the case of arbitrary order. Some other works used it to study the dynamics of complex networks [14,15,16]. In [14] the authors studied the dynamics of information and the uses in complex networks. Coupling dynamics of an epidemic spreading with information diffusion is analyzed in [15]. The events that determine spreading dynamics and the information transmission through internal and external influences are considered in [16].
Fractional calculus and its applications to real life problems is found extensively in the literature, for example [17,18,19,20,21]. In all these mentioned papers the focus is to eliminate the infection from the community and it is proven that the fractional models have the ability to model such epidemic disease efficiently and provide reasonable results for the case of non-integer. It is shown that the fractional models are useful for the data fitting [22]. The results suggest in [22] that fractional models are efficient to study disease dynamics well. Therefore, motivated with the above applications, we aim to study an Ebola disease model in the fractional order. We consider three different fractional operators, that is, the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives. According to the authors’ knowledge no one has applied the three operators to an epidemic model. So, this work is a useful study to analyze Ebola disease with different fractional operators. The rest of the work on Ebola disease is categorized as follows: The fractional background material are shown in Section 2. A mathematical model on Ebola disease is presented in Section 3 with basic mathematical results. In Section 4, a mathematical model in the frame of the Caputo derivative and their numerical results, the Caputo–Fabrizio derivative is used to formulate the model and their relevant results are presented, and we further consider the Atangana–Baleanu model for Ebola disease and discuss its existence and uniqueness and a useful numerical scheme for their solution, and lastly in this section, the comparison results for these operators with various fractional order parameters are shown. The Ebola disease models and their fractional results are summarized in Section 5.

2. Fundamental Concepts

Here, we recall the fundamental concepts regarding the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivative.
Definition 1.
For a function f : R + R , then the fractional integral of order α > 0 is given by
I t α ( f ( t ) ) = 1 Γ ( α ) 0 t ( t z ) α 1 f ( z ) d z .
where Γ shows the Gamma function and α is the fractional order parameter.
Definition 2.
For a function f C n , then the Caputo derivative with order α is defined as
C D t α ( f ( t ) ) = I n α D n f ( t ) = 1 Γ ( n α ) 0 t f n ( z ) ( t z ) α + n 1 d z ,
that is defined for the absolute continuous functions and n 1 < α < n N . Obviously, C D t α ( f ( t ) ) tends to f ( t ) as α 1 .
Definition 3.
[23]. Let z H 1 ( a , b ) , with b > a , and 0 α 1 , then the Caputo–Fabrizio derivative can be written as
D t α ( z ( t ) ) = K ( α ) 1 α a t z ( x ) exp α t x 1 α d x ,
the normalized function is shown by K ( α ) and it holds K ( 0 ) = K ( 1 ) = 1 . Consider the case for which z H 1 ( a , b ) then, we have the following:
D t α ( z ( t ) ) = α K ( α ) 1 α a t ( z ( t ) z ( x ) ) exp α t x 1 α d x .
Remark 1.
[24]. Let ν = 1 α α [ 0 , ) , α = 1 1 + ν [ 0 , 1 ] , then equation given by (2) can be expressed is as follows,
D t ν ( z ( t ) ) = K ( ν ) ν a t z ( x ) exp t x ν d x , K ( 0 ) = K ( ) = 1 .
Further,
lim ν 0 1 ν exp t x ν = φ ( x t ) .
Definition 4.
Consider α ( 0 , 1 ) , for a function z ( x ) then we can write the integral of fractional order α is as follows,
I t α ( z ( t ) ) = 2 ( 1 α ) ( 2 α ) K ( α ) g ( t ) + 2 α ( 2 α ) K ( α ) 0 t z ( s ) d s , t 0 .
Remark 2.
In Equation (4), the remainder of the Caputo type non-integer order integral of the function with order α ( 0 , 1 ) is a mean into z with integral of order 1. Thus, it requires,
2 2 K ( α ) α K ( α ) = 1 ,
implies that K ( α ) = 2 2 α , α ( 0 , 1 ) . Based on Equation (6), a new Caputo derivative is suggested with α ( 0 , 1 ) and is given by
D t α ( z ( t ) ) = 1 1 α 0 t z ( x ) exp α t x 1 α d x .
In the following we present the new derivative known as the Atangana–Baleanu derivatives having non-singular and non-local kernel [25].
Definition 5.
Let f H 1 ( a , b ) , b > a , α [ 0 , 1 ] , then in the Caputo sense the Atangana–Baleanu derivative is defined as:
a A B C D t α f ( t ) = K ( α ) 1 α a t f ( z ) E α α ( t z ) α 1 α d z .
Definition 6.
The fractional integral associated with the Atangana–Beleanu derivative is given by:
a A B C I t α f ( t ) = 1 α K ( α ) f ( t ) + α K ( α ) Γ ( α ) a t f ( z ) ( t z ) α 1 d z .
when the fractional order turns to zero, we can obtain the original function.
Theorem 1.
Consider the function f C [ a , b ] , then the following holds [25]:
a A B C D t α ( f ( t ) ) < K ( α ) 1 α f ( t ) , w h e r e f ( t ) = m a x a t b | f ( t ) | .
Further, for the newly derivative the Lipschitz condition can be easily satisfied [25]:
a A B C D t α f 1 ( t ) a A B C D t α f 2 ( t ) < ϖ 1 f 1 ( t ) f 2 ( t ) .
Theorem 2.
A given fractional differential equation:
a A B C D t α f ( t ) = s ( t ) ,
has the unique solution given by [25]:
f ( t ) = 1 α K ( α ) s ( t ) + α K ( α ) Γ ( α ) a t s ( z ) ( t z ) α 1 d z .

3. Model Formulation

We begin to formulate the Ebola epidemic disease by considering the human population in three compartments, that is, the susceptible individuals, S ( t ) , individuals infected with Ebola virus, I ( t ) , and the individuals recovered from the Ebola virus, R ( t ) . The individuals infected with Ebola and the deceased is D ( t ) and P ( t ) is the class for the Ebola virus pathogen in the environment. The model that describes the dynamics of Ebola disease modeled through differential equations is given by
d S d t = Λ λ S d S , d I d t = λ S ( d + δ + ϕ 1 ) I , d R d t = ϕ 1 I d R , d D d t = ( d + δ ) I ε D , d P d t = ϖ + ξ I + θ D κ P ,
where λ = β 1 I + β 2 D + ψ P , and the appropriate initial conditions are given by
S ( 0 ) = S 0 , I ( 0 ) = I 0 , R ( 0 ) = R 0 , D ( 0 ) = D 0 , P ( 0 ) = P 0 .
The birth rate of the susceptible individuals is recruited by the rate Λ , while the death rate is given by d. The susceptible individuals become infectious with the effective contact rate β 1 and β 2 with the deceased human individuals. The susceptible are able to attract the disease from the contaminated environment at a rate given by ψ . The death rate of the infected individuals due to Ebola virus is given by a rate δ , while the recovery from infection is ϕ 1 . The deceased people can be directly buried during funerals at rate ε . At a rate of ϖ the environment is contaminated by the Ebola virus. At rates of ξ and θ the infected and deceased individuals, respectively, shed the virus in the environment. The virus decay of the Ebola virus from the population is given by parameter κ .
The sum of the first three equations of the Ebola disease model Equation (14) is given by
d N d t = Λ d N δ I ,
where N = S + I + R denotes the total alive human population. It should be noted that ε ( d + δ ) , which is an appropriate condition for the compartment D for which the model becomes relevant, otherwise the deceased human individuals will disappear and the model would be irrelevant. Further, the model given by Equation (14) is well posed and biologically feasible in the region given by
Φ = M R + 5 : N ( t ) Λ d , D ( d + δ ) Λ ε d , P ( t ) = ε ( d ϖ + Λ ξ ) + θ Λ ( δ + d ) d κ ε ,
where M = ( S ( t ) , I ( t ) , R ( t ) , D ( t ) , P ( t ) ) .

4. Ebola Model in the Caputo Sense

The purpose of this section is to apply the proposed three operators on the Ebola disease model Equation (14). Initially, we will apply the Caputo derivative on the Ebola disease model, then, the Caputo–Fabrizio derivative, and finally the Atangana–Baleanu derivative. For each operator we will provide the solution procedure and the discussion on the graphical results in details. So, we start with the Caputo sense.

4.1. Ebola Model in the Caputo Sense

We can express the model given by Equation (14) in the Caputo derivative as follows:
0 C D t α S = Λ λ S d S , 0 C D t α I = λ S ( d + δ + ϕ 1 ) I , 0 C D t α R = ϕ 1 I d R , 0 C D t α D = ( d + δ ) I ε D , 0 C D t α P = ϖ + ξ I + θ D κ P ,
where λ = β 1 I + β 2 D + ψ P , and with the initial conditions, S ( 0 ) = S 0 , I ( 0 ) = I 0 , R ( 0 ) = R 0 , D ( 0 ) = D 0 , and P ( 0 ) = P 0 .

4.2. Equilibrium Points

For the Ebola disease model Equation (18) in the Caputo sense, there is no disease-free equilibrium when ϖ > 0 and we have the other equilibrium say, E = ( S , I , R , D , P ) , we have
S = Λ ( d + δ + ϕ 1 ) I d , R = ϕ 1 I d , D = ( d + δ ) I ϵ , P = θ I ( d + δ ) + I ξ ϵ + ϖ ε κ ε
Using these values in the second equation of the model Equation (18), we have
C 2 I 2 C 1 I C 0 = 0 ,
where
C 2 = d + δ + ϕ 1 β 1 κ ε + β 2 κ ( d + δ ) + ψ ( θ ( d + δ ) + ξ ε ) , C 1 = ( κ Λ β 1 ε + β 2 ( d + δ ) ψ ( θ Λ ( d + δ ) ϖ ε ( d + δ ) + Λ ξ ε ) + d κ ε ( d + δ ) ) ε ϕ 1 ( d κ + ϖ ψ ) , C 0 = Λ ϖ ψ ε .
We have from the coefficient C 1 ,
C 1 = d ε κ ( d + δ + ϕ 1 ) R 0 1 ψ ϖ d κ ,
where
R 0 = Λ β 1 d ( d + δ + ϕ 1 ) + Λ β 2 ( d + δ ) d ε ( d + δ + ϕ 1 ) + ψ Λ ( ε ξ + θ δ + θ d ) d ε κ ( d + δ + ϕ 1 ) .
Considering the case when ϖ = 0 , we have
S = Λ d R 0 , I = Λ ( R 0 1 ) ( d + δ + ϕ 1 ) R 0 , R = Λ ϕ 1 ( R 0 1 ) d ( d + δ + ϕ 1 ) R 0 , D = ( Λ ( d + δ ) ) ( R 0 1 ) ε ( d + δ + ϕ 1 ) R 0 , P = Λ ( ε ξ + θ ( d + δ ) ) ( R 0 1 ) ε κ ( d + δ + ϕ 1 ) R 0 ,
and we have a disease-free equilibrium,
E 0 = ( Λ d , 0 , 0 , 0 ) ,
known as Ebola virus-free equilibrium.

4.3. Numerical Procedure for the Ebola Disease Model in the Caputo Sense

In the present subsection, we present the numerical scheme for the solution of the fractional Ebola disease model in the Caputo sense Equation (18). The present scheme that we use for the solution of the fractional Caputo nonlinear ordinary differential equation has been presented in [26,27]. The following procedure is presented
0 C D t α z ( t ) = f ( t , z ( t ) ) .
Using the fundamental theorem on Equation (20), we obtain
z ( t ) z ( 0 ) = 1 Γ ( α ) 0 t f ( χ , z ( χ ) ) ( t χ ) α 1 d χ ,
thus, at t = t n + 1 , n = 0 , 1 , , the following is obtained
z ( t n + 1 ) z ( 0 ) = 1 Γ ( α ) 0 t n + 1 ( t n + 1 t ) α 1 f ( t , z ( t ) ) d t ,
and
z ( t n ) z ( 0 ) = 1 Γ ( α ) 0 t n ( t n t ) α 1 f ( t , z ( t ) ) d t .
From Equations (23) and (22), we have
z ( t n + 1 ) = z ( t n ) + 1 Γ ( α ) 0 t n + 1 ( t n + 1 t ) α 1 f ( t , z ( t ) ) d t A α , 1 1 Γ ( α ) 0 t n ( t n t ) α 1 f ( t , z ( t ) ) d t A α , 2 .
where
A α , 1 = 1 Γ ( α ) 0 t n + 1 ( t n + 1 t ) α 1 f ( t , z ( t ) ) d t ,
and
A α , 2 = 1 Γ ( α ) 0 t n ( t n t ) α 1 f ( t , z ( t ) ) d t .
Using the Lagrange approximation for the function f ( t , z ( t ) ) , we have
P ( t ) t t n 1 t n t n 1 f ( t n , z n ) + t t n t n 1 t n f ( t n 1 , z n 1 ) = f ( t n , z n ) h ( t t n 1 ) f ( t n 1 , z n 1 ) h ( t t n ) .
The use of the above expression leads to
A α , 1 = f ( t n , z n ) h Γ ( α ) 0 t n + 1 ( t n + 1 t ) α 1 ( t t n 1 ) d t f ( t n 1 , z n 1 ) h Γ ( α ) 0 t n + 1 ( t n + 1 t ) α 1 ( t t n ) d t .
We have, after further simplification
A α , 1 = f ( t n , z n ) h Γ ( α ) 2 h α t n + 1 α t n + 1 α + 1 α + 1 f ( t n 1 , z n 1 ) h Γ ( α ) h α t n + 1 α 1 α + 1 t n + 1 α + 1 .
Similarly,
A α , 2 = 1 Γ ( α ) 0 t n ( t n t ) α 1 [ f ( t n , z n ) h ( t t n 1 ) f ( t n 1 , z n 1 ) h ( t t n ) ] d t .
Further simplifying, we get
A α , 2 = f ( t n , z n ) h Γ ( α ) h α t n α t n α + 1 α + 1 + f ( t n 1 , z n 1 ) h Γ ( α ) 1 α + 1 t n α + 1 .
We have the final approximate solution for the fractional nonlinear ordinary differential equation by substituting the Equations (30) and (31) into (24), given by
z ( t n + 1 ) = z ( t n ) + f ( t n , z n ) h Γ ( α ) 2 h t n + 1 α α t n + 1 α + 1 α + 1 + h α t n α t n + 1 α + 1 α + 1 + f ( t n 1 , z n 1 ) h Γ ( α ) h α t n + 1 α + t n + 1 α + 1 α + 1 + t n α + 1 α + 1 .
The above scheme is used further for the solution of the Ebola disease model in the Caputo sense Equation (18) by considering the parameter values, d = 0.05 , δ = 0.05 , ϕ = 0.06 , Λ = 10 , ϵ = 0.008 , ξ = 0.004 , κ = 0.03 , ψ = 0.01 , β 1 = 0.006 , β 2 = 0.012 , ϖ = 1 , and θ = 0.004 , and with various values of the fractional order parameter α . We have the graphical results for the numerical solution of the Ebola disease model Equation (18) in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7. One can observe in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 by deceasing the value of α , the individuals infected with Ebola decreases while the population of infected individuals increases. We use this further to check the graphical results for the case when α = 0.3 , 0.1 , then, one can see that infection is almost on the steady state, see Figure 6 and Figure 7.

4.4. Ebola Model in the Caputo–Fabrizio Sense

We can express the model given by Equation (14) in Caputo–Fabrizio derivative as follows:
0 C F D t α S = Λ λ S d S , 0 C F D t α I = λ S ( d + δ + ϕ 1 ) I , 0 C F D t α R = ϕ 1 I d R , 0 C F D t α D = ( d + δ ) I ε D , 0 C F D t α P = ϖ + ξ I + θ D κ P ,
where λ = β 1 I + β 2 D + ψ P , and with the initial conditions, S ( 0 ) = S 0 , I ( 0 ) = I 0 , R ( 0 ) = R 0 , D ( 0 ) = D 0 , and P ( 0 ) = P 0 .

4.5. Numerical Solution for Caputo–Fabrizio Model

Here we present the numerical solution for the Caputo–Fabrizio model Equation (33) by using the scheme presented [27]. The following steps are taken as specified in one for the solution of Equation (33).
S ( t ) S ( 0 ) = ( 1 α ) B ( α ) F 1 ( t , S ) + α B ( α ) 0 t F 1 ( ζ , S ) d ζ .
For t = t n + 1 , n = 0 , 1 , 2 , , we obtain
S ( t n + 1 ) S 0 = 1 α B ( α ) F 1 ( t n , S n ) + α B ( α ) 0 t n + 1 F 1 ( t , S ) d t .
The successive terms difference is given as follows:
S n + 1 S n = 1 α B ( α ) { F 1 ( t n , S n ) F 1 ( t n 1 , S n 1 ) } + α B ( α ) t n t n + 1 F 1 ( t , S ) d t .
Over the close interval [ t k , t ( k + 1 ) ] , the function F 1 ( t , S ) can be approximated by the interpolation polynomial
P k ( t ) f ( t k , y k ) h ( t t k 1 ) f ( t k 1 , y k 1 ) . h ( t t k ) ,
where h = t n t n 1 . Calculating the integral in Equation (36) using above polynomial approximation we get
t n t n + 1 F 1 ( t , S ) d t = t n t n + 1 F 1 ( t n , S n ) h ( t t n 1 ) F 1 ( t n 1 , S n 1 ) h ( t t n ) d t = 3 h 2 F 1 ( t n , S n ) h 2 F 1 ( t n 1 , S n 1 ) .
Substituting Equation (38) in (36) and after simplification we obtain
S n + 1 = S n + 1 α B ( α ) + 3 h 2 B ( α ) F 1 ( t n , S n ) 1 α B ( α ) + α h 2 B ( α ) F 1 ( t n 1 , S n 1 ) .
In a similar way, for the rest of equations of system Equation (33) we obtain the recursive formula as below
I n + 1 = I 0 + 1 α B ( α ) + 3 h 2 B ( α ) F 2 ( t n , I n ) 1 α B ( α ) + α h 2 B ( α ) F 2 ( t n 1 , I n 1 ) , R n + 1 = R 0 + 1 α B ( α ) + 3 h 2 B ( α ) F 3 ( t n , R n ) 1 α B ( α ) + α h 2 B ( α ) F 3 ( t n 1 , R n 1 ) , D n + 1 = D 0 + 1 α B ( α ) + 3 h 2 B ( α ) F 4 ( t n , D n ) 1 α B ( α ) + α h 2 B ( α ) F 4 ( t n 1 , D n 1 ) , P n + 1 = P 0 + 1 α B ( α ) + 3 h 2 B ( α ) F 5 ( t n , P n ) 1 α B ( α ) + α h 2 B ( α ) F 5 ( t n 1 , P n 1 ) .
The numerical scheme presented above by considering the parameter values, d = 0.05 , δ = 0.05 , ϕ = 0.06 , Λ = 10 , ϵ = 0.008 , ξ = 0.004 , κ = 0.03 , ψ = 0.01 , β 1 = 0.006 , β 2 = 0.012 , ϖ = 1 , and θ = 0.004 , we have the graphical results for the Ebola disease model in Caputo–Fabrizio model Equation (33). Various graphical results considering the fractional order parameter α are presented, see Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14. In these figures, we obtain various graphical results for α , and we observe that by decreasing the value of α the infected compartments are decreasing well. Especially, when choosing α = 0.3 , 0.1 , we can see that the number of infected individuals decreases rapidly.

4.6. Ebola Model in the Atangana–Baleanu Sense

We can express the model given by Equation (14) in Atangana–Baleanu derivative as follows:
0 A B C D t α S = Λ λ S d S , 0 A B C D t α I = λ S ( d + δ + ϕ 1 ) I , 0 A B C D t α R = ϕ 1 I d R , 0 A B C D t α D = ( d + δ ) I ε D , 0 A B C D t α P = ϖ + ξ I + θ D κ P ,
where λ = β 1 I + β 2 D + ψ P .

4.7. Existence of Solutions for the Atangana–Baleanu Model

It is obvious that the given model Equation (14) shows the dynamics of Ebola disease, which is described by a nonlinear system of differential equations, so it is not possible to obtain their exact solution but the existence of an approximate solution can be very effective if we show that the solution for the Ebola disease model Equation (41) under some conditions exists. To do this, we follow the results of fixed point theory for the given Ebola disease model Equation (41). We write the Ebola disease model given by Equation (41) for simplification purposes as follows:
0 A B C D t α x ( t ) = F ( t , x ( t ) ) , x ( 0 ) = x 0 , 0 < t < T < ,
where x ( t ) = ( S , I , R , D , P ) represent the vector with state variables S , I , R , D , P and is a continuous vector function and can be defined as follows:
F = F 1 F 2 F 3 F 4 F 5 = Λ λ S d S λ S ( d + δ + ϕ 1 ) I ϕ 1 I d R ( d + δ ) I ε D ϖ + ξ I + θ D κ P .
The function F can be shown easily to satisfy the Lipschitz condition and can be represented as:
F ( t , x 1 ( t ) ) F ( t , x 2 ( t ) ) K x 1 ( t ) x 2 ( t ) .
Now we have the results for the existence and uniqueness for the Ebola disease model in the Atangana–Baleanu derivative sense. We state and prove the following theorem:
Theorem 3.
The Ebola disease model in the Atangana–Baleanu form Equation (41) can have a unique solution under some conditions if the following holds
( 1 α ) A B C ( α ) K + α T m a x α A B C ( α ) Γ ( α ) K < 1 .
Proof. 
The use of the Atangana–Baleanu fractional integration on model Equation (42) both sides, the following is obtained,
x ( t ) = x 0 + 1 α A B C ( α ) F ( t , x ( t ) ) + α A B C ( α ) Γ ( α ) 0 t ( t η ) α 1 F ( η , x ( η ) ) d η .
Suppose J = ( 0 , T ) and the operator Y : C ( J , R 5 ) C ( J , R 5 ) defined by
Y [ x ( t ) ] = x 0 + 1 α A B C ( α ) F ( t , x ( t ) ) + α A B C ( α ) Γ ( α ) 0 t ( t η ) α 1 F ( η , x ( η ) ) d η .
Then we can write Equation (46) as follows:
x ( t ) = Y [ x ( t ) ] .
We have, after applying the supremum norm on J,
x ( t ) J = s u p t J x ( t ) , x ( t ) C .
Obviously, C ( J , R 5 ) and the norm . J is a Banach space. Using the operator Equation (48), the following is presented
Y [ x 1 ( t ) ] Y [ x 2 ( t ) ] J ( 1 α ) A B C ( α ) ( F ( t , x 1 ( t ) ) F ( t , x 2 ( t ) ) + α A B C ( α ) Γ ( α ) × 0 t ( t η ) α 1 ( F ( η , x 1 ( η ) ) F ( η , x 2 ( η ) ) ) d η J .
Using the triangular inequality and Lipschitz condition presented in Equation (44) with some simplifications, we have
Y [ x 1 ( t ) ] Y [ x 1 ( t ) ] J ( 1 α ) K A B C ( α ) + α A B C ( α ) Γ ( α ) K T α x 1 ( t ) x 2 ( t ) J .
Finally, we have
Y [ x 1 ( t ) ] Y [ x 1 ( t ) ] J L x 1 ( t ) x 2 ( t ) J ,
where
L = ( 1 α ) M A B C ( α ) + α A B C ( α ) Γ ( α ) M T α .
If the condition given by Equation (45) holds then the operator Y will be a contraction. Thus, the Banach fixed point theorem ensures that a unique solution for the Ebola disease model in the Atangana–Baleanu form Equation (41) exists, Equation (42). □

4.8. Numerical Results for the Atangana–Baleanu Model and Simulation Results

In the present subsection we aim to obtain the numerical results for the Ebola disease model in the Atangana–Baleanu form given by Equation (41). First, we provide a numerical scheme in details and then show the graphical results for various values of the fractional order parameter α . The scheme given in [28] will be used to obtain the approximate solution of the Ebola disease model in the Atangana–Baleanu form Equation (41).
We write the model Equation (42) after using the fundamental theorem of fractional calculus:
w ( t ) w ( 0 ) = ( 1 α ) A B C ( α ) F ( t , w ( t ) ) + α A B C ( α ) × Γ ( α ) 0 t F ( ϕ , x ( ϕ ) ) ( t ϕ ) α 1 d ϕ .
At t = t n + 1 , n = 0 , 1 , 2 , , we have
w ( t n + 1 ) w ( 0 ) = 1 α A B C ( α ) F ( t n , w ( t n ) ) + α A B C ( α ) × Γ ( α ) 0 t n + 1 F ( ϕ , w ( ϕ ) ) ( t n + 1 ϕ ) α 1 d ϕ , = 1 α A B C ( α ) F ( t n , w ( t n ) ) + α A B C ( α ) × Γ ( α ) j = 0 n t j t j + 1 F ( ϕ , w ( ϕ ) ) ( t n + 1 ϕ ) α 1 d ϕ .
The function F ( ϕ , w ( ϕ ) ) can be approximated over [ t j , t j + 1 ] , using the interpolation polynomial
F ( ϕ , w ( ϕ ) ) F ( t j , w ( t j ) ) h ( t t j 1 ) F ( t j 1 , w ( t j 1 ) ) h ( t t j ) .
Substituting in Equation (54) we get
w ( t n + 1 ) = w ( 0 ) + 1 α A B C ( α ) F ( t n , w ( t n ) ) + α A B C ( α ) × Γ ( α ) j = 0 n ( F ( t j , w ( t j ) ) h t j t j + 1 ( t t j 1 ) ( t n + 1 t ) α 1 d t F ( t j 1 , w ( t j 1 ) ) h t j t j + 1 ( t t j ) ( t n + 1 t ) α 1 d t ) .
After some calculation, we obtain the following:
w ( t n + 1 ) = w ( t 0 ) + 1 α A B C ( α ) F ( t n , w ( t n ) ) + α A B C ( α ) j = 0 n ( h α F ( t j , w ( t j ) ) Γ ( α + 2 ) ( ( n + 1 j ) α ( n j + 2 + α ) ( n j ) α ( n j + 2 + 2 α ) ) h α F ( t j 1 , w ( t j 1 ) ) Γ ( α + 2 ) ( ( n + 1 j ) α + 1 ( n j ) α ( n j + 1 + α ) ) ) .
For the Ebola disease model we have the following results:
S ( t n + 1 ) = S ( t 0 ) + 1 α A B C ( α ) F 1 ( t n , w ( t n ) ) + α A B C ( α ) j = 0 n ( h α F 1 ( t j , w ( t j ) ) Γ ( α + 2 ) ( ( n + 1 j ) α ( n j + 2 + α ) ( n j ) α ( n j + 2 + 2 α ) ) h α F 1 ( t j 1 , w ( t j 1 ) ) Γ ( α + 2 ) ( ( n + 1 j ) α + 1 ( n j ) α ( n j + 1 + α ) ) ) ,
I ( t n + 1 ) = I ( t 0 ) + 1 α A B C ( α ) F 2 ( t n , w ( t n ) ) + α A B C ( α ) j = 0 n ( h α F 2 ( t j , w ( t j ) ) Γ ( α + 2 ) ( ( n + 1 j ) α ( n j + 2 + α ) ( n j ) α ( n j + 2 + 2 α ) ) h α F 2 ( t j 1 , w ( t j 1 ) ) Γ ( α + 2 ) ( ( n + 1 j ) α + 1 ( n j ) α ( n j + 1 + α ) ) ) ,
R ( t n + 1 ) = R ( t 0 ) + 1 α A B C ( α ) F 3 ( t n , w ( t n ) ) + α A B C ( α ) j = 0 n ( h α F 3 ( t j , w ( t j ) ) Γ ( α + 2 ) ( ( n + 1 j ) α ( n j + 2 + α ) ( n j ) α ( n j + 2 + 2 α ) ) h α F 3 ( t j 1 , w ( t j 1 ) ) Γ ( α + 2 ) ( ( n + 1 j ) α + 1 ( n j ) α ( n j + 1 + α ) ) ) ,
D ( t n + 1 ) = D ( t 0 ) + 1 α A B C ( α ) F 4 ( t n , w ( t n ) ) + α A B C ( α ) j = 0 n ( h α F 4 ( t j , w ( t j ) ) Γ ( α + 2 ) ( ( n + 1 j ) α ( n j + 2 + α ) ( n j ) α ( n j + 2 + 2 α ) ) h α F 4 ( t j 1 , w ( t j 1 ) ) Γ ( α + 2 ) ( ( n + 1 j ) α + 1 ( n j ) α ( n j + 1 + α ) ) ) ,
P ( t n + 1 ) = P ( t 0 ) + 1 α A B C ( α ) F 5 ( t n , w ( t n ) ) + α A B C ( α ) j = 0 n ( h α F 5 ( t j , w ( t j ) ) Γ ( α + 2 ) ( ( n + 1 j ) α ( n j + 2 + α ) ( n j ) α ( n j + 2 + 2 α ) ) h α F 5 ( t j 1 , w ( t j 1 ) ) Γ ( α + 2 ) ( ( n + 1 j ) α + 1 ( n j ) α ( n j + 1 + α ) ) ) .
We using the above scheme for the numerical solution of the Ebola disease model Equation (41) and obtain the graphical results shown in Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20 and Figure 21 by considering different values of the fractional order parameter α . In these Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20 and Figure 21, by decreasing the values of the fractional order parameter α , the population of infected compartments decreases more efficiently for the cases of α = 0.3 , 0.1 . One can see that the numerical results in the form of graphs obtained through the Atangana–Baleanu operator in comparison to the Caputo and Caputo–Fabrizio operator decrease the infection faster. So, from the above graphical results, it is suggested that the Atangana–Baleanu operator is more useful for infection elimination by decreasing the value of α .

4.9. Graphical Comparison of the Operators

Here, we provide comparison plots for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu operators. We considered various values of the fractional order parameter α = 1 , 0.9 , 0.7 , 0.5 , 0.3 , 0.1 and presented the graphical results for comparison (see Figure 22, Figure 23, Figure 24, Figure 25, Figure 26 and Figure 27). By decreasing the value of α it can be seen that the number of infected individuals decreases well, compared to the Caputo–Fabrizio and Caputo derivatives. Especially, for the cases when α = 0.3 , 0.1 , the Atangana–Baleanu derivative provides useful results for the infection elimination in comparison to the Caputo and the Caputo–Fabrizio operators.

5. Conclusions

We presented the dynamics of an Ebola disease model in the framework of fractional calculus. We applied three fractional operators, which are the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu models. Initially, we proposed an epidemic model available literature for Ebola disease and then applied the proposed operators. The Ebola disease model with the Caputo derivative is presented and an effective numerical scheme for the numerical solution was provided. We used many values for the fractional order parameters and obtained the graphical results. The same model is used further and applied to the Caputo–Fabrizio derivative and we then presented a numerical solution for their solution. The solution was obtained and presented in graphical shape with the use of various fractional order parameter values. Then the newly established derivative known as the Atangana–Baleanu derivative was successfully applied to the Ebola disease model. The Ebola disease model in the Atangana–Baleanu sense is used and the uniqueness and existence were presented. Then, we presented a numerical scheme for the solution and presented various graphical results for α . Comparisons of the proposed three operators for various values of the fractional order parameter α = 1 , 0.9 , 0.7 , 0.5 , 0.3 , 0.1 are presented and discussed. The comparison results show that the Atangana–Baleanu derivative is more helpful for disease elimination by decreasing the value of α , since the population of infected individuals decreased well. The use of three different fractional operators on the Ebola disease model suggests that the fractional order parameter greatly affects disease elimination for the non-integer case when decreasing α . Therefore, we suggest that the application of the various fractional derivatives on the present disease model shows the greater effectiveness of the arbitrary order derivative than that of the integer order model for the case of fractional order parameters.

Author Contributions

Conceptualization, K.M.A. and A.A.; methodology, A.A.; software, K.M.A.; validation, K.M.A. and A.A.; formal analysis, A.A.; investigation, A.A.; resources, K.M.A.; data curation, A.A.; writing–original draft preparation, A.A. and K.M.A.; writing–review and editing, K.M.A.; visualization, A.A.; supervision, A.A.; project administration, A.A.; funding acquisition, A.A.

Funding

No source of funding.

Acknowledgments

The authors are thankful to the anonymous reviewers and handling editor for the careful reading and suggestions.

Conflicts of Interest

No conflict of interests exists regarding the publishing of this work.

References

  1. Ebola (Ebola Virus Disease). The Centers for Disease Control and Prevention. Available online: http://www.cdc.gov/ebola/resources/virus-ecology.html (accessed on 1 August 2014).
  2. Bibby, K.; Casson, L.W.; Stachler, E.; Haas, C.N. Ebola virus persistence in the environment: State of the knowledge and research needs. Environ. Sci. Technol. Lett. 2015, 2, 2–6. [Google Scholar] [CrossRef]
  3. Piercy, T.J.; Smither, S.J.; Steward, J.A.; Eastaugh, L.; Lever, M.S. The survival of floviruses in liquids, on solid substrates and in a dynamic aerosol. J. Appl. Microbiol. 2010, 109, 1531–1539. [Google Scholar] [PubMed]
  4. Leroy, E.M.; Rouquet, P.; Formenty, P.; Souquière, S.; Kilbourne, A.; Froment, J.M.; Bermejo, M.; Smit, S.; Karesh, W.; Swanepoel, R.; et al. Multiple Ebola virus transmission events and rapid decline of central African wildlife. Science 2004, 303, 387–390. [Google Scholar] [CrossRef] [PubMed]
  5. Leroy, E.M.; Kumulungui, B.; Pourrut, X.; Rouquet, P.; Hassanin, A.; Yaba, P.; Délicat, A.; Paweska, J.T.; Gonzalez, J.-P.; Swanepoel, R. Fruit bats as reservoirs of Ebola virus. Nature 2005, 438, 575–576. [Google Scholar] [CrossRef] [PubMed]
  6. Althaus, C. Estimating the reproduction number of Ebola (EBOV) during outbreak in West Africa. PLoS Curr. 2014. [Google Scholar] [CrossRef] [PubMed]
  7. Chowell, G.; Hengartner, N.W.; Castillo-Chavez, C.; Fenimore, P.W.; Hyman, J.M. The basic reproductive number of Ebola and the e?ects of public health measures: The cases of Congo and Uganda. J. Theor. Biol. 2004, 229, 119–126. [Google Scholar] [CrossRef] [PubMed]
  8. Fasina, F.O.; Shittu, A.; Lazarus, D.; Tomori, O.; Simonsen, L.; Viboud, C.; Chowell, G. Transmission dynamics and control of Ebola virus disease outbreak in Nigeria, July to September 2014. Euro Surveill. 2014, 19, 20920. Available online: http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=20920 (accessed on 1 August 2018). [CrossRef] [Green Version]
  9. Fisman, D.; Khoo, E.; Tuite, A. Early epidemic dynamics of the Western African 2014 Ebola outbreak: Estimates derived with a simple two Parameter model. PLoS Curr. 2014. [Google Scholar] [CrossRef] [PubMed]
  10. Ivorra, B.; Ngom, D.; Ramos, A.M. Be-CoDiS: A mathematical model to predict the risk of human diseases spread between countries-validation and application to the 2014–2015 ebola virus disease epidemic. Bull. Math. Biol. 2015, 77, 1668–1704. [Google Scholar] [CrossRef]
  11. Wang, X.-S.; Zhong, L. Ebola outbreak in West Africa: Real-time estimation and multiplewave prediction. Math. Biosci. Eng. 2015, 12, 1055–1063. [Google Scholar] [CrossRef]
  12. Lekone, P.E.; Finkenstädt, B.F. Statistical inference in a stochastic epidemic SEIR model with control intervention: Ebola as a case study. Biometrics 2006, 62, 1170–1177. [Google Scholar] [CrossRef] [PubMed]
  13. Berge, T.; Lubuma, J.M.-S.; Moremedi, G.M.; Morris, N.; Kondera-Shava, R. A simple mathematical model for Ebola in Africa. J. Biol. Dyn. 2017, 11, 42–74. [Google Scholar] [CrossRef]
  14. Zhang, Z.; Liu, C.; Zhan, X.; Lu, X.; Zhang, C.; Zhang, Y. Dynamics of Information Diffusion and Its Applications on Complex Networks. Phys. Rep. 2016, 651, 1–34. [Google Scholar] [CrossRef]
  15. Zhan, X.; Liu, C.; Zhou, G.; Zhang, Z.; Sun, G.; Zhu, J.J.H.; Jin, Z. Coupling dynamics of epidemic spreading and information diffusion on complex networks. Appl. Math. Comput. 2018, 332, 437–448. [Google Scholar] [CrossRef]
  16. Liu, C.; Zhan, X.; Zhang, Z.; Sun, G.; Hui, P.M. How events determine spreading patterns: Information transmission via internal and external influences on social networks. New J. Phys. 2015, 17, 113045. [Google Scholar] [CrossRef]
  17. Khan, M.A.; Ullah, S.; Okosun, K.O.; Shah, K. A fractional order pine wilt disease model with Caputo–Fabrizio derivative. Adv. Differ. Equ. 2018, 410. [Google Scholar] [CrossRef]
  18. Khan, M.A.; Ullah, S.; Chaos, M.F. A new fractional model for tuberculosis with relapse via Atangana–Baleanu derivative. Chaos Solitons Fractals 2018, 116, 227–238. [Google Scholar] [CrossRef]
  19. Ullah, S.; Khan, M.A.; Farooq, M. A fractional model for the dynamics of TB virus Chaos Solitons Fractals 2018, 116, 63–71. 116.
  20. Ullah, S.; Khan, M.A.; Farooq, M. Modeling and analysis of the fractional HBV model with Atangana–Baleanu derivative. Eur. Phys. J. Plus 2018, 133, 313. [Google Scholar] [CrossRef]
  21. Ullah, S.; Khan, M.A.; Farooq, M. A new fractional model for the dynamics of the hepatitis B virus using the Caputo–Fabrizio derivative. Eur. Phys. J. Plus 2018, 133, 237. [Google Scholar] [CrossRef]
  22. Diethelm, K. A fractional calculus based model for the simulation of an outbreak of dengue fever. Nonlinear Dyn. 2013, 71, 613–619. [Google Scholar] [CrossRef]
  23. Caputo, M.; Fabrizio, M. On the notion of fractional derivative and applications to the hysteresis phenomena. Meccanica 2017, 52, 3043–3052. [Google Scholar] [CrossRef]
  24. Losada, J.; Nieto, J.J. Properties of the new fractional derivative without singular Kernel. Progr. Fract. Differ. Appl. 2015, 1, 87–92. [Google Scholar]
  25. Atangana, A.; Baleanu, D. New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model. arXiv, 2016; arXiv:1602.03408. [Google Scholar]
  26. Khan, M.A. Neglecting nonlocality leads to unrealistic numerical scheme for fractional differential equation: Fake and manipulated results. Chaos 2019, 29, 013144. [Google Scholar] [CrossRef]
  27. Atangana, A.; Owolabi, K.M. New numerical approach for fractional differential equations. Math. Model. Nat. Phenom. 2018, 13. [Google Scholar] [CrossRef]
  28. Toufik, M.; Atangana, A. New numerical approximation of fractional derivative with non-local and non-singular kernel: Application to chaotic models. Eur. Phys. J. Plus 2017, 132, 444. [Google Scholar] [CrossRef]
Figure 1. The graphical results show the dynamics of the Caputo derivative model (18), when α = 1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 1. The graphical results show the dynamics of the Caputo derivative model (18), when α = 1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g001
Figure 2. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.95 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 2. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.95 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g002
Figure 3. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.9 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 3. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.9 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g003
Figure 4. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.85 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 4. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.85 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g004
Figure 5. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.5 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 5. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.5 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g005
Figure 6. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.3 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 6. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.3 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g006
Figure 7. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 7. The graphical results show the dynamics of the Caputo derivative model (18), when α = 0.1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g007
Figure 8. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 8. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g008
Figure 9. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.95 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 9. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.95 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g009
Figure 10. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.9 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 10. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.9 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g010
Figure 11. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.85 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 11. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.85 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g011
Figure 12. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.5 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 12. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.5 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g012
Figure 13. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.3 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 13. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.3 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g013
Figure 14. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 14. The graphical results show the dynamics of the Caputo–Fabrizio model (33), when α = 0.1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g014
Figure 15. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 15. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g015
Figure 16. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.95 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 16. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.95 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g016
Figure 17. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.9 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 17. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.9 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g017
Figure 18. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.85 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 18. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.85 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g018
Figure 19. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.5 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 19. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.5 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g019
Figure 20. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.3 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 20. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.3 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g020
Figure 21. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 21. The graphical results show the dynamics of the Atangana–Baleanu model (41), when α = 0.1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g021
Figure 22. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 22. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g022
Figure 23. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.9 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 23. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.9 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g023
Figure 24. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.7 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 24. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.7 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g024
Figure 25. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.5 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 25. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.5 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g025
Figure 26. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.3 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 26. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.3 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g026
Figure 27. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Figure 27. Comparison graphs for the Caputo, Caputo–Fabrizio, and the Atangana–Baleanu derivatives when α = 0.1 , where (a) Susceptible individuals, (b) Infected individuals, (c) recovered individuals, (d) deceased individuals, (e) Environment pathogens.
Entropy 21 00303 g027

Share and Cite

MDPI and ACS Style

Muhammad Altaf, K.; Atangana, A. Dynamics of Ebola Disease in the Framework of Different Fractional Derivatives. Entropy 2019, 21, 303. https://doi.org/10.3390/e21030303

AMA Style

Muhammad Altaf K, Atangana A. Dynamics of Ebola Disease in the Framework of Different Fractional Derivatives. Entropy. 2019; 21(3):303. https://doi.org/10.3390/e21030303

Chicago/Turabian Style

Muhammad Altaf, Khan, and Abdon Atangana. 2019. "Dynamics of Ebola Disease in the Framework of Different Fractional Derivatives" Entropy 21, no. 3: 303. https://doi.org/10.3390/e21030303

APA Style

Muhammad Altaf, K., & Atangana, A. (2019). Dynamics of Ebola Disease in the Framework of Different Fractional Derivatives. Entropy, 21(3), 303. https://doi.org/10.3390/e21030303

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