key: cord-0689161-ub6zij26
authors: Zarin, Rahat; Khan, Amir; Yusuf, Abdullahi; Abdel‐Khalek, Sayed; Inc, Mustafa
title: Analysis of fractional COVID‐19 epidemic model under Caputo operator
date: 2021-03-25
journal: Math Methods Appl Sci
DOI: 10.1002/mma.7294
sha: 1c10b37f21e5da5a8d93b956b589485ebb04de43
doc_id: 689161
cord_uid: ub6zij26

The article deals with the analysis of the fractional COVID‐19 epidemic model (FCEM) with a convex incidence rate. Keeping in view the fading memory and crossover behavior found in many biological phenomena, we study the coronavirus disease by using the noninteger Caputo derivative (CD). Under the Caputo operator (CO), existence and uniqueness for the solutions of the FCEM have been analyzed using fixed point theorems. We study all the basic properties and results including local and global stability. We show the global stability of disease‐free equilibrium using the method of Castillo‐Chavez, while for disease endemic, we use the method of geometrical approach. Sensitivity analysis is carried out to highlight the most sensitive parameters corresponding to basic reproduction number. Simulations are performed via first‐order convergent numerical technique to determine how changes in parameters affect the dynamical behavior of the system.

In December 2019, Wuhan City (China) gains worldwide attention, as an unknown virus starts infected human, resulting in a huge number of casualties coupled with unparalleled financial loss and social disruption. The virus was then recognized as a coronavirus (virus) and is then referred to as an extreme virus (SARS-CoV-2). The government tried their best to control the spread of the virus, but it took the world by storm and it got perforated to other countries by different means. Initial emergency control measures were put in place, in addition to the emphasis of the World Health Organization (WHO) on research funding to investigate its causes, get accurate diagnostics, and find a cure for the disease. These factors are still in the evolving stage, and very recently, vaccination efforts have born fruit after successful tests, which lead to the vaccination of chosen groups of population in a developed nation.

The duration of the virus is 11 to 14 days as per initial medical studies. When the comorbid disorders are noticed in patients, aged individuals are particularly susceptible to COVID-19. The rate of transmission of COVID-19 is very high, and between 2.2 and 3.58 is the basic reproduction number. That is why it is quickly spreading across the globe and impacting 213 countries. Thus, COVID-19 was announced as one of the global pandemics by WHO on January 30, 2020. 1 After China, the country with highly exposed individuals were Iran and Italy. Iran is Pakistan's neighboring country, and thousands of people visit Iran for religious and commercial reasons every year. Some of these individuals were infected while coming back from Iran, and they became the cause of the spread of the virus in Pakistan. In spite of the border closure, the first case in Pakistan was officially stated on February 26, 2020, after the infected person returned from Iran from Karachi. As a result, the government acted swiftly by imposing smart lock down quarantine of the infected persons at home. Mathematical modeling is regarded as an efficient method to explain the infection complex behavior. [2] [3] [4] [5] [6] For decades, mathematicians have often used mathematical modeling methods for various natural calamities. The study of infectious diseases for the purposes of disease planning and control has been vague. [7] [8] [9] [10] [11] [12] It is not practical to use a single model to characterize the actual process distribution. Increasing complexity, but at the expense of complicated mathematics, will make models more precise and important. Simple models such as the one described in the current study are also useful where there is no formal vaccine or careful monitoring of treatment. 13 The generalization of classical calculus is fractional calculus. Fractional derivative mathematical models provide more insights into the disease under consideration. [14] [15] [16] [17] [18] [19] In Srivastava et al 20 and Singh, 21 separate fractional operators with singular and nonsingular kernels were proposed. Latest literature and references to previous works [22] [23] [24] can be found in the implementations of these fractional operators. Very few COVID-19 models have recently been proposed based on fractional-order operators. Bonyah and Zarin 25 found a fractional derivative mathematical model for cancer and hepatitis co-dynamics and tested its effects. Currently, in this work, inspired by the above discussion, we work the dynamic analysis of the COVID-19 model given in Mandal et al. 26 The remainder of the organization of the paper is as follows: Section 2 introduces the model formulation for the fractional-order derivative. Section 3 includes, for the Caputo model, the presence and uniqueness of solutions, the basic reproductive number, the invariant region, and the points of equilibrium. Section 4 presents the local stability of the proposed model at the corona-free and corona-present equilibrium points of the fractional-order derivative. In Section 5, we review the global stability of the proposed model at the equilibrium points of corona free and corona current. The relevant parameters are calculated in Section 6, and their effect on the basic reproduction number is presented through sensitivity analysis. Section 7 describes the graphical outcomes by biological debate. The work has been summarized in Section 8 with valuable suggestions about the management of diseases in society.

Consider the recruitment constant rate to susceptible (S) class be A and to depict the rate of transmission. Now, when S mingles with an exposed (E) ones, the virus gets a significant increase; therefore, one can attain 1 (0 < 1 < 1) some part of S human with cautious precaution and 2 (0 < 2 < 1) some part of the E individual with cautious precaution for the transmission of disease.

The portion of the population class E with positive COVID-19 is also deemed to be contaminated, and they are believed to be hospitalized. Let the parts of the E class which moves to I and Q classes, respectively, be and b 2 . Let the cQ parts of communities move to the I level among the Q groups of populations, and the b 1 Q part changes to I.

Surmise and be the R rate of I and E hospitalized infected populations, respectively. Assume d is the normal death rate to all groups of populations and the death rate caused by COVID-19 is . It is also statistically observed that an individual who recovered from COVID-19, has a little risk of becoming infected again. We therefore presume that no part of R is going back to the disadvantaged class. We also recognize that the individuals of the hospitalized infected class will not transmit the disease or spread an eligible amount of disease in the formulation of a statistical model since they are separated entirely from the S individuals. Further, sufficient governmental initiatives were introduced by the numerous governmental and nongovernmental organizations to monitor the COVID-19 pandemic. This policy may therefore be regarded as one of the most powerful control methods, and this policy will primarily support the vulnerable population of COVID-19 cases.

Following the aforementioned assumptions and using the steps described, 26 the Caputo variant of the model is as comes next:

We consider

and

Summarizing that S(t) and S * (t) are couple functions, we reach

Taking into account

Continuing in the same way, one gets

Then, we get

Additionally, by using Equations (3.5) and (3.6) and considering that

(3.9)

Proof. The K 1 , K 2 , K 3 , K 4 , and K 5 satisfy the Lipchitz condition. Thus, we obtain

with r i = (M( )∕ )b i < 1 by hypothesis. Thus, uniform convergent is achieved. 27 Applying the theorem on limit in Equation (3.8) as n → ∞ affirms the solution of uniqueness (2.1). Hence, the unique solution for Equation (2.1) exists on the basis of (3.10).

Surmise that E 0 represents the disease-free equilibrium point for (2.1) so that

. (3.13) Let (E, Q, I) be the infected compartment. Then, we obtain

(3.14)

Then, we reach

) .

) . R 0 is given by

We define

Then, we get

When t → ∞, we reach N ≤ (b / ). We study (2.1) as

The endemic equilibrium point for (2.1) is denoted by E * = (S * , E * , Q * , I * , R * ) for which the disease is endemic in the population (i.e., at least one of E * , Q * , and I * is nonzero). The equations of (2.1) are reset to attain S * , E * , Q * , I * , and R * .

,

For R 0 > 1, the positivity of the above equilibrium point (3.15) is assured.

We establish the local stability of the (2.1) in this section at corona-free point E 0 as well as at corona-present equilibrium point E * . Proof. Jacobian matrix of the system (2.1) at E 0 is

. Then, we obtain

Therefore, the virus-free equilibrium point is asymptotically stable. 13

Now we study the local asymptotic stability of the endemic equilibrium E * .

Theorem 3. If R 0 > 1 and = 1, then the corona-present equilibrium E * of the system (2.1) is locally asymptotically stable.

Proof. Suppose R 0 > 1; then the existence of the endemic equilibrium point is assured. The Jacobian matrix of the system (2.1) at E * is

The two eigenvalues of J | * | are negative; that is, 1 = −d and 2 = −( + d + ). Additionally, we can have

)

.

Then, we get

It is observed that here A 1 , A 2 , A 3 and A 1 A 2 − A 3 are all positive for any parametric value. Therefore, (2.1) is locally asymptotically stable around its endemic equilibrium point E * .

We use the method of Castillo-Chavez et al 28 to establish the global stability for disease-free equilibrium, whereas for the global stability of endemic equilibrium, the generalization of Lyapunov theory 29 is used.

According to Castillo-Chavez et al, 28 decompose (2.1) into two subsystems and put the fractional parameter equal to one i.e

