In this work, we present a model that uses the fractional order Caputo derivative for the novel Coronavirus disease 2019 (COVID-19) with different hospitalization strategies for severe and mild cases and incorporate an awareness program. We generalize the SEIR model of the spread of COVID-19 with a private focus on the transmissibility of people who are aware of the disease and follow preventative health measures and people who are ignorant of the disease and do not follow preventive health measures. Moreover, individuals with severe, mild symptoms and asymptomatically infected are also considered. The basic reproduction number (ℛ0) and local stability of the disease-free equilibrium (DFE) in terms of ℛ0 are investigated. Also, the uniqueness and existence of the solution are studied. Numerical simulations are performed by using some real values of parameters. Furthermore, the immunization of a sample of aware susceptible individuals in the proposed model to forecast the effect of the vaccination is also considered. Also, an investigation of the effect of public awareness on transmission dynamics is one of our aim in this work. Finally, a prediction about the evolution of COVID-19 in 1000 days is given. For the qualitative theory of the existence of a solution, we use some tools of nonlinear analysis, including Lipschitz criteria. Also, for the numerical interpretation, we use the Adams-Moulton-Bashforth procedure. All the numerical results are presented graphically.
COVID-19 is an infectious viral disease brought about by a virus known as severe acute syndrome Coronavirus [1]. This virus was first reported at the end of 2019 in Wuhan, China. Then, it was reported in other places in China surprisingly [2]. In three months, the disease was transmitted in many countries of the world. Therefore, in April 2020, the World Health Organization (WHO) announced it as an outbreak. Its symptoms involve breathing trouble, exhaustion, fever, dry hack, sleepiness, conjunctivitis, chest torment, loss of discourse, loose bowels, and painful throat. Furthermore, in serious cases, pneumonia, multiorgan letdown, intense respiratory pain condition, septic shock, abnormal heart rhythm, myocarditis, blood clumps, cardiovascular breakdown, encephalitis, stroke, and Guillain Barré disorder, reciprocal lung entrance have been reported [3,4]. Additionally, a few patients might experience the ill effects of looseness of the bowels, loss of craving, taste, or smell, with no indications of breathing problem [5–7].
Investigating infectious diseases through various procedures is an attractive area of research in recent times. One of the important areas in this regard is devoted to the mathematical modeling of infectious diseases. Mathematical modeling is a powerful tool to study infectious diseases for analysis of the dynamic of diseases and backing control systems [8]. Plenty of research work has been published by researchers. The mentioned work is devoted to investigating disease transmission dynamics and their control. In the same fashion to understand the dynamics of COVID-19 and the tracking down—of the reasonable result of the outbreak is helpful for public health initiatives and consciousness programs, for which significant work has been published, and we refer to a few such as [9–14].
Fractional calculus can more precisely describe natural phenomena as compared to ordinary calculus, where derivatives and integrations have integer-orders. Numerous scientists have given more attention to investigating the dynamics of fractional-order models regarding the COVID-19 pandemic [15–38]. There are many papers showing the effectiveness of fractional calculus (see [39–45]). Further, the dynamics of a stochastic epidemic model has been studied, for which we refer [46,47] and mathematical modeling of COVID-19 pandemic using the Caputo-Fabrizio fractional derivative (see [48]). Also, a model of COVID-19 disease has been investigated under the stochastic concept in [49].
In this research work, an SEIR-type model that concentrates on the transmission dynamics of the COVID-19 outbreak are considered. Here, with a special spotlight on the transmissibility of people with mild, severe, and without side effects like the existence of people who tested positive for sharp, mild, or asymptomatic manifestations, and separating infectious compartment into two fundamental compartments of hospitalized people with mild symptoms and those in dense care units are involved in our proposed model.
This work is organized as: Part 2 is devoted to the formulation of the proposed model for COVID-19. In Parts 3 and 4, we analyze the proposed model with fractional order. Some qualitative analysis, computation of ℛ0, and local stability for the disease-free equilibrium (DFE) in terms of ℛ0, of the proposed model are obtained in Part 5. The stability analysis is studied in Part 6. Part 7 is dedicated to data fit and explanation of our model via numerical simulation and comparison with real data for Morocco. The conclusion is given in the last section.
The Proposed Epidemic Model
Here in this section, the proposed model is formulated. The flowchart of the model is given in Fig. 1 to understand the evolution of the mentioned disease. Let the total human population at time ϑ be denoted by N(ϑ). The proposed model divides the human population at time ϑ into eleven compartments which are described in Table 1. We generalize the SEIR model of the spread of COVID-19 with a specific focus on the transmissibility of people who are Individuals who follow precautionary health measures Sa and Individuals who do not follow precautionary health measures Su. We separate hospitalized cases into severe Hs and mild cases Hm and incorporate awareness programs. Moreover, individuals with severe Iss, mild symptoms Ims and asymptomatically infected Ia are also considered. Furthermore, the immunization of a sample of aware susceptible individuals in the proposed model to forecast the effect of vaccination V is also considered.
Flowchart of model <xref ref-type="disp-formula" rid="eqn-1">(1)</xref>
Compartments symbols and description
Symbol
Description
Sa
Individuals follow precautionary health measures
Su
Individuals do not follow precautionary health measures
E
Exposed class
Iss
Severe symptoms infectious individuals
Ims
Mild severe symptoms infectious individuals
Ia
Infectious but asymptomatic individuals
Hm
Hospitalized class
Hs
Hospitalized individuals with sharp symptoms
Icu
Intensive cure unit class
R
Recovery with immunity class
D
Dead class
θ
The rate of conversion from exposed to infectious
β
The human-to-human transmission coefficient of each unit
p1
The rate from exposed class to symptomatic infectious includes Iss
p2
The rate of exposed class to infectious class Ims
1−p1−p2
Rate of exposed persons to asymptomatic class Ia
h
The rate at which a person leaves the class Iss
τ
Hospitalization rates from Iss to Hm
ρ
Hospitalization rates from Iss to Hs
ηk,k=a,i,Hm,Hs
Recovery rates
δi,δHs
COVID-19 induced death rates
ψI1
Rate at which Iss to Ims
φI1
Rate at which Ims to Iss
ψI2
The rate of hospitalized class from temperate to sharp isolation
φI2
The rate of hospitalized class from sharp to temperate isolation
Λ
Requitement rate
Thus, we have
N(ϑ)=Sa(ϑ)+Su(ϑ)+E(ϑ)+Iss(ϑ)+Ims(ϑ)+Ia(ϑ)+Hm(ϑ)+Hs(ϑ)+Icu(ϑ)+R(ϑ)+D(ϑ).
From Fig. 1, we formulate our proposed model as
{Sa′(ϑ)=Λ−βSa(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))+vSu(ϑ),Su′(ϑ)=−βSu(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−vSu(ϑ),E′(ϑ)=βSa(ϑ)+Su(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−θE(ϑ),Iss′(ϑ)=θp1E(ϑ)+φ1Ims(ϑ)−(ψ1+h)Iss(ϑ),Ims′(ϑ)=θp2E(ϑ)+ψ1Iss(ϑ)−(ηi+φ1)Ims(ϑ),Ia′(ϑ)=θ(1−p1−p2)E(ϑ)−ηaIa(ϑ),Hm′(ϑ)=hτIss(ϑ)+ψ2Hs(ϑ)−(φ2+ηHm)Hm(ϑ),Hs′(ϑ)=hρIss(ϑ)+φ2Hm(ϑ)−(ψ2+ηHs+δHs)Hs(ϑ),Icu′(ϑ)=h(1−τ−ρ)Iss(ϑ)−Icu(ϑ),R′(ϑ)=ηiIms(ϑ)+ηaIa(ϑ)+ηHsHs(ϑ)+ηHmHm(ϑ)+(1−δi)Icu(ϑ),D′(ϑ)=δHsHs(ϑ)+δiIcu(ϑ),with initial conditions as
{Sa(0)=Sa,0,Su(0)=Su,0,E(0)=E0,Iss(0)=Iss,0,Ims(0)=Ims,0,Ia(0)=Ia,0,Hm(0)=Hm,0,Hs(0)=Hs,0,Icu(0)=Icu,0,D(0)=D0,R(0)=R0.
Fractional Order Model
In this part, we recall some definitions from fractional calculus which are needed throughout the paper.
Definition 3.1. [50] The Caputo fractional derivative of f of order μ is described as
f(μ)(ϑ)=1Γ(n−μ)∫0ϑf(n)(s)(ϑ−s)μ−n+1ds,where n−1<μ≤n and Γ(x)=∫0∞e−ϑϑx−1dϑ, which is called the Euler’s Gamma function.
The fractional integral [51] having of order μ>0 is defined by
Iμ(f(ϑ))=1Γ(μ)∫0ϑf(s)(ϑ−s)1−μds,where ϑ>0, satisfies:
{(Iμf(ϑ))(μ)=f(ϑ),Iμ(f(ϑ)(μ))=f(ϑ)−f(0).
Now, let us present the fractional order model of the COVID-19 by means of the Caputo derivative as follows:
{Sa(μ)(ϑ)=Λ−βSa(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))+vSu(ϑ),Su(μ)(ϑ)=−βSu(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−vSu(ϑ),E(μ)(ϑ)=βSa(ϑ)+Su(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−θE(ϑ),Iss(μ)(ϑ)=θp1E(ϑ)+φ1Ims(ϑ)−(ψ1+h)Iss(ϑ),Ims(μ)(ϑ)=θp2E(ϑ)+ψ1Iss(ϑ)−(ηi+φ1)Ims(ϑ),Ia(μ)(ϑ)=θ(1−p1−p2)E(ϑ)−ηaIa(ϑ),Hm(μ)(ϑ)=hτIss(ϑ)+ψ2Hs(ϑ)−(φ2+ηHm)Hm(ϑ),Hs(μ)(ϑ)=hρIss(ϑ)+φ2Hm(ϑ)−(ψ2+ηHs+δHs)Hs(ϑ),Icu(μ)(ϑ)=h(1−τ−ρ)Iss(ϑ)−Icu(ϑ),R(μ)(ϑ)=ηiIms(ϑ)+ηaIa(ϑ)+ηHsHs(ϑ)+ηHmHm(ϑ)+(1−δi)Icu(ϑ),D(μ)(ϑ)=δHsHs(ϑ)+δiIcu(ϑ).
Existence Theory
In this part, we study the existence of a solution for the fractional order model Eq. (4). The existence theory is an important consequence of applied analysis and has many applications. It provides us with whether the model we study exists or not and if it exists whether it has a solution or not. If a problem has a solution, whether it will be unique or multiple. It also helps us in investigating the qualitative behavior of the solution to the problem.
Consider the following subsystem:
{Su(μ)(ϑ)=−βSu(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−vSu(ϑ),E(μ)(ϑ)=βSa(ϑ)+Su(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−θE(ϑ),Iss(μ)(ϑ)=θp1E(ϑ)+φ1Ims(ϑ)−(ψ1+h)Iss(ϑ),Ims(μ)(ϑ)=θp2E(ϑ)+ψ1Iss(ϑ)−(ηi+φ1)Ims(ϑ),Ia(μ)(ϑ)=θ(1−p1−p2)E(ϑ)−ηaIa(ϑ),Hm(μ)(ϑ)=hτIss(ϑ)+ψ2Hs(ϑ)−(φ2+ηHm)Hm(ϑ),Hs(μ)(ϑ)=hρIss(ϑ)+φ2Hm(ϑ)−(ψ2+ηHs+δHs)Hs(ϑ),Icu(μ)(ϑ)=h(1−τ−ρ)Iss(ϑ)−Icu(ϑ).
Since
{Sa(μ)(ϑ)=Λ−βSa(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))+vSu(ϑ),R(μ)(ϑ)=ηiIms(ϑ)+ηaIa(ϑ)+ηHsHs(ϑ)+ηHmHm(ϑ)+(1−δi)Icu(ϑ),D(μ)(ϑ)=δHsHs(ϑ)+δiIcu(ϑ).
Therefore, to study the existence and uniqueness of solutions of the system Eq. (4), we focus to study Eq. (5).
Assume that
{K1(ϑ,Su)=−βSuN(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−vSu,K2(ϑ,E)=βSa(ϑ)+Su(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−θE,K3(ϑ,Iss)=θp1E(ϑ)+φ1Ims(ϑ)−(ψ1+h)Iss,K4(ϑ,Ims)=θp2E(ϑ)+ψ1Iss(ϑ)−(ηi+φ1)Ims,K5(ϑ,Ia)=θ(1−p1−p2)E(ϑ)−ηaIa,K6(ϑ,Hm)=hτIss(ϑ)+ψ2Hs(ϑ)−(φ2+ηHm)Hm,K7(ϑ,Hs)=hρIss(ϑ)+φ2Hm(ϑ)−(ψ2+ηHs+δHs)HsK8(ϑ,Icu)=h(1−τ−ρ)Iss−Icu.
Applying the fractional integral to Eq. (5), it follows from initial conditions that
{Su(ϑ)=Su,0+1Γ(μ)∫0ϑ(ϑ−s)μ−1K1(s,Su(s))ds,E(ϑ)=E0+1Γ(μ)∫0ϑ(ϑ−s)μ−1K2(s,E(s))ds,Iss(ϑ)=Iss,0+1Γ(μ)∫0ϑ(ϑ−s)μ−1K3(s,Iss(s))ds,Ims(ϑ)=Ims,0+1Γ(μ)∫0ϑ(ϑ−s)μ−1K4(s,Ims(s))ds,Ia(ϑ)=Ia,0+1Γ(μ)∫0ϑ(ϑ−s)μ−1K5(s,Ia(s))ds,Hm(ϑ)=Hm,0+1Γ(μ)∫0ϑ(ϑ−s)μ−1K6(s,Hm(s))ds,Hs(ϑ)=Hs,0+1Γ(μ∫0ϑ(ϑ−s)μ−1K7(s,Hs(s))ds,Icu(ϑ)=Icu,0+1Γ(μ)∫0ϑ(ϑ−s)μ−1K8(s,Icu(s))ds.
Assume that ‖E(ϑ)‖≤c1, ‖Iss(ϑ‖≤c2, ‖Ims(ϑ)‖≤c3, where ci,i=1,…,3, are some positive constants. Denote L1=βN(c1+c2+c3)+v,L2=θ,L3=ψ1+h,L4=φ1+ηi,L5=ηa,L6=φ1+ηHm,L7=ψ2+ηHs+δHs.
Theorem 4.1. The kernels Ki(i=1,…,8), satisfy the Lipschitz and contractive conditions if
0≤Li<1hold.
Proof. Let S∗ and S∗∗ be two functions, then
‖K1(ϑ,S∗)−K1(ϑ,S∗∗)‖=‖(−βN(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−v)(S∗∗−S∗)‖≤[βN(‖Iss‖+‖Ims‖+‖Ia‖)+v]‖S∗(ϑ)−S∗∗(ϑ)‖≤L1‖S∗(ϑ)−S∗∗(ϑ)‖.
Thus, Ki satisfy Lipschitz condition, for i=1,….,8. From Eq. (9), Ki (i=1,…,8) are contraction.
Theorem 4.2. The model Eq. (4) has a unique solution provided
1μΓ(μ)rμLi<1, for i=1,2,…,8,when r∈[0,ϑ].
Proof. The proof for the above result is the same as the proof of Theorem 1 in [52].
Reproductive Number and Stability AnalysisReproductive Number
We follow the references [53,54] to calculate ℛ0. In order to calculate ℛ0, we consider the subsystem Eq. (5).
The DFE point is given by ℰ0=(N,0,0,0,0,0,0,0,0,0,0). From the system (4), we can write
F=(0βSa+SuN(Iss+Ims+Ia)000000),and V associated with the net rate outside of the corresponding compartments
v=(−βSuN(Iss+Ims+Ia)−vSu−θEθp1E+φ1Ims−(ψ1+h)Issθp2E+ψ1Iss−(ηi+φ1)Imsθ(1−p1−p2)E−ηaIahτIss+ψ2Hs−(φ2+ηHm)HmhρIss+φ2Hm−(ψ2+ηHs+δHs)Hsh(1−τ−ρ)Iss−Icu)
It follows that the Jacobian matrices of F and V at ℰ0 respectively are
ℱ=(0000000000βββ000000000000000000000000000000000000000000000000000),and𝒱=(−v00000000−θ0000000θp1−(h+ψ1)φ100000θp2ψ1−(ηi+φ1)00000θ(1−p1−p2)00−ηa00000hτ00−(φ2+ηHm)ψ2000hρ00φ2−(ψ2+ηHm+δHs)000h(1−τ−ρ)0000−1).
Thus, the ℛ0 is obtained as the spectral radius of −ℱ𝒱−1 as
ℛ0=maxλ∈Sp(−ℱ⋅𝒱−1)|λ|=Sp(−ℱ⋅𝒱−1)=β(1−p1−p2)ηa−β(p2(h+ψ1)+p1ψ1+p1(ηi+φ1)+p2φ1)hηi+hφ1+ψ1ηi.
Stability Analysis
The reported result below is given in Theorem 2 of [54].
Theorem 5.1. The DFE ℰ0 of model Eq. (4), is locally asymptotically stable if ℛ0<1 and unstable if ℛ0>1.
Fractional Order Model with Vaccination
In this part, we discuss the fractional order model with vaccination which is overseeing and controlling contagious diseases by giving security to powerless and susceptible individuals. We suppose that a specified proportion p of people in the mindful susceptible class are vaccinated. For this situation, inoculated people are moved to another compartment V(ϑ). Due to the vaccine does not supply impunity to all vaccine recipients, vaccinated people might become infected yet a lower rate than unvaccinated. In this situation, let σv∈[0,1] such that (1−σv) be the vaccine efficacy. The flowchart of the fractional order model is as given in Fig. 2.
Flowchart of the fractional model <xref ref-type="disp-formula" rid="eqn-10">(10)</xref>
We can interpret this diagram into fractional differential equations as follows:
{Sa(μ)(ϑ)=Λ−pSa(ϑ)−βSa(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))+vSu(ϑ),Su(μ)(ϑ)=−βSu(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−vSu(ϑ),E(μ)(ϑ)=βSa(ϑ)+Su(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))+σVβV(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ))−θE(ϑ),Iss(μ)(ϑ)=θp1E(ϑ)+φ1Ims(ϑ)−(ψ1+h)Iss(ϑ),Ims(μ)(ϑ)=θp2E(ϑ)+ψ1Iss(ϑ)−(ηi+φ1)Ims(ϑ),Ia(μ)(ϑ)=θ(1−p1−p2)E(ϑ)−ηaIa(ϑ),Hm(μ)(ϑ)=hτIss(ϑ)+ψ2Hs(ϑ)−(φ2+ηHm)Hm(ϑ),Hs(μ)(ϑ)=hρIss(ϑ)+φ2Hm(ϑ)−(ψ2+ηHs+δHs)Hs(ϑ),Icu(μ)(ϑ)=h(1−τ−ρ)Iss(ϑ)−Icu(ϑ),R(μ)(ϑ)=ηiIms(ϑ)+ηaIa(ϑ)+ηHsHs(ϑ)+ηHmHm(ϑ)+(1−δi)Icu(ϑ),D(μ)(ϑ)=δHsHs(ϑ)+δiIcu(ϑ),V(μ)(ϑ)=pSa(ϑ)−σVβV(ϑ)N(Iss(ϑ)+Ims(ϑ)+Ia(ϑ)),with Sa(0)=(1−p)Sa,0,Su(0)=(1−p)Su,0,E(0)=E0,Iss(0)=Iss,0,Ims(0)=Ims,0,Ia(0)=Ia,0, Hm(0)=Hm,0,Hs(0)=Hs,0,Icu(0)=Icu,0,D(0)=D0,R(0)=R0,V(0)=pN. Then, the fractional order model (10) has one equilibrium point
ℰv=((1−p)N,0,0,0,0,0,0,0,0,0,0,pN).
The components of infection in this model are, Su,E,Iss,Ims,Ia,Hm,Hs, and Icu.
Theorem 6.1. The reproduction number ℛv for the vaccinated fractional order model (10) is given by
ℛv=(1−(1−σv)p)ℛ0.
Proof. Let F and V be given by
F=(0βSa+SuN(Iss+Ims+Ia)+βσvVN(Iss+Ims+Ia)000000),and
V=(−βSuN(Iss+Ims+Ia)−vSu−θEθp1E+φI1Ims−(ψI1+h)Issθp2E+ψI1Iss−(ηi+φI1)Imsθ(1−p1−p2)E−ηaIahτIss+ψI2Hs−(φI2+ηHm)HmhρIss+φI2Hm−(ψI2+ηHs+δHs)Hsh(1−τ−ρ)Iss−Icu).
By the differential of F and V with respect to E,Iss,Ims,Ia,Hm, and Icu and evaluating at the ℰv, respectively, we get
ℱ:=JF(ℰv)=(1−(1−σv)p)(0000000000βββ000000000000000000000000000000000000000000000000000)and 𝒱:=JV(ℰv)=(−v00000000−θ0000000θp1−(h+ψI1)φI100000θp2ψI1−(ηi+φI1)00000θ(1−p1−p2)00−ηa00000hτ00−(φI2+ηHm)ψI2000hρ00φI2−(ψI2+ηHm+δHs)000h(1−τ−ρ)0000−1).
Thus, the ℛv for the vaccinated fractional order model is:
ℛv=β(1−(1−σv)p)((1−p1−p2)ηa−(p2(h+ψI1)+p1ψI1+p1(ηi+φI1)+p2φI1)hηi+hφI1+ψI1ηi)=(1−(1−σv)p)ℛ0.
Theorem 6.2. The DFE ℰv of model Eq. (10) will be asymptotically stable if ℛv<1, and unstable if ℛv>1.
Proof. The proof is given in the Theorem 2 of [54].
The critical percentage of the population pc needed to fulfill herd immunity is
pc=(11−σv)(1−1ℛ0).
pc represents the rate for which ℛ0 under the vaccination ℛv is equal to 1.
Numerical Simulations
In the numerical analysis of traditional differential equations, various numerical methods have been established. Among various numerical procedures, Adams’s methods represent one of the most used and studied classes of implicit (Adams-Moulton) and explicit (Adams-Bashforth) linear multistep methods. The aforementioned methods have gotten much popularity among researchers due to their good stability properties, reasonable computational cost, and simplicity of implementation. Therefore, several authors have extended the aforementioned methods. In particular, the mentioned procedures have also been extended to deal with fractional order differential equations numerically. A predictor-corrector algorithm has been extended to fractional differential equations with the corresponding Adams-Moulton-Bashforth procedure for traditional differential equations. The numerical method used to solve the model is the predictor-corrector method for fractional differential equations in [55]. Consider, the problem with fractional order
{y(μ)(ϑ)=f(ϑ,y(ϑ)),0<μ≤1,y(ϑ0)=y0,where f:[ϑ0,T]×D→R,D⊂R. The proposed scheme for (11) can be written as
{ynp=y0+hμ∑j=0n−1Λn−j−1fjyn=y0+hμκn,0f0+hμ∑j=0n−1Δn−jfj+hμΔ0f(ϑn,ynp),where Λn=(n+1)μ−nμΓ(μ+1),κn,0=(n−1)μ−1−nμ(n−μ−1)Γ(μ+2) and
Δn={1Γ(μ+2),n=0,(n−1)μ+1−2nμ+1+(n+1)μ+1Γ(μ+2),n=1,2,…
We extend the afore given scheme in Eq. (12) for our considered model Eq. (4) to simulate the results. The numerical simulations for fractional dynamic of the considered model Eqs. (4) and (10) are performed as the case study of Morocco. The parameters values are given in the Table 2 for the numerical interpretation. Also we compute the value of the reproduction number as ℛ0=1.0575 and ℛv=1.0574.
Values of the model parameters
Name
Value
Name
Value
Name
Value
Name
Value
Sa,0
699205 [56]
R0
0 [56]
Ia,0
3 [56]
Iss,0
100 [56]
Su,0
3000000 [56]
D0
0 [56]
Ims,0
100 [56]
Hs,0
0 (Estimated)
E0
185 [56]
Icu,0
37 [56]
Hm,0
5 (Estimated)
β
0.45 [57]
θ
1/5.1 [6]
ηa
0.4 [58]
ηi
0.4 [58]
τ
0.1259 [59]
ρ
0.13266 [59]
δi
0.6 [56]
ηHm
0.11624 [59]
ηHs
0.155 [59]
δHs
1−ηHm−ηHs
ψI1
0.169055 [59]
ϕI1
0.0341 [59]
ψI2
0.169055 (Estimated)
ϕI2
0.0041 (Estimated)
v
0.0173768 [59]
Λ
0.008 N (Estimated)
h
0.5 [56]
p1
0.3 [56]
p2
0.4 [56]
(1−σv)
0.9 [56]
p
4/1000 (Estimated)
In Fig. 3, the confirmed infected and death cases are described daily from July 01, 2020, to January 01, 2021, which is a total of 185 days. We implement numerical simulations to compare the results of our model with the real data in Fig. 3. The forecasted evolution of the outbreak of COVID-19 without and with vaccination in Morocco can be seen in Figs. 4 and 5, respectively.
Evolution of COVID-19 in Morocco per day
The evolution of the epidemic predicted by the fractional model <xref ref-type="disp-formula" rid="eqn-4">(4)</xref> with <inline-formula id="ieqn-143"><mml:math id="mml-ieqn-143"><mml:mi>μ</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>
The evolution of the epidemic predicted by the fractional model vaccinated with <inline-formula id="ieqn-144"><mml:math id="mml-ieqn-144"><mml:mi>μ</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula>
From the analysis of the acquired graphs, we observe that the population of the infected class decreases significantly by decreasing the fractional order of the derivative and the endemic state of the disease goes to the disease-free state in all the cases deliberated. Besides, we demonstrate that our fractional order model well describes the actual data of daily confirmed, recovered, and death cases. In Fig. 5, we note that the plot is much flatter than that in Fig. 4. As well, the curve for asymptomatic individuals is almost identical to the x-axis, and this indicates the importance of the vaccination approach to overcome the pandemic.
In Fig. 6, we give the graphical results of the suggested model (4) to analyze the effect of the fractional order.
The evolution of the epidemic predicted by the model <xref ref-type="disp-formula" rid="eqn-4">(4)</xref>
In Fig. 7, we give the graphical results of the suggested vaccination model (10) to analyze the influence of the fractional order.
Epidemic evolution predicted by the fractional vaccination model <xref ref-type="disp-formula" rid="eqn-10">(10)</xref>
We represent the spread of infection without vaccination in Fig. 8 and with vaccination in Fig. 9, for a period of 1000 days.
Epidemic evolution predicted by the fractional model
Epidemic evolution predicted by the fractional vaccination model
Conclusion
In this manuscript, we have given an analysis of the fractional order model of COVID-19, with theoretical calculations. We have established sufficient conditions for the existence theory using some tools of nonlinear analysis. The existence theory is an important area of research in recent times. Also, we have computed the basic reproductive number ℛ0. From the number mentioned above, we have predicted the proposed model’s stability. Further, we have extended the predictor-corrector method to numerically interpret the proposed model. We have given several graphical presentations using the numerical procedure. We have used some real data. It should be kept in mind that the model under consideration has described the evolution of COVID-19 in Morocco. We have given a graphical presentation of approximation solutions of various compartments in Figs. 3 and 4) and a simulation of this model under vaccination in Fig. 5). On another side, we have compared the numerical simulations with various values of the fractional-order μ. Moreover, our model gives a prediction about the evolution of COVID-19 in the next 1000 days. It seems from Figs. 6, and 7 that we must commit to preventive measures with the necessity of vaccination and adaptation to control and reduce the infection from further spreading in communities. In the future, the aforementioned analysis can be established for a more complex dynamical system of COVID-19 and the effect of vaccination on its reduction in the community. In our future work, we intend to study the model using the newly generalized partial derivative (see [60]), which overcomes some of the problems associated with conformational derivatives and some other fractional order derivatives.
The authors Kamal Shah, and Thabet Abdeljawad would like to thank Prince Sultan University for the support through the TAS Research Lab.
Funding Statement
The authors Kamal Shah, and Thabet Abdeljawad would like to thank Prince Sultan University for paying the APC.
Conflicts of Interest
The authors declare that they have no conflicts of interest to report regarding the present study.
ReferencesLai, C. C., Shih, T. P., Ko, W. C., Tang, H. J., Hsueh, P. R. (2020). Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges. International Journal of Antimicrobial Agents,55(3),105924. DOI 10.1016/j.ijantimicag.2020.105924.Baba, I. A., Ahmed, I., Al-Mdallal, Q. M., Jarad, F., Yunusa, S. (2022). Numerical and theoretical analysis of an awareness COVID-19 epidemic model via generalized Atangana-Baleanu fractional derivative. Journal of Applied Mathematics and Computational Mechanics,21(1),7–18. DOI 10.17512/jamcm.2022.1.01.Bedford, J., Enria, D., Giesecke, J., Heymann, D. L., Ihekweazu, C.et al. (2020). COVID-19: Towards controlling of a pandemic. The Lancet,395(10229),1015–1018. DOI 10.1016/S0140-6736(20)30673-5.Memon, Z., Qureshi, S., Memon, B. R. (2021). Assessing the role of quarantine and isolation as control strategies for COVID-19 outbreak: A case study. Chaos, Solitons & Fractals,144,110655. DOI 10.1016/j.chaos.2021.110655.Guo, Y. R., Cao, Q. D., Hong, Z. S., Tan, Y. Y., Chen, S. D.et al. (2020). The origin, transmission and clinical therapies on coronavirus disease 2019 (COVID-19) outbreak–an update on the status. Military Medical Research,7(1),1–10. DOI 10.1186/s40779-020-00240-0.Lauer, S. A., Grantz, K. H., Bi, Q., Jones, F. K., Zheng, Q.et al. (2020). The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: Estimation and application. Annals of Internal Medicine,172(9),577–582. DOI 10.7326/M20-0504.Liu, J., Liao, X., Qian, S., Yuan, J., Wang, F.et al. (2020). Community transmission of severe acute respiratory syndrome coronavirus 2, Shenzhen, China, 2020. Emerging Infectious Diseases,26(6),1320–1323. DOI 10.3201/eid2606.200239.Momoh, A. A., Ibrahim, M. O., Uwanta, I. J., Manga, S. B. (2013). Mathematical model for control of measles epidemiology. International Journal of Pure and Applied Mathematics,87(5),707–718. DOI 10.12732/ijpam.v87i5.4.Wu, P., Hao, X., Lau, E. H., Wong, J. Y., Leung, K. S.et al. (2020). Real-time tentative assessment of the epidemiological characteristics of novel coronavirus infections in Wuhan, China, as at 22 January 2020. Eurosurveillance,25(3),2000044. DOI 10.2807/1560-7917.ES.2020.25.3.2000044.Tang, B., Wang, X., Li, Q., Bragazzi, N. L., Tang, S.et al. (2020). Estimation of the transmission risk of the 2019-nCoV and its implication for public health interventions. Journal of Clinical Medicine,9(2),462. DOI 10.3390/jcm9020462.Peter, O. J., Qureshi, S., Yusuf, A., Al-Shomrani, M., Idowu, A. A. (2021). A new mathematical model of COVID-19 using real data from Pakistan. Results in Physics,24,104098. DOI 10.1016/j.rinp.2021.104098.Abdullahi, I., Baba, B. A., Nasidi, B. (2021). Optimal control model for the transmission of novel COVID-19. Computers, Materials, & Continua,66(3),3089–3106. DOI 10.32604/cmc.2021.012301.Mahmoudi, M. R., Baleanu, D., Band, S. S., Mosavi, A. (2021). Factor analysis approach to classify COVID-19 datasets in several regions. Results in Physics,25,104071. DOI 10.1016/j.rinp.2021.104071.Qureshi, S., Yusuf, A. (2019). Mathematical modeling for the impacts of deforestation on wildlife species using caputo differential operator. Chaos, Solitons & Fractals,126,32–40. DOI 10.1016/j.chaos.2019.05.037.Sindhu, T. N., Shafiq, A., Al-Mdallal, Q. M. (2021). On the analysis of number of deaths due to COVID-19 outbreak data using a new class of distributions. Results in Physics,21,103747. DOI 10.1016/j.rinp.2020.103747.Khan, M. A., Atangana, A. (2020). Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative. Alexandria Engineering Journal,59(4),2379–2389. DOI 10.1016/j.aej.2020.02.033.Shafiq, A., Lone, S. A., Sindhu, T. N., El Khatib, Y., Al-Mdallal, Q. M.et al. (2021). A new modified Kies Fréchet distribution: Applications of mortality rate of COVID-19. Results in Physics,28,104638. DOI 10.1016/j.rinp.2021.104638.Thabet, S., Abdo, M. S., Shah, K. (2021). Theoretical and numerical analysis for transmission dynamics of COVID-19 mathematical model involving Caputo-Fabrizio derivative. Advances in Difference Equations,2021(1),1–17. DOI 10.1186/s13662-021-03316-w.Abdulwasaa, M. A., Abdo, M. S., Shah, K., Nofal, T. A., Panchal, S. K.et al. (2021). Fractal-fractional mathematical modeling and forecasting of new cases and deaths of COVID-19 epidemic outbreaks in India. Results in Physics,20,103702. DOI 10.1016/j.rinp.2020.103702.Redhwan, S. S., Abdo, M. S., Shah, K., Abdeljawad, T., Dawood, S.et al. (2020). Mathematical modeling for the outbreak of the coronavirus (COVID-19) under fractional nonlocal operator. Results in Physics,19,103610. DOI 10.1016/j.rinp.2020.103610.Thabet, S. T., Abdo, M. S., Shah, K., Abdeljawad, T. (2020). Study of transmission dynamics of COVID-19 mathematical model under ABC fractional order derivative. Results in Physics,19,103507. DOI 10.1016/j.rinp.2020.103507.Sindhu, T. N., Shafiq, A., Al-Mdallal, Q. M. (2021). Exponentiated transformation of gumbel type-II distribution for modeling COVID-19 data. Alexandria Engineering Journal,60(1),671–689. DOI 10.1016/j.aej.2020.09.060.Akgül, A., Ahmed, N., Raza, A., Iqbal, Z., Rafiq, M.et al. (2021). New applications related to COVID-19. Results in Physics,20,103663. DOI 10.1016/j.rinp.2020.103663.Gao, W., Veeresha, P., Baskonus, H. M., Prakasha, D. G., Kumar, P. (2020). A new study of unreported cases of 2019-nCOV epidemic outbreaks. Chaos, Solitons & Fractals,138,109929. DOI 10.1016/j.chaos.2020.109929.Singh, H., Srivastava, H. M., Hammouch, Z., Nisar, K. S. (2021). Numerical simulation and stability analysis for the fractional-order dynamics of COVID-19. Results in Physics,20,103722. DOI 10.1016/j.rinp.2020.103722.Alshomrani, A. S., Ullah, M. Z., Baleanu, D. (2021). Caputo SIR model for COVID-19 under optimized fractional order. Advances in Difference Equations,2021(1),1–17. DOI 10.1186/s13662-021-03345-5.Habenom, H., Aychluh, M., Suthar, D. L., Al-Mdallal, Q., Purohit, S. D. (2022). Modeling and analysis on the transmission of COVID-19 pandemic in Ethiopia. Alexandria Engineering Journal,61(7),5323–5342. DOI 10.1016/j.aej.2021.10.054.Gao, W., Baskonus, H. M., Shi, L. (2020). New investigation of bats-hosts-reservoir-people coronavirus model and application to 2019-nCoV system. Advances in Difference Equations,2020(1),1–11. DOI 10.1186/s13662-020-02831-6.Naik, P. A., Yavuz, M., Qureshi, S., Zu, J., Townley, S. (2020). Modeling and analysis of COVID-19 epidemics with treatment in fractional derivatives using real data from Pakistan. The European Physical Journal Plus,135(10),1–42. DOI 10.1140/epjp/s13360-020-00819-5.Srivastava, H. M., Jan, R., Jan, A., Deebani, W., Shutaywi, M. (2021). Fractional-calculus analysis of the transmission dynamics of the dengue infection. Chaos: An Interdisciplinary Journal of Nonlinear Science,31(5),053130. DOI 10.1063/5.0050452.Boulaaras, S., Jan, R., Khan, A., Ahsan, M. (2022). Dynamical analysis of the transmission of dengue fever via Caputo-Fabrizio fractional derivative. Chaos, Solitons & Fractals: X,8,100072. DOI 10.1016/j.csfx.2022.100072.Baleanu, D., Fernandez, A. (2019). On fractional operators and their classifications. Mathematics,7(9),830. DOI 10.3390/math7090830.Gao, W., Veeresha, P., Prakasha, D. G., Baskonus, H. M. (2020). Novel dynamic structures of 2019-nCoV with nonlocal operator via powerful computational technique. Biology,9(5),107. DOI 10.3390/biology9050107.Jajarmi, A., Ghanbari, B., Baleanu, D. (2019). A new and efficient numerical method for the fractional modeling and optimal control of diabetes and tuberculosis co-existence. Chaos: An Interdisciplinary Journal of Nonlinear Science,29(9),093111. DOI 10.1063/1.5112177.Acay, B., Bas, E., Abdeljawad, T. (2020). Fractional economic models based on market equilibrium in the frame of different type kernels. Chaos, Solitons & Fractals,130,109438. DOI 10.1016/j.chaos.2019.109438.Uçar, S., Uçar, E., Özdemir, N., Hammouch, Z. (2019). Mathematical analysis and numerical simulation for a smoking model with Atangana-Baleanu derivative. Chaos, Solitons & Fractals,118,300–306. DOI 10.1016/j.chaos.2018.12.003.Khan, M. A., Atangana, A., Alzahrani, E. (2020). The dynamics of COVID-19 with quarantined and isolation. Advances in Difference Equations,2020(1),1–22. DOI 10.1186/s13662-020-02882-9.Wrapp, D., Wang, N., Corbett, K. S., Goldsmith, J. A., Hsieh, C. L.et al. (2020). Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science,367(6483),1260–1263. DOI 10.1126/science.abb2507.Kumar, P., Erturk, V. S., Nisar, K. S., Jamshed, W., Mohamed, M. S. (2022). Fractional dynamics of 2019-nCOV in Spain at different transmission rate with an idea of optimal control problem formulation. Alexandria Engineering Journal,61(3),2204–2219. DOI 10.1016/j.aej.2021.07.005.Arfan, M., Alrabaiah, H., Rahman, M. U., Sun, Y. L., Hashim, A. S.et al. (2021). Investigation of fractal-fractional order model of COVID-19 in Pakistan under Atangana-Baleanu Caputo (ABC) derivative. Results in Physics,24,104046. DOI 10.1016/j.rinp.2021.104046.Boudaoui, A., El hadj Moussa, Y., Hammouch, Z., Ullah, S. (2021). A fractional-order model describing the dynamics of the novel coronavirus (COVID-19) with nonsingular kernel. Chaos, Solitons & Fractals,146,110859. DOI 10.1016/j.chaos.2021.110859.Zeb, A., Kumar, P., Erturk, V. S., Sitthiwirattham, T. (2022). A new study on two different vaccinated fractional-order COVID-19 models via numerical algorithms. Journal of King Saud University-Science,34(4),101914. DOI 10.1016/j.jksus.2022.101914.Alalyani, A., Saber, S. (2022). Stability analysis and numerical simulations of the fractional COVID-19 pandemic model. International Journal of Nonlinear Sciences and Numerical Simulation,2022,1–14. DOI 10.1515/ijnsns-2021-0042.Allegretti, S., Bulai, I. M., Marino, R., Menandro, M. A., Parisi, K. (2021). Vaccination effect conjoint to fraction of avoided contacts for a Sars-Cov-2 mathematical model. Mathematical Modelling and Numerical Simulation with Applications,1(2),56–66. DOI 10.53391/mmnsa.2021.01.006.Özköse, F., Yavuz, M., Senel, M. T., Habbireeh, R. (2022). Fractional order modelling of omicron SARS-CoV-2 variant containing heart attack effect using real data from the United Kingdom. Chaos, Solitons & Fractals,157,111954. DOI 10.1016/j.chaos.2022.111954.Alsakaji, H. J., Rihan, F. A., Hashish, A. (2022). Dynamics of a stochastic epidemic model with vaccination and multiple time-delays for COVID-19 in the UAE. Complexity,2022. DOI 10.1155/2022/4247800.Rihan, F. A., Alsakaji, H. J. (2021). Dynamics of a stochastic delay differential model for COVID-19 infection with asymptomatic infected and interacting people: Case study in the UAE. Results in Physics,28,104658. DOI 10.1016/j.rinp.2021.104658.Pandey, P., Gómez-Aguilar, J. F., Kaabar, M. K., Siri, Z., Abd Allah, A. M. (2022). Mathematical modeling of COVID-19 pandemic in India using Caputo-Fabrizio fractional derivative. Computers in Biology and Medicine,145,105518. DOI 10.1016/j.compbiomed.2022.105518.Mahrouf, M., Boukhouima, A., Zine, H., Lotfi, E. M., Torres, D. F.et al. (2021). Modeling and forecasting of COVID-19 spreading by delayed stochastic differential equations. Axioms,10(1),18. DOI 10.3390/axioms10010018.Kilbas, A. A., Marichev, O. I., Samko, S. G. (1993). Fractional integrals and derivatives (Theory and applications). Gordon and Breach, Switzerland: Elesvier.Podlubny, I., Magin, R. L., Trymorush, I. (2017). Niels Henrik Abel and the birth of fractional calculus. Fractional Calculus and Applied Analysis,20(5),1068–1075. DOI 10.1515/fca-2017-0057.Yusuf, A., Acay, B., Mustapha, U. T., Inc, M., Baleanu, D. (2021). Mathematical modeling of pine wilt disease with Caputo fractional operator. Chaos, Solitons & Fractals,143,110569. DOI 10.1016/j.chaos.2020.110569.Diekmann, O., Heesterbeek, J. A. P., Metz, J. A. (1990). On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. Journal of Mathematical Biology,28(4),365–382. DOI 10.1007/BF00178324.Van den Driessche, P., Watmough, J. (2002). Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Mathematical Biosciences,180(1–2),29–48. DOI 10.1016/S0025-5564(02)00108-6.Garrappa, R. (2010). On linear stability of predictor-corrector algorithms for fractional differential equations. International Journal of Computer Mathematics,87(10),2281–2290. DOI 10.1080/00207160802624331.Amouch, M., Karim, N. (2021). Modeling the dynamic of COVID-19 with different types of transmissions. Chaos, Solitons & Fractals,150,111188. DOI 10.1016/j.chaos.2021.111188.Zine, H., Boukhouima, A., Lotfi, E. M., Mahrouf, M., Torres, D. F.et al. (2020). A stochastic time-delayed model for the effectiveness of Moroccan COVID-19 deconfinement strategy. Mathematical Modelling of Natural Phenomena,15,50. DOI 10.1051/mmnp/2020040.Kada, D., Labzai, A., Balatif, O., Rachik, M., Labriji, E. H. (2020). Spread of COVID-19 in Morocco discrete mathematical modeling: Optimal control strategies and cost-effectiveness analysis. Journal of Mathematical and Computational Science,10(5),2070–2093.Musa, S. S., Qureshi, S., Zhao, S., Yusuf, A., Mustapha, U. T.et al. (2021). Mathematical modeling of COVID-19 epidemic with effect of awareness programs. Infectious Disease Modelling,6,448–460. DOI 10.1016/j.idm.2021.01.012.Abu-Shady, M., Kaabar, M. K. (2021). A generalized definition of the fractional derivative with applications. Mathematical Problems in Engineering,2021. DOI 10.1155/2021/9444803.