# The Impact of Media Coverage and Curfew on the Outbreak of Coronavirus Disease 2019 Model: Stability and Bifurcation

Authors: Afrah K. S. Al-Tameemi; Raid K. Naji
Journal: International Journal of Differential Equations
 (2021)

---

## Abstract

In this study, the spreading of the pandemic coronavirus disease (COVID-19) is formulated mathematically. The objective of this study is to stop or slow the spread of COVID-19. In fact, to stop the spread of COVID-19, the vaccine of the disease is needed. However, in the absence of the vaccine, people must have to obey curfew and social distancing and follow the media alert coverage rule. In order to maintain these alternative factors, we must obey the modeling rule. Therefore, the impact of curfew, media alert coverage, and social distance between the individuals on the outbreak of disease is considered. Five ordinary differential equations of the first-order are used to represent the model. The solution properties of the system are discussed. The equilibria and the basic reproduction number are computed. The local and global stabilities are studied. The occurrence of local bifurcation near the disease-free equilibrium point is investigated. Numerical simulation is carried out in applying the model to the sample of the Iraqi population through solving the model using the Runge–Kutta fourth-order method with the help of Matlab. It is observed that the complete application of the curfew and social distance makes the basic reproduction number less than one and hence prevents the outbreak of disease. However, increasing the media alert coverage does not prevent the outbreak of disease completely, instead of that it reduces the spread, which means the disease is under control, by reducing the basic reproduction number and making it an approachable one.


---

## Body

## 1. Introduction

The most considerable species in habitat are viruses. Since they cannot reproduce themselves, they are known as parasites. Some of these viruses, when they are replicated, cause dangerous diseases in the host, whether humans or animals. One of these species is the coronavirus (COVID-19). It is a disease caused by a new coronavirus called SARS-CoV-2. An outbreak of pneumonia caused by a novel coronavirus (COVID-19) began on 31 December 2019, in Wuhan, China. This pandemic that rapidly spread over 215 countries as of March 2020 continues to hit heavy public health and socio-economic burden in many parts of the world, including in Iraq. From January 22, 2020, up to March 3, 2021, it was reported that over 115 million confirmed cases and about 2,565,310 deaths around the world [1].The experience in China showed that the use of complete curfew and the speed diagnosis of the infected individuals have a strong impact on the dynamics of the epidemic [2]. It is well known that there is a long history of applying mathematical models in epidemiology to understand and control the outbreak of disease. Carcione et al. [3] implemented an SEIR model to compute the COVID-19 infected population and the number of casualties of this epidemic. Eikenberry et al. [4] developed a compartmental model for assessing the community-wide impact of masks used in the US states of New York and Washington by the general, asymptomatic public, a portion of which may be asymptomatically infectious. Maji [5] formulated a mathematical model for the COVID-19 outbreak by introducing a quarantine class with media-induced fear in the disease transmission rate to analyze the dynamic behavior of this epidemic. Asamoah et al. [6] developed a deterministic model to study the transmission dynamics of the disease with two categories of the susceptibles: Immigrant Susceptibles and Local Susceptible. Peter et al. [7], proposed a new mathematical model to investigate the recent outbreak of coronavirus disease (COVID-19) in Pakistan. Das et al. [8] formulated a mathematical model with comorbidity to study the transmission dynamics as well as an optimal control-based framework to diminish COVID-19. Musa et al. [9] employed an SEIR model to study the transmission dynamics of SARS-CoV-2 outbreaks in Nigeria; the model incorporates a different group of populations, namely, high and moderate risk populations. They used the model to investigate the influence of each population on the overall transmission dynamics. Recently, Paul et al. [10] proposed a simple population dynamics solution based on the incidence-fitness relationship in predicting that a plateau or steady-state of SARS-CoV-2 will be reached using the basic concept of geometry. Many other researchers have been proposed and studied mathematical models to describe the dynamic of COVID-19 (see [11–14] and the references therein).In view of the above, in this paper, however, a mathematical model that describes the spreading of COVID-19 is proposed and studied. It is assumed that the population is divided into five compartments including susceptible, exposed, infected, hospitalized infected, and recovered. The impact of curfew, media alert coverage, and social distance on the dynamical behavior of the model is considered. The paper is organized as follows. Section2 treats the formulation of the model with standard properties of the solution of the model. Section 3 interests in the existence of equilibria and computing the basic reproduction number. Section 4 studies local stability. While global stability is discussed in Section 5. The local bifurcation analysis is investigated in Section 6. However, Section 7 contained the numerical simulation for the model. Finally, Section 8 gives the discussion and conclusion of this study.


## 2. Formulation of the Model

In this section, the dynamics of COVID-19 are mathematically modeled using first-order ordinary differential equations. It is known that COVID-19 is transmitted mainly from humans to humans through direct contact between them. Accordingly, the population is divided into a number of separated compartments depending on the rules of the health section. Therefore, the population can be divided into five main compartments; these are the susceptible individuals at timet denoted by St in which the individuals are at risk of infection with COVID-19. The exposed individuals at time t are denoted by Et, in which the individuals have contracted the COVID-19 but are not yet infectious. The infected individuals at time t are represented as It in which the individuals are capable of spreading the disease, and they are having mild symptoms or not. The hospitalized infected individuals at time t are denoted by Iht, in which the individuals have severe symptoms or quarantined due to getting mild symptoms, and they put in the hospital to get treatment. Finally, the recovered individuals at time t are denoted by Rt in which the individuals have recovered from the disease and are no more infectious. Therefore, the mathematical model that describes the dynamics of COVID-19 can be described as follows:(1)dSdt=A−1−cb1−b2Im+IIS−d0S,dEdt=1−cb1−b2Im+IIS−qpE−1−qpE−d0E,dIdt=qpE−r1I−v1−r1I−d0+d1I,dIhdt=1−qpE+v1−r1I−r2Ih−d0+d2Ih,dRdt=r1I+r2Ih−d0R.All the parameters are nonnegative real numbers, and their descriptions are explained in Table1. Note that bI=b1−b2I/m+I is the contact rate after the media alert (see [15]). The choice of this function to describe the media alert is due to the following reasons. Most human beings will protect themselves from infection using all possible ways as soon as infected individuals are reported by media coverage. This reduces the transmission rate more or less. Accordingly, as the infected individuals increase, the susceptible individuals reduce their contact with other individuals to avoid the infection. Moreover, the term b2I/m+I refers to the reduced value of the contact (transmission) rate in the case of reporting of the infectious individuals. Clearly, as the infected individuals increase without bound, the term b2I/m+I approaches its maximum that is given by b2. However, the value of b2I/m+I equals half of the maximum b2 when the reported infective number arrives at m. Finally, since the media alert report cannot prevent the disease from spreading completely, it is assumed that b1≥b2.Table 1
The parameters description.
ParameterDescriptionAThe recruitment ratecThe probability of individuals responding to the curfewb1The direct contact rate betweenS and I in the absence of media alertsb2The direct contact rate betweenS and I in the existence of media alertsmNonresponse rate of the individuals to the media alertspThe rate of incubation termination and transition from uninfected to infectious individualsqThe fraction of the exposed individuals move to the infected compartment, while the rest1−q moves to the hospitalized infected populationr1The recovery rate due to natural immune responser2The recovery rate due to treatmentvThe fraction of the nonrecovered individuals who start suffering severe symptoms move to the hospitalized infected populationd0The natural death rated1The death rate of disease in the infectious compartmentd2The death rate due to disease in the hospitalized infected populationRecall that the interaction functions between the compartments given by the right-hand side of the system (1) are continuous and have continuous partial derivatives. Therefore, they are Lipschitz functions, and hence system (1) has a unique solution.Theorem 1.
All solutions of system (1) which initiate in R+5 are uniformly bounded.Proof.
LetSt,Et,It,Iht,Rt be any solution of system (1) that initiates at the nonnegative value S0,E0,I0,Ih0,R0.
Consider the functionNt=St+Et+It+Iht+Rt.
By taking the derivative ofNt with respect to the time along with the solution of system (1), it is obtained that(2)dNdt=dSdt+dEdt+dIdt+dIhdt+dRdt.
Then, using direct computation gives that(3)dNdt=A−d0S−d0E−d0+d1I−d0+d2Ih–d0R.
Accordingly, it is obtained that(4)dNdt≤A−d0N,which implies that Nt≤A/d0, as t⟶∞. Hence, all the solutions of system (1) are uniformly bounded, and therefore we finished the proof.
According to system (1), the feasible region of it can be written as follows:(5)Ω=S,E,I,Ih,R∈R+5:0≤S+E+I+Ih+R≤Ad0.


## 3. Existence of the Equilibrium Points and Basic Reproduction Number

Note that since the variableR, which represents the recovered individuals, dose not appear in the first four equations of system (1), hence one can solve the following system instead of system (1) and then substitute the solution values of I and Ih in the fifth equation of system (1) to solve it separately as a linear differential equation with respect to the variable R. Straightforward computation gives that the solution of the fifth equation as t⟶∞ can be written as(6)Rt=r1I^+r2I^hd0,where I^,I^h represents the solution values of system (7).Accordingly, the following system will be studied instead of system (1):(7)dSdt=A−1−cb1−b2Im+IIS−d0S,dEdt=1−cb1−b2Im+IIS−qpE−1−qpE−d0E,dIdt=qpE−r1I−v1−r1I−d0+d1I,dIhdt=1−qpE+v1−r1I−r2Ih−d0+d2Ih.Clearly, system (7) has two equilibrium points: the disease-free equilibrium point and the endemic equilibrium point. The disease-free equilibrium point is given by ρ0=S0,0,0,0, where(8)S0=Ad0.Since the basic reproduction number is computed at this equilibrium point, hence in the following the computation of this number is carried out, and then the endemic equilibrium point is determined.Now, in order to determine the basic reproduction number, the “next-generation method” or “Spectral Radius method” is used (see [16, 17]). Consider an epidemic model having n different compartments from which m compartments contained infected individuals with the disease, then the next-generation matrix (operator) is given by FV−1, where(9)F=∂Fiρ0∂xjandV=∂Viρ0∂xj,where i,j=1,2,…,m, ρ0 is the disease-free equilibrium point, while xi denotes the number of individuals in the ith infected compartment. However, Fixi is the rate of appearance of new infections in the compartment i, while Vixi=Vi−xi−Vi+xi where Vi+xi represents the rate of the shift of individuals into ith compartment by all other means and Vi−xi denotes the rate of the shift of individuals out of the ith compartment. The difference Fixi−Vixi gives the rate of change of xi. Note that Fi should include only infections that are newly arising but does not include terms that describe the shift of infectious individuals from one infected compartment to another. Finally, the basic reproduction number, that is denoted by R0, is given by the spectral radius (dominant eigenvalue) of the matrix FV−1. It is well known that the basic reproduction number R0 is one of the most crucial quantities in infectious diseases as R0 measures how contagious a disease is. For R0<1, the disease is expected to stop spreading, but for R0=1, an infected individual can infect on an average 1 person; that is, the spread of the disease is stable. The disease can spread and become epidemic if R0>1.Accordingly, regarding system (7), it is obtained that(10)F=01−cb1S00000000,V=B100−qpB20−1−qp−v1−r1B3,where B1=p+d0, B2=r1+v1−r1+d0+d1, and B3=r2+d0+d2. Then, we get(11)V−1=1B1B2B3B2B300qpB3B1B30qpv1−r1+1−qpB2v1−r1B1B1B2.Therefore,(12)FV−1=1B1B2B3qp1−cb1S0B31−cb1S0B1B30000000.Consequently, the basic reproduction number of system (7) is determined as(13)R0=qp1−cb1S0B3B1B2B3.Now, the endemic equilibrium point of system (7) can be determined by ρ1=S∗,E∗,I∗,Ih∗, where(14)S∗=A−B1E∗d0,I∗=qpE∗B2,Ih∗=pB2–qpB4E∗B2B3,where B4=r1+d0+d1, while E∗ is a positive root of the second-order polynomial equation:(15)α1E2+α2E+α3=0,where α1=−1−cb1–b2B1qp2<0, α2=qp1−cb1–b2Aqp−1−cb1mB1B2−B1B2, and α3=mB21−cb1qpA−B1B2d0=B1B2d0mB2R0−1.Therefore, equation (11) has a unique positive root, given by E∗=−α2/2α1−1/2α1α22−4α1α3, if and only if the following condition holds:(16)R0>1.Moreover, in order to obtain a positive endemic equilibrium pointρ1 for system (7), we should have(17)E∗<AB1.


## 4. Local Stability Analysis

This section treats the local stability of system (7) using the linearization technique. The Jacobian matrix for system (7) at the point S,E,I,Ih can be written as follows:(18)J=−1−cb1–b2Im+II−d00−1−cb1–b22m+IIm+I2S01−cb1–b2Im+II−B11−cb1–b22m+IIm+I2S00qp−B2001−qpv1−r1−B3,where the symbols Bi,i=1,2,3, are given in equation (10). Therefore, the local stability analysis for the disease-free equilibrium point can be studied in the following theorem.Theorem 2.
The disease-free equilibrium pointρ0 of system (7) is locally asymptotically stable provided the following condition holds:(19)R0<1.Proof.
According to equation (18), the Jacobian matrix at the point ρ0 is given by(20)Jρ0=−d00−1−cb1S000−B11−cb1S000qp−B2001−qpv1−r1−B3.
Direct computations show that this Jacobian matrix has the following characteristic equation:(21)−d0–λ−B3–λλ2−Tλ+D=0,where T=−B1+B2<0 and D=B1B2−1−cqpb1S0=B1B21−R0. Therefore, all the eigenvalues of the matrix Jρ0 are real and negative if and only if condition (19) holds. Moreover, since the other two eigenvalues λ1=−d0 and λ4=−B3 are negative, hence the disease-free equilibrium point ρ0 of system (7) is locally asymptotically stable under condition (19). However, if condition (19) is revised, then ρ0 is an unstable point. Thus, this completes the proof.
According to theorem (6), the basic reproduction number represents the vital value for controlling the outbreak of COVID-19, so that the disease will disappear provided that the basic reproduction number is less than one while the COVID-19 will outbreak when the basic reproduction number is greater than one.Theorem 3.
The endemic equilibrium pointρ1 of system (7) is locally asymptotically stable provided that(22)b22b33>b23b32,where bij,i,j=1,2,3, are given in the proof.Proof.
Substituting the endemic equilibrium pointρ1 in the Jacobian matrix given by equation (18) gives that(23)Jρ1=bij4×4,where b11=−1−cb1–b2I∗/m+I∗I∗−d0, b13=−1−cb1–b22m+I∗I∗/m+I∗2S∗, b21=1−cb1–b2I∗/m+I∗I∗, b22=−B1, b23=1−cb1–b22m+I∗I∗/m+I∗2S∗, b32=qp, b33=−B2, b42=1−qp, b43=v1−r1, and b44=−B3.
Hence, the characteristic equation of this Jacobian matrix can be written as(24)b44–λλ3+A1λ2+A2λ+A3=0,where A1=−b11+b22+b33, A2=b11b22+b11b33+b22b33–b23b32, and A3=−b11b22b33+b13b21b32−b11b23b32.
WhileΔ=A1A2−A3=A1b11b22+b22b33−b23b32−b11b33b11+b33, accordingly, by substituting the values of bij and then using condition (22) along with the existence condition (16), we obtain that A1,A3, and Δ are positive. Therefore, all the roots (eigenvalues) of the third-order polynomial equation in (24) have negative real parts. Furthermore, since the fourth eigenvalue λ4=b44=−B3<0, then the proof is complete.
From theorem (7), it concluded that the COVID-19 will be endemic and the disease is out of control if in addition to the basic reproduction number greater than one, condition (22) is satisfied.


## 5. Global Stability

In this section, the stability analysis of system (7) is studied using the Lyapunov function in order to obtain the global stability conditions or the basin of attraction for each equilibrium point. Therefore, the following theorem discusses the global stability of the free-disease equilibrium point.Theorem 4.
Assume that free-disease equilibrium pointρ0=A/d0,0,0,0 is locally asymptotically stable, then it is global asymptotically stable.Proof.
Consider the following positive definite real valued function:(25)V0S,E,I,Ih=S−S0−S0lnSS0+E+I+Ih.
Clearly,V0:R+4⟶R, V0ρ0=0, and V0x>0 where x≠ρ0. Moreover, direct computation gives that(26)dV0dt=∂V0∂sdsdt+dEdt+dIdt+dIhdt=−1−cb1−b2Im+IIS+1−cb1−b2Im+IIS0−d0SS−S02+1−cb1–b2Im+IIS−d0E−r1+d0+d1I−r2+d0+d2Ih≤−d0SS−S02−d0E−r1+d0+d1−1−cb1S0I−r2+d0+d2Ih.
Clearly, we havedV0/dt≤0 if and only if the local stability condition (19) for the free-disease equilibrium point holds. Hence, V0 is negative definite for any initial point that belongs to R+4, and then the proof is complete.Theorem 5.
Assume that endemic equilibrium pointρ1=S∗,E∗,I∗,Ih∗ is locally asymptotically stable, then it is a global asymptotically stable provided that the following conditions hold:(27a)π122<23π11B1,(27b)π132<23π11B2,(27c)π232<49B1B2,(27d)π242<23B1B3,(27e)π342<23B2B3,where all the symbols πij are given in the proof.Proof.
Consider the positive definite function(28)V1=S−S∗22+E−E∗22+I−I∗22+Ih−Ih∗22.
Clearly,V1:R+4⟶R, V1ρ1=0, and V1x>0 where x≠ρ1. Moreover, direct computation gives that(29)dV1dt=∂V1∂sdsdt+∂V1∂EdEdt+∂V1∂IdIdt+∂V1∂IhdIhdt,=S−S∗−1−cb1−b2Im+IIS−d0S−S∗+1−cb1−b2I∗m+I∗I∗S∗,+E−E∗1−cb1–b2Im+IIS−B1E−E∗−1−cb1–b2I∗m+I∗I∗S∗,+I−I∗qpE−E∗−B2I−I∗,+Ih−Ih∗1−qpE−E∗+v1−r1I−I∗−B3Ih−Ih∗,=−π11S−S∗2−π12S−S∗E−E∗−B1E−E∗2,−π13S−S∗I−I∗−B2I−I∗2,−π23E−E∗I−I∗−π24E−E∗Ih−Ih∗,+π34I−I∗Ih−Ih∗−B3Ih−Ih∗2,where π11=d0+1−cI∗/m+Im+I∗mb11+I+b1−b2m+II∗, π12=1−cI∗/m+Im+I∗mb11+I+b1−b2m+II∗, π13=1−cS/m+Im+I∗mb1+b1−b2I+b1−b2m+II∗, π23=1−cS/m+Im+I∗mb1+b1−b2I+b1−b2m+II∗+qp, and π24=1−qp, π34=v1−r1.
Accordingly, by using conditions (27a)–(27e), we obtain that(30)dV1dt<−12π11S−S∗+23B1E−E∗2−12π11S−S∗+23B2I−I∗2,−23B1E−E∗+23B2I−I∗2−1223B1E−E∗+B3Ih−Ih∗2,−1223B2I−I∗−B3Ih−Ih∗2.
Clearly, the derivativedV1/dt is negative definite for any initial point belonging to R+4, and then the proof is complete.


## 6. The Local Bifurcation Analysis

Rewrite system (7) in the form as follows:(31)dXdt=FX,X=S,E,I,IhT,F=f1,f2,f3,f4T,where the functions fi,i=1,2,3,4, are the interaction functions given on the right-hand side of system (7). Accordingly, the second derivative of F with respect to X can be written as(32)D2FXV,V=−21−cb1–b22m+IIm+I2v1v3+1−c2m2b2m+I3Sv3221−cb1–b22m+IIm+I2v1v3−1−c2m2b2m+I3Sv3200,where V=v1,v2,v3,v4T be any nonzero vector. Therefore, the local bifurcation condition at the disease-free equilibrium point is established by the following theorem.Theorem 6.
System (7) exhibits a transcritical bifurcation near the disease-free equilibrium point ρ0, where R0=1 or b1=B1B2/qp1−cS0≡b1∗.Proof.
Straightforward computation shows that the Jacobian matrix of system (7) at ρ0,b1∗ can be written as(33a)J0=Jρ0,b1∗=−d00−1−cb1∗S000−B11−cb1∗S000qp−B2001−qpv1−r1−B3.
Therefore, the eigenvalues ofJ0 can be written as(33b)λ01=−d0,λ02=−B1+B2,λ03=0,λ04=−B3.
Accordingly, the disease-free equilibrium point becomes the nonhyperbolic point. Now, defineU1=x1,x2,x3,x4T as the eigenvector corresponding to the eigenvalue λ03=0. Direct computation gives that(34)U1=η1x3,η2x3,x3,η3x3T,where η1=−B1B2/qpd0<0, η2=B2/qp>0, and η3=1/B31−qB2/q+v1−r1>0, with x3 be any nonzero real number.(35)U2=0,η4y3,y3,0T,where η4=qp/B1>0, with y3 be any nonzero real number.
DefineU2=y1,y2,y3,y4T as the eigenvector corresponding to the eigenvalue λ03=0 of the matrix J0T. Direct computation gives that
Hence,(36)∂F∂b1=Fb1=−1−cSI1−cSI00⇛Fb1ρ0,b1∗=0000.
Then, it is obtained thatU2TFb1ρ0,b1∗=0. Furthermore, it is observed that(37)DFb1ρ0,b1∗=00−1−cS00001−cS0000000000.
Therefore,U2TDFb1ρ0,b1∗U1=1−cS0η4x3y3≠0. However, substituting the values of U1 in equation (32) instead of V and then computing D2Fρ0,b1∗ give that(38)U2TD2Fρ0,b1∗U1,U1=21−cη4−B1B22qp2d01−cS0−b2mS0x32y3.
Clearly,U2TD2Fρ0,b1∗U1,U1≠0. Hence, all the transcritical bifurcation conditions in the sense of Sotomayor’s theorem are satisfied. Therefore, the proof is complete.
It is observed that although system (7) becomes a nonhyperbolic near the endemic equilibrium point ρ1 as p=−b11b22b33/qb13b21−b11b23≡p∗ where bij,ij=1,2,3,4 are the elements if Jρ1 that is given in equation (23), it does not exhibit any types of local bifurcation in the sense of Sotomayor. This is can be verified directly by substituting the values of ρ1 and p∗ in the value of A3 that is given in equation (24) and then computing the eigenvectors corresponding to the resulting zero eigenvalues of Jρ1,p∗ and their transpose, respectively.
It is obtained that the eigenvectors corresponding to the zero eigenvalue ofJρ1,p∗ and their transpose are, respectively, given by(39)H1=μ1z3,μ2z3,z3,μ3z3T,H2=0,0,w3,0T,where μ1=−ξ2/ξ1+d0<0, μ2=ξ2d0/ξ1+d0p∗+d0>0, and μ3=1/B31−qp∗ξ2d0/ξ1+d0p∗+d0+v1−r1>0, with z3 and w3 are any nonzero real numbers. Moreover, it is obtain that(40)H2TFpρ1,p∗=qE∗w3≠0,H2TD2Fρ1,p∗H1,H1=0.
Accordingly, any type of the local bifurcations in the sense of Sotomayor cannot occur.


## 7. Numerical Simulation

To complete the analysis of the proposed system, it should be simulated numerically with real and reasonable parameter values and then the obtained results are validated with those obtained in the previous sections. It is well known that the model parameters affect the dynamical behavior of the model greatly. This occurs due to either error in their measurement or missing the parameter values. Accordingly, the impact of varying the model’s parameters on the behavior of the solution of the model is important. Therefore, the parameters that are having vital roles in the dynamics of the model will be given clear attention to understand their effects on the control of the spread of disease. Most of the parameter value ranges used in this section are taken from the World Health Organization reports and adopted by many references as stated in Table2. Since the population of Iraq is nonhomogenous distributed, therefore, we will confine ourselves to apply the proposed model on a sample that exists in a specific region to get more accurate results.Table 2
The baseline values of the parameters of system (1).
ParameterRangeValue usedSourceA30–1000 day−150[6]c[0,1]0.5Assumedb16.038 × 10−8–0.8158 day −10.001[6, 11, 12]b2(0,b1) day−10.0001[15]m[0,Ad0]10[15]p[1/14–1/3] day−10.1[4, 13, 14]q(0,1)0.5Assumedr1[1/14–1/3] day−10.08[4, 13, 14]r20.05–0.9 day−10.1[6]v0.01–0.5 day−10.3Assumedd00.000042–0.011 day−10.01[6, 11, 12]d10.0044–0.1850 day−10.001[6, 11, 12]d20.0044–0.1850 day−10.01[6, 11, 12]For the parameter values used in Table2, system (1) is solved, and the trajectories are obtained as shown in Figure 1.Figure 1
The trajectory of system (1) approaches to endemic equilibrium point.According to Figure1, system (1) is settled at the endemic point when the basic reproduction number is R0=2.4771, which means that the disease is out of control. However, for the data used in Table 1 with different values of parameter m, that m=0,100,200,300, respectively, the trajectories of system (1) are drawn in Figure 2.Figure 2
The trajectory of system (1), for the data given in Table 1 with different values of m, approaches to endemic point. The blue color for m=0, green color for m=100, red color for m=200, and cyan color for m=300.It is clear from Figure2 that with the increase in the value of m, the trajectory still settled at endemic but the populations of exposed, infected, and hospitalized infected increased too.Now, for the data given in Table2 with c=0.9, it is observed that the basic reproduction number decreases to R0=0.6193 and the trajectory of system (1) is shown in Figure 3. It is clear from Figure 3 that when R0<1, system (1) approaches asymptotically to disease-free equilibrium point and that confirms the theoretical results in previous sections.Figure 3
The trajectory of system (1) approaches to disease-free equilibrium point.On the other hand, increasing the parameterb1 from b1=0.001 to b1=0.01 with data used in Figure 3 makes the trajectory of system (7) return back to settle at the endemic equilibrium point as shown in Figure 4.Figure 4
The trajectory of system (1) approaches to endemic equilibrium point for the data in Table 2 with c=0.9 and b1=0.01.Clearly, from Figure4, the bifurcation occurs in system (1) and the trajectory of system (1) transfers from disease-free point to endemic point as the bifurcation parameter b1 increases, passing through the value b1∗ and then becomes b1=0.01. It is noted that the basic reproduction number increases from the value R0=0.6193 to the value R0=6.192717 passing the bifurcation value R0=1. Now, for the data given in Table 2 with reducing the value of q from the value q=0.5 to the value q=0.2, the basic reproduction number decreases from the value R0=2.4771 to the value R0=0.9908, and the trajectory of system (1) is shown in Figure 5.Figure 5
The trajectory of system (1) approaches to disease-free equilibrium point for the data in Table 2 with q=0.2.Again, Figure5 ensures the obtained results regarding the stability analysis of the disease-free equilibrium point. Finally, other parameters are also varying to study their effects on the dynamical behavior of system (1), and it is observed that quantitative effects exist rather than qualitative one so that the population size of the different compartments in the system changes slightly.


## 8. Discussion and Conclusion

In this paper, an epidemic model to describe the spread of COVID-19 is proposed using the first-order ordinary differential equations. The model takes into account the impact of the curfew and media coverage on the control of the spreading of disease. It is observed that model (1) has two equilibrium points: the disease-free point that is stable, which means the disease is under control from a biological point of view, when the basic reproduction number satisfiesR0<1 and unstable for R0>1, and the endemic equilibrium point that is stable, which means the disease out of control from a biological point of view, whenever the disease-free becomes unstable. All the properties of the solution of the model are studied. The local bifurcation near the disease-free point is investigated too. An application of the proposed model on the sample of Iraqi population of size 5000 has been done using baseline values of the parameters of system (1) as suggested by World Health Organization reports and adopted by many researchers, and the following observations are obtained.The model has only one type of attractor, approaches either to disease-free equilibrium point or to the endemic equilibrium point, which means either stop the outbreak of COVID-19 or spread the COVID-19, depending onR0<1 and R0>1, respectively. The model is solved numerically using the Runge–Kutta fourth-order method with the help of Matlab. The model is applied to a sample of the Iraqi population using the parameter value ranges as suggested by the World Health Organization reports and adopted by many references as stated in Table 2. It is observed that increasing the probability of an individual’s response to the curfew leads to reducing the basic reproduction number below one and then controlling the spread of COVID-19, and hence system (1) approaches asymptotically to the disease-free point. However, increasing the direct contact rate between S and I in the absence of media alerts or decreasing the social distance between the population individuals leads to an increase in the basic reproduction number above one, and then an increase in the outbreak of the disease occurred. Hence, the trajectory of system (1) approaches the endemic point. Although increasing the response rate of the individuals to the media alerts (that means decreasing m) cannot prevent the disease from spreading completely, it leads to reduce it up to a minimum rate. A similar result is observed by increasing the parameter b2, which increases the term b2I/m+I that reduces the value of the transmission rate when infectious individuals appear and are reported. Finally, decreasing the fraction of the exposed individuals moving to the infected compartment, which leads to increasing the number of individuals moving to hospitalized infected populations, can prevent the disease from spreading and the system approaches a disease-free equilibrium point.



---

*Source: 1892827-2021-11-24.xml*