Effect of Awareness Programs on the Epidemic Outbreaks with Time Delay

and Applied Analysis 3 to unaware class. The constants ], d represent the recovery rate and the natural death rate, respectively. All the above constants are assumed to be positive. Using the fact that X(t) + X m (t) + Y(t) = 1, the system (2) becomes as follows:


Introduction
Infectious diseases that cause mortality, disability, and social and economic disruption are a major threat to mankind, which are responsible for a quarter of all deaths annually in the world [1][2][3].Once an infectious disease appears and spreads in a region, the Center of Disease Control and Prevention will do its best to stop the propagation of the disease [4].One of the measures is to tell people the appropriate preventive knowledge of diseases as soon as possible through media and education which make people take precautions to reduce their chances of being infected [5].As the awareness disseminates, people will change their behaviors to alter their susceptibility.There is much evidence that media coverage can play an important role in the spread and control of infectious diseases [6][7][8][9][10][11].In [8] Cui and others established a framework of transmission coefficient with media coverage which is a decreasing function of the number of the infective individuals and they observed a classic threshold-type behavior, with the disease becoming extinct when  0 < 1 and going to a globally asymptotically stable equilibrium when  0 > 1, and they concluded that media coverage is critical in disease eradication.
The models studying the spread of infectious diseases are very useful in evaluating strategies to control the diseases in populations.Recently some authors studied the impacts of media coverage and education on the spread of infectious diseases in a given region [12][13][14].In [12] Misra et al. proposed a nonlinear mathematical model for the effects of awareness programs on the spread of infectious diseases and assumed the growth rate of awareness programs impacting the population is proportional to the number of infective individuals.The model analysis showed that the spread of an infectious disease can be controlled by using awareness programs but the disease remains endemic due to immigration.Yorke and London [15] proposed an SIS type compartmental model for sexually transmitted infections with the assumption that the whole population is aware of risk but only a certain proportion choose to respond by limiting their contacts with the infective population.As a result the spread of infection is controlled, leading to a reduction of the number of individuals becoming infected.A nonlinear mathematical model with delay to capture the dynamics of effect of awareness programs on the prevalence of any epidemic is proposed and analyzed [16], which assumed it increases at a rate proportional to the number of the deaths of infective individuals.Some scholars focus on the contact rate and most of them assume that media campaigning will aid in modifying the contact rate between susceptible and infective individuals [17][18][19][20][21][22][23][24][25][26][27].To prevent the unboundedness of the contact rate, Capasso [17] used a saturated incidence rate of the form /(1 + ); Liu and coworkers [23,24] used a nonlinear incidence rate given by   /(1 +  ℎ ) to incorporate the effect of behavioral changes.In [25], the authors study how media coverage influences the dynamics of infectious disease by using SIRS model with the contact rate ( 1 −  2 /( + )), where  1 ,  2 , and  are positive constants.Tchuenche and Bauch used an exponentially decreasing function  − ()  to reveal the force of media [27].Cui et al. used a similar function as [13] and developed an  model using incidence rate  (−)  to investigate the impact of media coverage on the transmission [28].Stability analysis of the model has shown that Hopf bifurcation can occur which causes oscillatory phenomena when  is sufficiently small.Numerical simulations suggested that the media impact was stronger when the basic reproduction number  0 > 1. Liu et al. [29] described the impact of media coverage using the transmission coefficient  (− 1 − 2 − 3 ) , where  is the number of hospitalized individuals.And this impact leads to the change of avoidance and contact patterns at both individual and community levels.Liu and Xiao [30] introduce a segmented function to describe the saturated media impact  −  when formulating an epidemic model.A Filippov epidemic model with media coverage is proposed to describe the real characteristics of media impact in the spread of infectious diseases by incorporating a piecewise continuous transmission rate  (−)  in [31].Mathematical and bifurcation analysis with regard to the local, global stability of equilibria and local sliding bifurcations are performed.Bhunu et al. [32] and Tchuenche and Bauch [27] focused on the different types of population in their work.
In order to better describe population mixed conditions, some authors study infectious diseases models on the different networks [33,34].Funk et al. have overlaid the model of information spread of a contagious disease on two, not necessarily identical networks, with more informed individuals acting to reduce their susceptibility [35].Liu et al. took into the random perturbation [36].In [37], the authors extended the classical SIRS epidemic model from a deterministic framework to a stochastic differential equation, and then they gave the conditions of existence of unique positive solution and the stochastic extinction and discussed the exponential p-stability and global stability.
Most of the articles, such as [12], assume that, due to awareness programs, driven by media, some susceptible individuals will avoid their contacts with the infectious individuals resulting in formation of a new class in the population and this newly formed aware class may contract infection only if they lose awareness.But we regard it is unreasonable.In fact sometimes even if persons are conscious of diseases, they will also be infected.Therefore we propose a mathematical model for predicting the future course of any epidemic by considering this newly formed aware class into the model.
The rest of this paper is organized as follows.In the next section, to capture the dynamics of effect of cumulative density of awareness programs on the prevalence of any epidemic a mathematical model is proposed and analyzed.Then we analyze the local and global stability of the disease-free and the unique endemic equilibrium in Section 3. Furthermore, in Section 4 we perform some numerical examples to validate the analytical findings in Section 3 and introduce the importance of the dissemination rate and implementation rate in disease control.In Section 5 we discuss the above content.

