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

In this study, the spreading of the pandemic coronavirus disease (COVID-19) is formulated mathematically. .e 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. .erefore, 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. .e solution properties of the system are discussed. .e equilibria and the basic reproduction number are computed. .e local and global stabilities are studied. .e 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.


Introduction
e 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. is 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]. e 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  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. ey 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][12][13][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. e impact of curfew, media alert coverage, and social distance on the dynamical behavior of the model is considered. e paper is organized as follows. Section 2 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. e 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.

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. erefore, the population can be divided into five main compartments; these are the susceptible individuals at time t denoted by S(t) in which the individuals are at risk of infection with COVID-19. e exposed individuals at time t are denoted by E(t), in which the individuals have contracted the COVID-19 but are not yet infectious. e infected individuals at time t are represented as I(t) in which the individuals are capable of spreading the disease, and they are having mild symptoms or not. e hospitalized infected individuals at time t are denoted by I h (t), 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 R(t) in which the individuals have recovered from the disease and are no more infectious. erefore, the mathematical model that describes the dynamics of COVID-19 can be described as follows: All the parameters are nonnegative real numbers, and their descriptions are explained in Table 1. Note that b(I) � (b 1 − b 2 I/m + I) is the contact rate after the media alert (see [15]). e 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.
is 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 (b 2 I/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 (b 2 I/m + I) approaches its maximum that is given by b 2 . However, the value of (b 2 I/m + I) equals half of the maximum b 2 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 b 1 ≥ b 2 .
Recall that the interaction functions between the compartments given by the right-hand side of the system (1) are continuous and have continuous partial derivatives. erefore, they are Lipschitz functions, and hence system (1) has a unique solution. Proof. Let (S(t), E(t), I(t), I h (t), R(t)) be any solution of system (1) that initiates at the nonnegative value (S(0), E(0), I(0), I h (0), R(0)).
By taking the derivative of N(t) with respect to the time along with the solution of system (1), it is obtained that en, using direct computation gives that Accordingly, it is obtained that which implies that N(t) ≤ A/d 0 , 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: □

Existence of the Equilibrium Points and Basic Reproduction Number
Note that since the variable R, 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 I h 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 where (I, I h ) represents the solution values of system (7). Accordingly, the following system will be studied instead of system (1): Clearly, system (7) has two equilibrium points: the disease-free equilibrium point and the endemic equilibrium point. e disease-free equilibrium point is given by 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 where i, j � 1, 2, . . . , m, ρ 0 is the disease-free equilibrium point, while x i denotes the number of individuals in the i th infected compartment. However, F i (x i ) is the rate of e direct contact rate between S and I in the absence of media alerts b 2 e direct contact rate between S and I in the existence of media alerts m Nonresponse rate of the individuals to the media alerts p e rate of incubation termination and transition from uninfected to infectious individuals q e fraction of the exposed individuals move to the infected compartment, while the rest (1 − q) moves to the hospitalized infected population r 1 e recovery rate due to natural immune response r 2 e recovery rate due to treatment v e fraction of the nonrecovered individuals who start suffering severe symptoms move to the hospitalized infected population d 0 e natural death rate d 1 e death rate of disease in the infectious compartment d 2 e death rate due to disease in the hospitalized infected population International Journal of Differential Equations 3 gives the rate of change of x i . Note that F i 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 R 0 , is given by the spectral radius (dominant eigenvalue) of the matrix FV − 1 . It is well known that the basic reproduction number R 0 is one of the most crucial quantities in infectious diseases as R 0 measures how contagious a disease is. For R 0 < 1, the disease is expected to stop spreading, but for R 0 � 1, an infected individual can infect on an average 1 person; that is, the spread of the disease is stable. e disease can spread and become epidemic if R 0 > 1. Accordingly, regarding system (7), it is obtained that where . en, we get erefore, Consequently, the basic reproduction number of system (7) is determined as Now, the endemic equilibrium point of system (7) can be determined by where while E * is a positive root of the second-order polynomial equation: where . erefore, equation (11) has a unique positive root, given if and only if the following condition holds: Moreover, in order to obtain a positive endemic equilibrium point ρ 1 for system (7), we should have

Local Stability Analysis
is section treats the local stability of system (7) using the linearization technique. e Jacobian matrix for system (7) at the point (S, E, I, I h ) can be written as follows: where the symbols B i , i � 1, 2, 3, are given in equation (10). erefore, the local stability analysis for the disease-free equilibrium point can be studied in the following theorem.

Theorem 2.
e disease-free equilibrium point (ρ 0 ) of system (7) is locally asymptotically stable provided the following condition holds: Proof. According to equation (18), the Jacobian matrix at the point (ρ 0 ) is given by Direct computations show that this Jacobian matrix has the following characteristic equation: where erefore, 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 � − d 0 and λ 4 � − B 3 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. us, 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. e endemic equilibrium point (ρ 1 ) of system (7) is locally asymptotically stable provided that where b ij , 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 Hence, the characteristic equation of this Jacobian matrix can be written as where , accordingly, by substituting the values of b ij and then using condition (22) along with the existence condition (16), we obtain that A 1 , A 3 , and Δ are positive. erefore, all the roots (eigenvalues) of the third-order polynomial equation in (24) have negative real parts. Furthermore, since the fourth eigenvalue λ 4 � b 44 � − B 3 < 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.

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. erefore, the following theorem discusses the global stability of the free-disease equilibrium point.
Clearly, we have dV 0 /dt ≤ 0 if and only if the local stability condition (19) for the free-disease equilibrium point holds. Hence, V 0 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 * , I * h ) is locally asymptotically stable, then it is a global asymptotically stable provided that the following conditions hold: where all the symbols π ij are given in the proof.

The Local Bifurcation Analysis
Rewrite system (7) in the form as follows: where the functions f i , 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 where V � (v 1 , v 2 , v 3 , v 4 ) T be any nonzero vector. erefore, 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 ρ
Proof. Straightforward computation shows that the Jacobian matrix of system (7) at (ρ 0 , b 1 * ) can be written as erefore, the eigenvalues of J 0 can be written as International Journal of Differential Equations Accordingly, the disease-free equilibrium point becomes the nonhyperbolic point. Now, define U 1 � (x 1 , x 2 , x 3 , x 4 ) T as the eigenvector corresponding to the eigenvalue λ 03 � 0. Direct computation gives that where where η 4 � qp/B 1 > 0, with y 3 be any nonzero real number. Define U 2 � (y 1 , y 2 , y 3 , y 4 ) T as the eigenvector corresponding to the eigenvalue λ 03 � 0 of the matrix J T 0 . Direct computation gives that Hence, However, substituting the values of U 1 in equation (32) instead of V and then computing Clearly, U T 2 D 2 F(ρ 0 , b * 1 )(U 1 , U 1 ) ≠ 0. Hence, all the transcritical bifurcation conditions in the sense of Sotomayor's theorem are satisfied. erefore, the proof is complete.
It is observed that although system (7) becomes a nonhyperbolic near the endemic equilibrium point ρ 1 as where b ij , 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. is is can be verified directly by substituting the values of ρ 1 and p * in the value of A 3 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 of J(ρ 1 , p * ) and their transpose are, respectively, given by where r 1 )] > 0, with z 3 and w 3 are any nonzero real numbers. Moreover, it is obtain that Accordingly, any type of the local bifurcations in the sense of Sotomayor cannot occur.

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.
is 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. erefore, 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 Table 2. 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.
For the parameter values used in Table 2, system (1) is solved, and the trajectories are obtained as shown in Figure 1.
According to Figure 1, system (1) is settled at the endemic point when the basic reproduction number is R 0 � 2.4771, which means that the disease is out of control. However, for the data used in Table 1 with different values of 8 International Journal of Differential Equations parameter m, that m � 0, 100, 200, 300, respectively, the trajectories of system (1) are drawn in Figure 2. It is clear from Figure 2 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 Table 2 with c � 0.9, it is observed that the basic reproduction number decreases to R 0 � 0.6193 and the trajectory of system (1) is shown in Figure 3. It is clear from Figure 3 that when R 0 < 1, system (1) approaches asymptotically to disease-free equilibrium point and that confirms the theoretical results in previous sections.
On the other hand, increasing the parameter b 1 from b 1 � 0.001 to b 1 � 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.
Clearly, from Figure 4, the bifurcation occurs in system (1) and the trajectory of system (1) transfers from diseasefree point to endemic point as the bifurcation parameter b 1 increases, passing through the value b * 1 and then becomes b 1 � 0.01. It is noted that the basic reproduction number increases from the value R 0 � 0.6193 to the value R 0 � 6.192717 passing the bifurcation value R 0 � 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 R 0 � 2.4771 to the value R 0 � 0.9908, and the trajectory of system (1) is shown in Figure 5.
Again, Figure 5 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   (1), for the data given in Table 1 with different values of m, approaches to endemic point. e blue color for m � 0, green color for m � 100, red color for m � 200, and cyan color for m � 300. observed that quantitative effects exist rather than qualitative one so that the population size of the different compartments in the system changes slightly.

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. e 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 satisfies R 0 < 1 and unstable for R 0 > 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. e 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  Table 2 with c � 0.9 and b 1 � 0.01.  Table 2 with q � 0.2.

International Journal of Differential Equations
(1) as suggested by World Health Organization reports and adopted by many researchers, and the following observations are obtained. e 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 on R 0 < 1 and R 0 > 1, respectively. e model is solved numerically using the Runge-Kutta fourth-order method with the help of Matlab. e 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 b 2 , which increases the term (b 2 I/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 diseasefree equilibrium point.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest.