key: cord-0920843-69aa8xwi
authors: Singh, Harendra
title: Analysis for fractional dynamics of Ebola virus model
date: 2020-06-16
journal: Chaos Solitons Fractals
DOI: 10.1016/j.chaos.2020.109992
sha: a1091d9d4cabfd74ef5114daf6c3ecc25a1cbc3f
doc_id: 920843
cord_uid: 69aa8xwi

Ebola virus is very challenging problem of the world. The main purpose of this work is to study fractional Ebola virus model. An efficient computational method based on iterative scheme is proposed to solve fractional Ebola model numerically. Stability of proposed method is also discussed. Efficiency of proposed method is shown by listing CPU time. Proposed computational method will work for long time domain. Numerical results are presented graphically. The main reason for using this technique is low computational cost and high accuracy. It is also shown how the approximate solution varies for fractional and integer order Ebola virus model.

Ebola virus disease, which was first discovered near Ebola River in Africa in 1976, is a serious disease with rare outbreaks. The Ebola virus affects humans and animals like monkeys, gorillas, and chimpanzees. At the initial stage, this virus was affecting people and causing outbreaks in some African countries. The origin of the Ebola virus was not determine exactly. The bats or other animals are considered as origin for this type of virus. This virus infects the humans by direct contact with the blood in the starting [1 , 2] . Now many researchers are working in the area of fractional calculus. Many analytical and numerical approach have been developed to solve problems in fractional calculus [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] . Fractional calculus has many real life applications like in engineering, epidemiology, fluid mechanics, physics, chemistry, medical and health sciences. The mathematical model based on Ebola outbreak in Liberia was studied by Rachah et al. [14] and the results were simulated using existing data provided by the World Health Organization. Area et al. [15] investigated fractional SEIR Ebola outbreak model. In [16] , authors used two models for optimal control of Ebola disease. Atangana and Goufo [17] , studied generalised version of Ebola virus model and also discussed the endemic equilibrium points. In [18] , authors discussed this model using Atangana Baleanu fractional derivative and also discussed existence and uniqueness of the solutions. Some other method on Ebola disease can be found in [19] [20] [21] . Some recent developments in the area of numerical methods to solve this type of models can be found in [22 , 23] .

In this paper we will use an iterative scheme to solve fractional Ebola virus model as given in [18] [19] [20] . Proposed iterative scheme is based on the discretization of the domain. We are us-E-mail address: harendra059@gmail.com ing Grünwald-Letnikov (GL) fractional derivative for the fractional Ebola virus model. Coefficients in the approximation of GL derivative is very easy in computation and provides long time behaviour of solution in the fraction of second. The solution behaviour is shown in long time interval. Stability is also discussed for proposed iterative scheme. Figures and tables are used to show the numerical results. CPU time taken in iterative scheme for different values of iterations are listed. Behaviour of approximate solutions for integer and fractional order derivatives are shown using figures. Finally some concluding remarks are given.

The present model is based on the fact that total population at time t is denoted by N ( t ) and divided into four subgroups. We write these subgroups as follows:

A ( t ) is the subgroup of susceptible people; B ( t ) is the subgroup of infected people; C ( t ) denote the subgroup of recovered people; D ( t ) is the subgroup of died people in that region. These are speci- 

https://doi.org/10.1016/j.chaos.2020.109992 0960-0779/© 2020 Elsevier Ltd. All rights reserved. Also, N represents the total population in the region, a 1 is the rate of recovery, b 1 is the rate of infection from Ebola, c 1 is the amount of death rate, d 1 is the rate of death from Ebola and e 1 is the rate of susceptibility. The Ebola virus model strongly depends on initial conditions and integer order Ebola virus model can not explain perfectly this model. In addition to show the effect of above mentioned parameters numerically and for better understand of Ebola virus model it is require to replace integer order Ebola virus model to fractional order model. In this paper we will replace time derivative in Ebola virus model with fractional time derivative given by

With initial conditions

where, 0 ≤ p 1 , p 2 , p 3 , p 4 < 1.

In this section, we will describe some basic preliminaries for Grünwald-Letnikov derivative and Jacobian of matrix.