Mathematical Model and Equilibrium Analysis
In this paper due to the awareness programs about the disease, susceptible individuals are rarely in contact with the infective ones and form a different class, namely, aware susceptible class; thus the total population is divided into three classes, the susceptible population, the aware population, and the infective population, and the proportions of them in the total population are (), (), and   ().Assuming that at time  the cumulative density of awareness programs driven by media in that region is (), which increases at a rate proportional  to the infective population and consumes due to causes like inefficiency and psychological barriers at  0 , thus In fact people cannot take measures to protect themselves in time after the media reports the disease, so we introduce a time delay  that represents the interval between the report time and the time of taking measures.We assume that a proportion of infected individuals recover through treatment and, after recovery, a fraction  of recovered people will become aware and join the aware susceptible class whereas the remaining fraction  ( +  = 1) will become unaware susceptible.Keeping the above facts in mind, the dynamics of a model are governed by the following systems of nonlinear ordinary differential equations: Here (0) > 0, (0) ≥ 0, and   (0) ≥ 0. Assume diseases spread due to the direct contact between susceptible and infective individuals only.In the above model, the rate of immigration of susceptible population is . is the contact rate of aware susceptible with infective population and  ( < ) is the contact rate of unaware ones.The constant  represents the dissemination rate of awareness among unaware susceptible due to which they form a different class; then  0 denotes the rate of transfer of aware susceptible to unaware class.The constants ],  represent the recovery rate and the natural death rate, respectively.All the above constants are assumed to be positive.Using the fact that () +   () + () = 1, the system (2) becomes as follows: ( Now it is sufficient to study system (3) in detail rather than system (2).
For the analysis of system (3), we need the region of attraction which is given by the set: Ω = {(  , , ) ∈ R 3 + : 0 ≤   ,  ≤ 1, 0 ≤  < / 0 }, which attracts all solutions initiating in the interior of the positive orthant.
Using ( 5) and ( 6), we get Further, using ( 7) and ( 8) in ( 4), we obtain a quadratic equation in  * as where + ] 0 ( − ) +  0  ( + ) , Solving ( 9) we get We obtain  1 > 0 and  3 < 0 when  0 > 1 and get Remark.From the expression of  * , it is easy to note that ( * /) < 0 and ( * /) < 0, which shows that the equilibrium number of infective individuals decreases as the rate of dissemination and the implementation rate of awareness programs increase.

Stability Analysis
In this section we present the local and global stability of  0 and  * .
Proof.To establish the global stability of the disease-free equilibrium  0 , we use Lyapunov's method and consider the following positive definite function without : Now differentiating  with respect to , we get When  0 < 1, / ≤ 0. The largest compact invariant set in {(  , , ) ∈ Ω :   = 0} when  0 < 1 is the singleton { 0 }.Then LaSalle's invariance principle implies that  0 is globally asymptotically stable in Ω.