where 1 and 2 are the number of uninfected and infected individuals, respectively. Therefore, 1 = (S, R) ∈ R 2 and 2 = (E, Q, I) ∈ R 3 . For disease-free equilibrium, the existence of global stability depends on the following: 0) is an M-matrix having the positive off-diagonal entries and Δ represents the feasible region. Thus, the following statement holds. 

.

Then, we obtain

) .

For S = S 0 , R = R 0 , and G( 1 , 0) = 0 , we get

From Equation (5.4) as t → ∞, 1 → 0 1 . Thus, 1 = 0 1 is globally asymptotically stable. To show the second condition, that is, 

Thus,H( 1 , 2 ) is positive definite. Matrix B is given by

) .

As from the model (2.1), the total population is bounded by S 0 , that is, S, E, Q, I, R ≤ S 0 , so 2 ) is positive definite. Also from equation (5.8), it is clear that the matrix B is M-matrix that is the off diagonal element are non-negative. Thus, Conditions 1 and 2 are satisfied; so by Lemma 1, the disease-free equilibrium point E 0 is globally asymptotically stable.

According to Li and Muldowney, 29 take into account . x = (x), (5.9) where U ⊂ R n is connected simply and f : U → R n is a function that f ∈ C 1 (U). Surmising that (x * ) = 0 depict the solution of (5.9) and to x(t, x 0 ), we consider the following:

3. A compact absorbing set K ∈ U exists. 4. Equation (5.9) has a unique equilibrium. Now for asymptotically global and local stability, consider the following lemma.

Take into consideration the following:

.

(5.10) Equation (5.10) is a matrix valued function on U. Suppose that P −1 exist and is continuous for x ∈ K. Now describe a quantityq, such thatq x(s, x 0 ) ))]ds, (5.11) where J [2] is the second additive compound matrix of J; that is, J(x) = U (x) and B = P P −1 + PJ [2] P −1 . Let (B) be the Lozinski measure of the matrix B with respect to the norm ||.|| in R n 30 presented as follows:

is globally asymptotically stable in U, ifq < 0. 

(5.13)

The Jacobian matrix of system (5.13) is

The second additive compound matrix is

The function P( ) = P(S, E, Q) = diag

, and then

) .

,

Matrix B = PJ [2] P −1 + P P −1 can be expressed in the matrix form

where B 11 = −k 11 + k 22 ,

)

.

Suppose (k 1 , k 2 , k 3 ) is a vector in R 3 , with the norm ||.|| described by

and let (B) be the Lozinski measure with respect to this norm; we choose

where

and |B 12 |, |B 21 | are matrix norms with respect to the l 1 vector norm, and 1 refers to the Lozinski measure with respect to this l 1 norm; then

Hence, we have g 1 = −k 11 + k 22 ,

(5.14)

Thus, we have

(5.15) By using S * =(b 2 + + +d )

(1− 1 )(1− 2 ) , we get

Moreover, we get

which implies q = 0. Therefore, the Bendixson criterion is verified. We prove that positive equilibrium (S * , E * , Q * ) is globally asymptotically stable.

Examine the subsystem of system (2.1) The integrating factors for the system is e t( +d + ) and e t(d ) , which are used to solve the system. So for t → ∞, I → I * and R → R * , which is sufficient to prove that the endemic equilibrium point E * is globally asymptotically stable. Following the methods in Khan et al, 13 the elasticity index along with parameter values have been obtained as shown in Table 1 . The sensitivity of various parameters with R 0 is depicted in Figures 1-5. 

We present the numerical simulations in Figures 1-13 by using Table 2 . Consider a general Cauchy problem of fractional order having autonomous nature

where = (a, b, c, w) ∈ R 4 + is a real-valued continuous vector function which satisfies the Lipchitz criterion given as

where M is a positive real Lipchitz constant.

Using the fractional-order integral operators, one obtains where J Ω 0,t is the fractional-order integral operator in Riemann-Liouville. Consider an equispaced integration intervals over [0, T] with the fixed step size h (= 10 −2 for simulation) = T∕n, n ∈ N. Suppose that y p is the approximation of y(t) at t = t p for p = 0, 1, … n. Consequently, our model becomes 