If we consider n = t−a h , where a is a real constant, the Grünwald-Letnikov definition is defined as

where [.] denote the integer part. The fractional-order linear system can be written as

Using Eq. (4) , the p − th order Grünwald-Letnikov derivative at the points kh ( k = 1 , 2 , . . . ) has the following form

where L is "memory length", t k = kh , h is the step size of calculation and the coefficients c 

Then the general approximate solution of the equation

is given as

Short memory principle can be used for the memory term expressed by a sum and for the determination of lower index. Using short memory the lower index in sum is given as

For the L , we can use

The fractional order systems given by

where aD q t is fractional GL derivative and q ∈ (0, 1]. The equilibrium points for system of differential equation is obtained by solving

The Jacobian matrix of system is given by 

Theorem 3.1. [25 , 26] . The system in Eq. (11) , is locally stable if all the eigenvalues of Jacobian matrix given in Eq. (14) , will satisfy the condition

In this section we will apply our proposed method to get the approximate values to the unknowns in Ebola virus model. The integer order Ebola virus model is given as

Taking integral on both side Eq. (16) , we get

In this paper we consider fractional corona virus model by replacing integer order derivative by fractional order GL derivative and given as 

Using GL definition in Eq. (19) , we get following iterative relations between unknowns:

In Eq. (20) , we have got an iterative relations to get the approximate solution for fractional Ebola virus model. Now we will iterate above relations as required number of iterations to achieve the desired accuracy. As the step size will minimize then the number of iterations will increase and the solution become more accurate. We have taken step-size h = 0 . 01 in all numerical simulation to calculate the approximate solution for the group of fractional Ebola virus model. 

In this section we will discuss stability of our mathematical model. We will use concept of Jacobian matrix and equilibrium points to show the stability of our model.

The equilibrium points for system (18) is given by

Using Eq. (13) , the Jacobian matrix of above system is given as

Since we are calculating Endemic equilibrium end points which are characterized by the existence of infected nodes ( i.e . B = 0).

Since two coordinates are negative therefore it avoid the existence of this type equilibrium point. Supposing the existence of this type of equilibrium point. For the parameters, N = 1 , 0 0 0 , 

The ( 50 , −1 . 67 , −16 . 67 , 0 ) .

The Jacobian matrix at Endemic equilibrium point −1 . 67 , −16 . 67 , 0 ) is given as

The Eigen values corresponding to matrix J 2 are λ 1 = 0 , λ 2 = 0 . 1937 , λ 3 = −0 . 1770 , λ 4 = −0 . 0233 . J 2 has two negative Eigen values and one positive Eigen value. Two Eigen values are negative therefore it has a stable manifold of dimension two and unstable manifold of dimension one at Endemic equilibrium point −1 . 67 , −16 . 67 , 0 ) . Since it has a two dimensional stable manifold so from dynamical point view it has very importance. Consequently, if Endemic equilibrium points exist then it is unstable but from stability discussion it is clear that it has two dimensional stable manifold and one dimensional unstable manifold so from dynamical point of view it has importance.

From the Eq. (21) , it is clear that if we take B = 0 (i.e. disease free equilibrium) then C = 0 and we can not calculate A . Therefore disease free equilibrium is not possible for this model. Suppose disease free equilibrium exist and is denoted by ( A 0 , 0, 0, 0) . There is only one infected compartments i.e. group of infected people ( B ). The group of recovered people ( C ) and susceptible people ( A ) can not considered in infected compartments. The group of failure of treatment i.e. group of died people ( D ) can not also considered to be new infections. Now following the work in [27] , where the necessary computation of the matrices F and V are given by,

The spectral radius γ ( F V −1 ) is the required basic reproduction number of the model (2) , which is given by

In Fig. 1 , we have plotted the group of susceptible people ( A ) with time. From Fig. 1 , it is observed that it decreases with time and in long-time behaviour it come close to zero i.e. the number of susceptible people become very less in society. In Fig. 2 , we have plotted the group of infected people ( B ) with time. From Fig. 2 , it is observed that it increases suddenly with time and in long time it come close to zero i.e. the number of infected people become very less in society. In Fig. 3 , we have plotted the group of recovered people ( C ) with time. From Fig. 3 , it is observed that initially it increases i.e. the numbers of recovered people are increasing with time. In Fig. 4 , we have plotted the group of died people ( D ) with time. From Fig. 4 , it is observed that it increases suddenly with time and in long time it become constant. So our proposed model performances well because number of recovered people are increasing and dead people are decreasing with time. Fig. 5 , shows the dynamic of susceptible people; infected people; recovered people and dead people for integer order time. In Fig. 5 , the behaviour of A , B , C and D is shown simultaneously for integer order. From Fig. 5 , it is clear that the group of recovered people increases with time. In Fig. 6 , we have plotted the group of susceptible people A ( t ) with respect to time for different values of fractional order time derivatives. From Fig. 6 , it is observed that the group of susceptible people A ( t ) varies continuously with time derivative. In Fig. 7 , we have plotted infected people B ( t ) with respect to time for different values of fractional order time derivatives. From Fig. 7 , it is observed that group of infected people B ( t ) varies continuously with time derivative. In Fig. 8 , we have plotted group of recovered people C ( t ) with respect to time for different values of fractional order time derivatives.

From Fig. 8 , it is observed that group of recovered people C ( t ) varies continuously with time derivative. In Fig. 9 , we have plotted asymptotically died people D ( t ) with respect to time for different values of fractional order time derivatives. From Fig. 9 , it is observed that asymptotically died people D ( t ) varies continuously with time derivative.

In Tables 1 and 2 , we have listed CPU time taken for different values of t and n . From tables it is clear that our proposed method is excellent and time saving.

The proposed method is very easy to implement. The efficiency and accuracy of the proposed method is very attractive and clear from tables. From figures it is clear that results are very interest-ing showing the continuous behaviour of approximate solution for fractional order. Using the concept of Jacobian matrix and equilibrium points the stability of the model was performed. From the numerical solution and stability analysis discussion it is clear that dynamic of Ebola virus model for a given total number of population, T , in society at a time t , depends on its parameters. Their values are key to increase the number of recovered people in the society. The approximate solutions obtained by proposed method is effective for a long time so these solutions will be very useful in understanding Ebola virus model properly. The reduction in the size of the computational domain shows wider applicability of this method.

Author contributed in every section of the paper in best possible ways as well as read and approved the paper for the submission.

None.

The pathogenesis of Ebola virus disease

Historical analysis of the Ebola virus: prospective implications for primary care nursing today

Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications

Emerging science publishers

Applications of fractional calculus in physics

Theory and applications of fractional differential equations, north-holland mathematical studies

Jacobi collocation method for the approximate solution of some fractional-order Riccati differential equations with variable coefficients

Stable numerical approach for fractional delay differential equations

Error estimate for the numerical solution of reaction-sub diffusion process based on a mesh less method

A reliable numerical algorithm for fractional advection-dispersion equation arising in contaminant transport through porous media

A reliable numerical algorithm for the fractional vibration equation

A collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation

A new stable algorithm for fractional Navier-Stokes equation in polar coordinate

Mathematical modelling, simulation, and optimal control of the 2014 Ebola outbreak in West Africa

On a fractional order Ebola epidemic model

Predicting and controlling the Ebola infection

On the mathematical analysis of Ebola hemorrhagic fever: deathly infection disease in West African countries

Modelling the spread of Ebola virus with Atangana-Baleanu fractional operators

Optimal control application to an Ebola model

A fractional order model for Ebola virus with the new caputo fractional derivative without singular kernel

A case of Ebola virus infection

A fractional-order delay differential model for Ebola infection and CD8 + T-cells response: stability analysis and Hopf bifurcation

A novel method for a fractional derivative with non-local and non-singular kernel

Modeling and numerical analysis of fractional-order Bloch equations

Stability results for fractional differential equations with applications to control processing

Fractional dynamical system and its linearization theorem

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

The author is very grateful to the referees for their constructive comments and suggestions for the improvement of the paper.