The Stability of the
where We have the disease-free equilibrium   0 (0, 0, 0) and the endemic equilibrium   * ( * ,  * ,  * ), where the stability of  * about system (3) is corresponding with   0 .The characteristic equation of the above system at   0 is where  is the eigenvalue and The stability of the endemic equilibrium  * of system (3) is stated in the following theorems.
Proof.When  = 0, the characteristic equation is of the form It is easy to see that ( From  3 −  4 > 0, we can get  * + ( − ) *  > 0; thus  1  2 − ( 3 −  4 ) > 0. According to Hurwitz criterion, we can know all the 's have negative real parts; then   0 is locally asymptotically stable.
When  > 0, notice that (18) does not have nonnegative real roots.If it has roots with nonnegative real parts they must be complex and should have been obtained from a pair of complex conjugate roots which cross the imaginary axis.Consequently, (18) must have a pair of purely imaginary solutions for some  > 0. Assume that  =  ( > 0) is a root of (18) without loss of generality.That is the case if and only if  satisfies the equation Separating the real and imaginary parts, we have the following system, satisfied by : To eliminate the trigonometric functions we square both sides of each equation above and we add the squared above equations to obtain the following forth order equation in : To reduce this fourth order equation in  to a quadratic equation let  =  2 and denote the coefficients by where So all the coefficients of ( 25) are positive numbers.Then according to Lemma 3.3.1 in [38], (25) has no positive real roots; that is, we may not get any positive value of , which satisfy the transcendental equation (18).So all the 's have negative real parts for all values of the delay  ≥ 0; then   0 is locally asymptotically stable.Thus when  ≥ 0,  * is locally asymptotically stable if  0 > 1.
Proof.Using Lyapunov's method, we consider the following positive function: The derivative of  along the system is given by The largest compact invariant set when   = 0 is the singleton {  0 }.Then LaSalle's invariance principle implies that   0 is globally asymptotically stable; that is,  * is globally asymptotically stable in Ω.

Numerical Simulations and Results
To check the feasibility of our analysis of  > 0, we present some numerical computations in this section using Matlab by choosing the following set of parameter values:  = 0.35,  = 0.08,  0 = 0.02,  = 0.2,  = 0.002, ] = 0.43,  = 0.15,  = 0.85,  = 0.002,  0 = 0.02, and  = 1 when  0 < 1.Let  = 0.5; it may be checked that the condition  0 > 1 The basic reproduction number  0 , for the above set of parameter values, is found to be 1.157.
For the above parameter values, we select five sets of different initial starts; then the computer generated graphs of aware population, infective population, and cumulative density of awareness programs, respectively, have been drawn in Figure 1, which shows that all the trajectories initiating inside the region of attraction approach towards  0 and  * , respectively.Both of the equilibria  0 and  * are locally asymptotically stable for given set of parameter values that numerical simulations support the analysis given in Section 3. In fact, they are globally asymptotically stable in Ω as we have proved.
In the following, we research the relationships of  *  ,  * and the dissemination rate , the implementation rate  separately.We make ,  change from 0 to 0.4 and get the performances of   () and (), the trajectories of which with respect to time  for different  and  are shown in Figures 2 and 3, respectively.And there are no awareness programs when  and  are equal to zero.As shown in Figures 2 and   3,  *  both increase and  * both reduce as the increase of  and , which proves the conclusions of the remark.And ,  are greater influence on  * than  *  , which state awareness programs have a positive effect on prevention of diseases.In addition the reason why  *  ( * ) has a similar trend as the variations of  and  is that  and  have a similar influence on  *  ( * ).From the figure we obtain that  can postpone the time of the balance of equilibrium; thus we can have more time to formulate measures to prevent diseases.There really is an effort here to make it clear that  and  (awareness programs) play a key role in the prevention and control of diseases.

Discussion
The media is widely acknowledged as a key tool for influencing people's behaviors towards the disease to devise proper policies for controlling the epidemic.Awareness programs through media make people be aware about the disease and take various precautions to reduce their chances of being infected.In this paper, we propose and analyze a mathematical model to study the effect of awareness programs driven by media and the delay on the prevalence of an infectious disease.It is assumed that pathogens are transmitted via direct contact between the susceptible and the infective populations.Assume further that cumulative density of awareness programs increases at a rate proportional to the infective population.The model exhibits two equilibria; the disease-free equilibrium has been shown to be stable for basic reproduction number  0 < 1.For  0 > 1, it becomes unstable, which leads to the existence of an endemic equilibrium.The endemic equilibriums are globally asymptotically stable.The delay  has no effect on the stability of the system.The numerical simulations and results that prove the stability of equilibria suggest that if we want to reduce the proportions of the infective population and increase the aware population, we can increase the dissemination rate  and implementation rate .They are conducive to controlling the spread of diseases.

Figure 1 :Figure 2 :
Figure 1: The stability of  0 and  * with different initial values.

Figure 3 :
Figure 3: The stability of   ,  with variational .