key: cord-0844381-gz3sdh0q
authors: Awais, Muhammad; Alshammari, Fehaid Salem; Ullah, Saif; Khan, Muhammad Altaf; Islam, Saeed
title: Modeling and Simulation of the Novel Coronavirus in Caputo Derivative
date: 2020-11-17
journal: Results Phys
DOI: 10.1016/j.rinp.2020.103588
sha: f1e8efb72139d38a2388fd07ff6d31b029ab75ad
doc_id: 844381
cord_uid: gz3sdh0q

The Coronavirus disease or COVID-19 is an infectious disease caused by a newly discovered coronavirus. The COVID-19 pandemic is an inciting panic for human health and economy as there is no vaccine or effective treatment so far. Different mathematical modeling approaches have been suggested to analyze the transmission patterns of this novel infection. this paper, we investigate the dynamics of COVID-19 using the classical Caputo fractional derivative. Initially, we formulate the mathematical model and then explore some the basic and necessary analysis including the stability results of the model for the case when [Formula: see text]. Despite the basic analysis, we consider the real cases of coronavirus in China from January 11, 2020 to April 9, 2020 and estimated the basic reproduction number as [Formula: see text]. The present findings show that the reported data is accurately fit the proposed model and consequently, we obtain more realistic and suitable parameters. Finally, the fractional model is solved numerically using a numerical approach and depicts many graphical results for the fractional order of Caputo operator. Furthermore, some key parameters and their impact on the disease dynamics are shown graphically.

Coronavirus disease (COVID-19) is a highly contagious viral disease that affects million of humans and caused reasonable death cases around the world. The infection was named by the World Health Organization (WHO) as COVID- 19 . The infection for the first time COVID-19 cases and deaths that have accounted due to coronavirus infection. Recently, the role of non-pharmaceutical intervention on COVID-19 dynamics by considering the data of Pakistan has been studied in [9] . The authors used the real data of coronavirus of Pakistan and implemented a mathematical model with optimal control analysis. The fractional epidemic model is another strong tool to explore the dynamics of infectious disease in a better way than the ordinary integer-order models. These models are based on fractional order differential operators, which generalizes the classical integerorder derivatives. Epidemic models with fractional derivatives give a greater degree of accuracy and provide a better fit to the real data in compression with classical integerorder models [10, 11] . A diverse variety of fractional operators were introduced time to time in the literature with a different kernel. Some of the frequently used fractional orders (FO) operators are Caputo [12] , Caputo-Fabrizio (CF) [13] and Atangana-Baleanu operator (AB) [14] . Although most of the COVID-19 models developed so far are based on classical integer-order derivative but only a few can be found with fractional operators. For example, in [15] Khan and Atangana formulated a fractional model of COVID-19 using Atangana-Bleanu in Caputo sense (ABC) operator and provided a better fit to the reported case in Wuhan. Recently, in [16] Abdo et al. extended a classical COVID-19 model in the literature to fractional order using Mittag-Leffler kernel and carried out the Ulam stability analysis and simulation results. Atangana studied a new COVID-19 transmission model using different fractal-fractional operators [17] . The results are obtained both mathematically and statistically considering the real cases. Further, the lockdown effect in the model has been analyzed. The authors in [18] formulated a mathematical model for coronavirus from natural to human host. Some recent mathematical studied on corona and other infectious diseases are studied by authors see [19, 20, 21, 22] . A novel results for the corona virus model has been discussed in [19] . The application of fractal-fractional differential equations to partial differential equations has been analyzed in [20] . The authors in [21] presented a theoretical results for fractal fractional differential equations. Fractal fractional differential equations and its solutions in detailed with novel analysis is presented in [22] . Inspired by the above discussion, in the present study, we reformulate the model [15] using the fractional operator in the Caputo sense. In [15] , they considered the reported cases for a short period of time from January 21, 2020 to January 28, 2020, in Wuhan city China. In this study, we consider the reported case in the mainland of China from January 21, 2020 to April 9, 2020 and explore a better theoretical and graphical analysis of the model. The data is taken from [23] . Using the reported data we will obtain the real parameter values for the model and will perform the simulation to predict the dynamics of coronavirus in the community. We will show that the model provides better fitting to the reported cases. Some important parameters that can help the disease under control will be explored in detail. The rest of the results are organized is as follows: The concepts and basic results regarding fractional and Caputo fractional derivative are recalled in Section 2. Formulation of the model, the integer case, the model, and their fitting to the coronavirus cases and parameters are discussed in Section 3. The COVID-19 model with Caputo fractional derivative, its basic mathematical properties as well as stability results are given in Section 4. The numerical iterative scheme and the simulation results of the fractional model are plotted and discussed in Section 5. Finally, a brief concluding remarks are given in Section 6.