The initial conditions and the parameters' values are used as described in the table above. The physical perspective of the model's individual state variables under the Caputo fractional operator is shown in Figures 6-10 . Since the symptomatic individuals are assumed to be more infectious than the asymptomatic, the transmission of COVID-19 through contacts in households, workplaces, schools, from foodstuffs, or during commute rises. This leads to a surge of the virus in environments such as workplace, school, foodstuffs, and public transport, and consequently, more cases of the coronavirus are confirmed. In Figures 11-15 , it can be noticed that the fractional-order = is varied for 1, 0.888, 0.666. One can easily see the robust nature of the Caputo operator than the integer variant of the model. For decreasing varying values of (disease transmission rate) as shown in Figures 16 and 17 , I(t) is also increasing. Similarly, in Figures 18 and 19 , the effect of mortality rate on I(t) has been shown. For increasing values of as shown in Figure 18 , an increasing pattern in I(t) is noticed. Similarly, for decreasing values of as shown in Figure 19 , a decreasing pattern in I(t) is noticed. Increasing-decreasing patterns are shown in this case. The model shows that quarantined policy is needed to avoid a large COVID-19 epidemic. There is a growing concern that this disease could continue to ravage the human population globally since many aspects of the COVID-19 are yet to be discovered, which also poses a challenge to the long-term mathematical modeling of the disease. 

The fractional derivative in Caputo sense is deployed to analyze the COVID-19 epidemic model. Local and global stability of the model along with the existence and positivity of the solution is studied. The sensitivity of different parameters corresponding to basic reproduction number has been highlighted graphically. The present finding shows that the effects of our model's predictions are in line with those reported in the literature and are therefore estimated to be more acceptable values of the parameters. At different stages, we varied the related parameters to their baseline values and the graphical results for the separate values of the fractional parameter of the Caputo derivative. The decrease in successful contacts (up to 35 %) and an increase in contact-tracing policy to quarantine the exposed individuals b 2 (up to 45%) to their baseline value have been found to significantly decrease the peak of contaminated curves. Furthermore, for smaller values of the fractional parameter , the reduction in the infected population becomes more important. Thus, we conclude that the findings obtained are accurate, practical, and more biologically feasible for the fractional scenario.

COVID-19 outbreak: current scenario of Pakistan

Modelling diseases with relapse and nonlinear incidence of infection: a multi group epidemic model

The stability analysis of an SVEIR model with continuous age-structure in the exposed and infection classes

Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An Introduction

A mathematical model of hepatitis B virus transmission and its application for vaccination strategy in China

Hepatitis B in a high prevalence New Zealand population: a mathematical model applied to infection control policy

Complete global stability for an SIRS epidemic model with generalized nonlinear incidence rate and vaccination

Global stability for the SEIR model in epidemiology

Stability analysis and optimal vaccination of an SIR epidemic model

Modeling the transmission dynamics and control of hepatitis B virus in China

Mathematical analysis of a cholera model with public health interventions

Dynamical behavior of a hepatitis B virus transmission model with vaccination

Stability analysis and optimal control of COVID-19 with convex incidence rate in Khyber Pakhtunkhawa (Pakistan)

Planar system-masses in an equilateral triangle: numerical study within fractional calculus

A new iterative method for the numerical solution of high-order non-linear fractional boundary value problems

A new adaptive synchronization and hyperchaos control of a biological snap oscillator

The fractional features of a harmonic oscillator with position-dependent mass

On the fractional optimal control problems with a general derivative operator

A hybrid functions numerical scheme for fractional optimal control problems: application to non-analytic dynamical systems

An efficient computational approach for a fractional-order biological population model with carrying capacity

Analysis of fractional blood alcohol model with composite fractional derivative

An efficient computational technique for local fractional Fokker Planck equation

A new analysis of fractional fish farm model associated with Mittag-Leffler type kernel

On the analysis of vibration equation involving a fractional derivative with Mittag-Leffler law

Mathematical modeling of cancer and hepatitis co-dynamics with non-local and non-singular kernal

A model based study on the dynamics of COVID-19: prediction and control

Introduction to Functional Analysis

Mathematical approaches for emerging and reemerging infectious diseases: an introduction

A geometric approach to global-stability problems

Logarithmic norms and projections applied to linear differential systems

This work was supported by Taif University Researchers Supporting Project number (TURSP-2020/154), Taif University, Taif, Saudi Arabia.

This work does not have any conflicts of interest.

There are no funders to report for this submission.

Rahat Zarin https://orcid.org/0000-0002-9759-2924 Mustafa Inc https://orcid.org/0000-0003-4996-8373