In the following, we first recall some basic details of fractional calculus [12] . Definition 1. Let X (t) ∈ C n ([0, ∞+], R) be a function then the fractional derivative in Caputo sense having order n − 1 < p ≤ n where n ∈ N is given as:

Clearly, C D p t (X (t)) approaches to X ′ (t) as p → 1.

A fixed point say, x * satisfying f (t, x * ) = 0, then it is called the equilibrium point of the model formulated via Caputo derivative given by:

As proved in [24] , the point x * will be locally asymptotically stable if all eigenvalues Π of the Jacobian matrix calculated around the equilibrium point satisfy the following condition:

For the global stability results of the fractional system in Caputo derivative considering the Lyapunov method, we need to have the following important results [25, 26] .

Theorem 2.1. Suppose x * be an equilibrium point of the system (1) and Ψ ∈ R n be a domain containing x * , and let L : [0, ∞) × Ψ → R be a continuously differentiable function satisfying

and

∀ 0 < p < 1 and x ∈ Ψ. Where M 1 (x), M 2 (x) and M 3 (x) denote the continuous positive definite functions on Ψ. Then, x * will be uniformly asymptotically stable equilibrium point of the (1).

In this section, we briefly explain the model formulation using classical integer-order derivative. To develop the model total human population denoted by N(t) is further classified into five mutually-exclusive epidemiological classes, which are the susceptible S(t), exposed E(t), infected (symptomatic) I(t), asymptotically infected showing no clinical symptoms A(t), and the recovered people R(t). The class M(t) denotes the COVID-19 in reservoir or the seafood place or market. The recruitment rate is shown by the Λ while the natural death rate is represented by µ. The symptomatic and asymptomatic infected people could export the virus into M at the rate ̺ and ̟, respectively. The virus in M will subsequently leave the M class at a rate of ν, where 1/ν accounts the lifetime of the COVID-19 virus. The susceptible individuals acquired infection after effective contacts with the people in I and A compartments at the rate η(I+ψA)S N , where the parameter η is the disease transmission coefficient and ψ is the transmissibility multiple of A to I and 0 ≤ ψ ≤ 1. Further, we assume that susceptible people could be infected after the interaction with M, given by η w SM where, η w is the transmission rate from M to S. The exposed people develop the COVID-19 infection after the completion of incubation period and then move to either class I or A at the transmission rate ω and ρ respectively. The parameter θ is the proportion of asymptomatic infection. The infectious period of symptomatic I and asymptomatic A individuals are defined as 1/τ and 1/τ a respectively after which they enter to recovered class R. The nonlinear differential equations that describes the dynamics of the COVID-19 disease are given by [15] :

The corresponding initial conditions are

For simplicity, let us denote

Then, the above model can be written as

Next, we estimate the parameters of the model and provide the method with brief explanation.

We investigating the parameters estimations and curve fitting for the reported data of coronavirus cases from January 21, 2020 to April 9, 2020, in the mainland of China to the model (5) . We utilized the approach first by selecting some of available parameters from the literature and then the rest of the parameters are fitted to the coronavirus cases. We consider the time unit per day . In the following, we provide a brief explanations of the parameters estimation:

(i) Natural mortality rate µ: The life span an average in China for the year 2019 is 76.79 years [27] , so, the estimated value of natural death rate is µ = 1/76.79 per year.

(ii) Recruitment rate Λ: Population of China for the year 2019 is about N h (0)=1.43 billion [27] , so, we can obtain the recruitment rate parameter by computing Λ/µ = N h (0), and also under the assumption that in the disease absence it is the limiting population, so we have Λ = 46381 per day.

Using the above parameters values and the rest of the parameters are simulated and fitted to the actual cases that is depicted in Figure 1 using the method of least-square curve fitting briefly discussed and utilized in [28, 29] . The best fitted model parameters are calculated by minimizing the error between the reported infected data and the infected cases obtained from the model (5) . The objective function implemented in the estimation procedure is as follows:Θ

where I t l denotes the reported cumulative COVID-19 confirmed case and I t l shows the model solution for infected class at time t l , while m is the number of available reported data points. The cumulative reported infected cases in the mainland of China are shown in Figure 1 and the resulting best fit by the proposed model to field data is shown in Figure 2 . One can observe that the model provide very accurate fitting to the reported cases of coronavirus and considered to be more reasonable in order to do prediction based on the fitted parameters. Table 1 contains the fitted and the estimated parameters and the approximated value of the basic reproduction number using these parameter values is R 0 ≈ 4.95.

We extend the classical integer model by replacing the derivative with the Caputo fractional-order derivative in order to provide a better understanding of the COVID-19 infection. The parameter p ∈ (0, 1] in the model is used for the order of Caputo derivative. Hence, the model is given through the following system of nonlinear fractional differential equations:

with the initial conditions given in (6).

To proceed the model analysis, first we explore the basic necessary mathematical properties. For the non-negativity of the model solution, let us construct the following set Reported cases in mainland China We follow the approach used in [30] in order to prove the desired results and provide the following theorem.

Lemma 4.1. [30] Consider, H(t) ∈ C[r 1 , r 2 ] and C r 1 D p t H(t) ∈ C(r 1 , r 2 ], then

Corollary 1. [30] Suppose that H(t) ∈ C[r 1 , r 2 ] and C r 1 D p t H(t) ∈ C(r 1 , r 2 ], where p ∈ (0, 1], then the following conclusions are drawn if

Theorem 4.1. For the fractional model (9) there exists a solution, say, y(t) which is unique and will remain in R 6 + and non-negative.

Proof. Utilizing the approach presented in [31] , it is easy to show the existence of the Caputo model (9) . Further, the uniqueness of the solution can be obtained by making use of the remark 3.2 in [31] for all values of t > 0. For the positivity of the solution, it is necessary it is necessary to show that on each hyperplane bounding the positive orthant, the vector field points to R 6 + . From the system (9), we deduced that

Hence, based on the above corollary, we concluded that the solution will remain in R 6 + , for all t ≥ 0.

Finally, to confirm the boundedness of the fractional model, we proceed by summing the first five equations of the model (9) , which give

Applying the Laplace transform of (10), we obtained

Applying Laplace inverse, we obtained

The Mittag-Leffler function describe by infinite power series i.e;

and laplace transform of Mittag-leffler function is Thus keeping the fact in mind that the Mittag-Leffler function has an asymptotic behavior [12] , we obtain

Hence, the biologically feasible region is constructed as:

Next, we investigate the equilibria and the associated basic reproduction number R 0 .

The equilibrium points of the Caputo fractional model (9) can be evaluated by setting

The disease free equilibrium (DFE) denoted by D 0 is obtained by solving the system (12) at disease free state D 0 = S 0 , 0, 0, 0, 0 0 = Λ/µ, 0, 0, 0, 0, 0 .

Further, using next generation approach we evaluate the following expression of the basic reproduction number:

The endemic equilibrium (EE) of the Caputo model (9) is represented by D * * and is given as follows:

and λ * * = η(I * * + ψA * * ) S * * + E * * + I * * + A * * + R * * + η w M * * , satisfying the following equation

The coefficients in (13) are as follows:

Clearly, m 1 > 0 and m 2 ≥ 0 if R 0 < 1, and λ * * = − m 2 m 1 ≤ 0. Thus, no endemic equilibrium will exist if R 0 < 1.

In order to confirm the local stability of DFE it is enough to show that the eigenvalues of the matrix of dynamics given (16) lie outside the closed angular sector as proved in [24] .

The following theorem deals with the desired result. 

Proof. The Jacobian of linearization matrix of model is given as:

The characteristic equation in term of λ for J D 0 is given below:

where,

From (17) the argument of the first two eigenvalues −µ, satisfy the condition given in (15) for all p ∈ (0, 1]. Further, it is clear that if R 0 < 1, then all a i > 0, and it is easy to show that a 1 a 2 a 3 > a 2 3 + a 2 1 a 4 . Thus, following [24, 25] , the arguments of all eigenvalues satisfy the necessary condition given in (15) . This completes the proof.

In this part, we prove the global asymptotical stability (GAS) of the fractional model (9) . For this purpose, the following result is presented. Proof. To present the proof, let us define the following suitable Lyapunov function:

where A j , for j= 1, · · · , 4, which are positive constants to be consider later. Taking the time Caputo fractional derivative of F (t) we obtain

Using system (9), we get

and then after some simplification, we have,

It is clear that when R 0 < 1 then C D p t F (t) is negative, because all the model parameters are non-negative. Thus, it follows from the results given in Theorem 2.1, that the DFE D 0 is GAS in Ω.

This section investigates the numerical solution of the fractional model (9) and to present the simulation results for various values of model parameters and p ∈ (0, 1]. In order to do this, we utilize the Euler's type approach presented discussed in the recent literature and references therein [32, 33] . To obtain the iterative scheme, let us express the fractional model (9) in the following simple form:

where g = (S, E, I, A, R, M) ∈ R 6 , G t, g(t) is used for a continuous real valued vector function, which additionally satisfies the Lipschitz condition and g 0 stands for initial state vector. Taking Caputo integral on both sides of (18) we get

In order to formulate an iterative scheme, we consider a uniform grid on [0, T ] with h = T −0 m is the step size and m ∈ N. Thus, the equation (19) gets the structure as follows after make use of the Euler method [34]     

Thus, utilizing the above scheme (20) , we deduced the following iterative formulae for the corresponding classes of the model (9)

where

Further, we used the estimated parameters to obtain the simulations for the Caputo model (9) by varying fractional order p. The parameter values used in graphical interpretation are given in Table 1 . In plots, the time level is taken up to 150 days. The behavior of the model (9) for different values of fractional order p is described in Figure 3 . By making decrease in the value of arbitrary order p of the Caputo operator, a decrease in the infected compartments while increase in the susceptible compartment is observed. Figures  4 and 5 respectively describe the impact of disease transmission rate η w on the symptomatic infected (I(t)) and asymptomatic infected (A(t)) people for various cases of fractional order p. It is observed that the population in both infected classes is decreasing significantly with the decrease in η w as shown in subplots (a-d) of Figures 4 and 5 . Further, the same behavior is observed for all values of fractional order p but the decrease in infected individuals is more significant for the smaller values of p. 

The dynamics of novel coronavirus (COVID-19) pandemic is presented using a mathematical model. The model was first analyzed in [15] with Atangana-Baleanu fractional operator. In this study, we reformulated the model using the Caputo fractional operator and provides a better theoretical and graphical analysis of the model. Additionally, in the present paper, we considered the reported COVID-19 confirmed cases from the early stage until the end of pandemic in the mainland China. The findings of the present study show a better agreement of the model prediction to the real reported infected cases and ultimately, a more suitable model parameters are estimated. Initially, the model is briefly explored with classical integer-order derivative and then the Caputo type fractional derivative is applied to reformulate the model to best described the disease dynamics. The basic mathematical results of the fractional model such as solution positivity, feasibility, equilibria and the most important biological threshold parameter basic reproduction number is analyzed. Further, we have shown that the COVID-19 model with Caputo operator is locally and globally stable if R 0 < 1. Besides the basic analysis of the model, we further estimated the biological parameters using the reported confirmed infected cases in the mainland of China for the period January, 21 to April 9, 2020. It is worth mentioning that the present model fits very well with the real data of daily confirmed cases, as shown in Figure 2 which reflects the reality of this infection in China. Based on the fitted parameters, the estimated value of the reproduction number is obtained as R 0 ≈ 4.95. Furthermore, using the real estimated parameters we carried out the numerical simulation and depicted the graphical results of the fractional model for various values of p. We also analyzed the influence of model parameters on the symptomatic and asymptomatic infected individuals for four different values of fractional order p. It is observed that by decreasing the parameter η w (the disease transmission coefficient) the infected individuals decrease significantly. The same behavior is observed for the fractional order p, but the decrease in the infected individuals is more significant for smaller values of fractional order of the Caputo operator p. We believe the work presented in this paper for the mathematical analysis of coronavirus with real cases will be more useful for the health authorities or the other decision making agencies in combating the disease. The authors wish to continue this study in the near future by extending this work using the control measures by applying the optimal control theory in order to get the possible elimination or decrease the infection of COVID-19 pandemic. 

Writing-Original draft preparation, Muhammad Altaf Khan: Conceptualization, Methodology, Writing-Original draft preparation, Supervision, Reviewing and Editing, Saeed Islam: Conceptualization, Methodology, Writing-Original draft preparation

Optimal control analysis of the effect of treatment, isolation and vaccination on hepatitis B virus

Modeling and scientific computing for the transmission dynamics of Avian Influenza with Half-Saturated Incidence

A COVID-19 epidemic model with latency period

Analysis and forecast of COVID-19 spreading in China

Short-term forecasting COVID-19 cumulative confirmed cases: Perspectives for Brazil

Mathematical Modeling of COVID-19 Transmission Dynamics with a Case Study of Wuhan

Modeling the impact of non-pharmaceutical interventions on thedynamics of novel coronavirus with optimal control analysis with acase study

A fractional model for the dynamics of TB virus

A fractional model for the dynamics of competition between commercial and rural banks in Indonesia

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

A new definition of fractional derivative without singular kernel

New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model

Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative

On a Comprehensive Model of the Novel Coronavirus (COVID-19) Under Mittag-Leffler Derivative

Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination?

A Mathematical Model of the Evolution and Spread of Pathogenic Coronaviruses from Natural Host to Human Host

Analysis coronavirus disease (COVID-19) model using numerical approaches and logistic model

Modelling and analysis of fractal-fractional partial differential equations: Application to reaction-diffusion model

Analysis of fractal fractional differential equations

On solutions of fractal fractional differential equations

Stability results for fractional differential equations with applications to control processing

Stability analysis of Caputo fractional-order nonlinear systems revisited

Volterra-type Lyapunov functions for fractional-order epidemic systems

Modeling the transmission dynamics of tuberculosis in Khyber Pakhtunkhwa Pakistan

Optimal control strategies for dengue transmission in Pakistan

Generalized Taylor's formula

Global existence theory and chaos control of fractional differential equations

A new fractional analysis on the interaction of HIV with CD4+ T-cells

Classical and contemporary fractional operators for modeling diarrhea transmission dynamics under real statistical data

Numerical methods for fractional calculus

The authors have declared no conflict of interest.

Muhammad Awais: Conceptualization, Methodology, Writing-Original draft preparation, Fehaid Salem Alshammari: Conceptualization, Methodology, Writing-Original