# The Impact of Awareness Programs with Recruitment and Delay on the Spread of an Epidemic

Authors: Lixia Zuo; Maoxing Liu; Jinqiang Wang
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

A compartment epidemic model with delay is given to discuss the impact of awareness programs on the spread and control of infectious diseases in a given region. It is assumed that there is a constant recruitment rate in the cumulative density of awareness programs, and further it is assumed that awareness programs can influence the susceptible to a limited extent. The system exhibits two equilibria: the disease-free equilibrium is stable if the basic reproduction number is less than unity for any delay and the unique endemic equilibrium exhibits Hopf-bifurcation under certain conditions. Numerical simulations prove the results of analysis and the significance of awareness programs in preventing and controlling the diseases, by investigating the relationship between the proportion of the infective and the dissemination rate and the implementation rate, respectively.


---

## Body

## 1. Introduction

Plenty of evidence shows that awareness programs, which can influence the susceptible to a limited extent due to some objective factors, play an important role in the spread and control of infectious diseases. For example, during the outbreak of SARS, H1N1 influenza pandemic, and HIV epidemic, public media had massive reports on the number of the infections and deaths per day, which had a great impact on the diseases control [1–4]. That is due to the fact that the spread of diseases is often accompanied by a rise in awareness of those in the social vicinity of infected individuals and a subsequent change in behavior, such as keeping social distancing, wearing protective masks, and vaccination [5, 6]. Such reactions can manifest themselves in lower susceptibility as people try to prevent themselves from catching the disease, but also in lower infectivity because of self-imposed quarantine or better hygiene, shorter durations of infectiousness, or longer immunity [7]. And once the infective are cured, they will be aware of the disease [8].Recently some scholars used mathematical models to discuss the impact of awareness programs on the diseases spreading and controlling in a given region [6, 8–11]. Joshi et al. [9] formulated a model to investigate the effect on the HIV epidemic in Uganda and compared their model with three types of the susceptible to a standard SIR model and then pointed out that the awareness programs in Uganda are successful in combating diseases. Li and Cui [1] analyzed a SIS epidemic model incorporating media coverage under constant and pulse vaccination; then they obtained the exact periodic infection-free solution which is globally asymptotically stable under some conditions. In order to better describe populations mixed condition, some authors studied infectious diseases models on different networks [12–14]. In [12], the authors introduced a constant to represent the density level of media coverage from other regions. Funk et al. [14] formulated and analyzed a mathematical model in a host population; then they put forward that awareness programs can result in a lower size of the outbreak but do not affect the epidemic threshold, and if the behavioral response is treated as a local effect arising in the proximity of an outbreak, it can completely stop a disease from spreading, although only if the infection rate is below a threshold.Some scholars consider some other factors in awareness program models. Liu and Wang et al. took into account the random perturbation [15, 16]. In [16], the authors extended a deterministic SIRS epidemic model to a stochastic differential equation, and then they discussed the exponentialp-stability and global stability of unique positive solution. Some scholars focused their attention at the contact rate, and most of them assumed that awareness programs will aid in modifying the contact rate between the susceptible and the infective [17–29]. Liu and Cui [17] used the contact rate (
β
1
-
β
2
I
/
(
m
+
I
)
) and found that both of the two equilibria are asymptotically stable. Mirsa et al. [20] used the constant k with the contact rate M
/
k
+
M to limit the effect of awareness programs on the susceptible and sought out the conditions that Hopf-bifurcation occurs. Tchuenche and Bauch [23] used an exponentially decreasing function e
-
M
(
t
) to affect the force of infection; then numerical results showed the potential short-term beneficial effect of awareness programs. Pang and Cui [26] used β
(
I
)
=
μ
1
-
μ
2
f
(
I
) to show the contact rate after awareness programs alerts and found that though it is not a determined fact to eradicate the infection of the diseases, the effective awareness programs can postpone the arrival of the infection peak. Elenbaas et al. [28] introduced a segmented function to describe the media impact e
-
m
I
c when they formulate an epidemic model. A Filippov epidemic model was proposed to describe the real characteristics of media impact on the spread of infectious diseases by incorporating a piecewise continuous transmission rate β
e
(
-
α
ϵ
I
)
S
I in [29]. Mathematical and bifurcation analyses with regard to the local and global stability of equilibria and local sliding bifurcations are performed.But what we regard as unreasonable is that most of the articles assume that the cured infective become unaware of the disease. In fact people have a certain consciousness about the disease once they get sick. Therefore we propose a delayed mathematical model for predicting the future course of any epidemic by considering some of the infect join the aware susceptible after recovery. We hold the opinion that awareness programs can influence the susceptible to a limited extent for some objective factors and then consider the interaction between the susceptible and awareness programs as Holling type-II functional response. In addition we make a constant to represent the density level of media coverage from other regions with the disease because other regions can also effect the region that we consider. In fact the results about global stability of other delayed systems could be further utilized for other related problems [30–32].The rest of this paper is organized as follows. In the next section, a mathematical model with delay has been proposed to capture the dynamics of the effect of awareness programs. It is assumed that diseases spread due to the contact between the susceptible and the infective only. Then we analyze the conditions of the stability of equilibria and the existence of Hopf-bifurcation in Section3. Furthermore, in Section 4 we perform some numerical examples to validate the analysis in Section 3 and then introduce the importance of the dissemination rate, implementation rate, and the delay in disease control. In Section 5 we discuss the above contents.


## 2. Mathematical Model and Equilibrium Analysis

In the region under consideration the rate of immigration of the susceptible isb. It is assumed that the disease spreads due to the direct contact between the susceptible and the infective only and due to awareness programs the susceptible avoid being in contact with the infective and form a different class with a proportion λ, named the aware susceptible. So the total population is divided into three classes: the susceptible, the aware susceptible, and the infective, the proportions of which at time t in the total population are X
(
t
), X
m
(
t
), and Y
(
t
). It is assumed that the aware susceptible may lose awareness with passage of time and become susceptible with the proportion of λ
0 again. The infective can be cured with a proportion ν and a fraction q of recovered people will become aware and join the aware susceptible class whereas the remaining fraction p
(
p
+
q
=
1
) will join the susceptible.Consider that the cumulative density of awareness programs driven by media in that region at timet is M
(
t
), which is related to the infective. We make the constant μ represent the executed rate of awareness programs. As the time passes, some campaigns lose their impact on people and lead to the diminution of the awareness programs, so we introduce μ
0 to denote the rate of their depletion. Moreover m
0 represents the density level of media coverage on the disease from other regions. Using the fact that there is a limited extent of the awareness programs influence on the susceptible due to some objective factors, we introduce k to limit the effect and consider the interaction between them as Holling type-II functional response. It is plausible that the policy makers need some time to gather the cases of the infective generally; then we introduce τ and consider that the cumulative density of awareness programs at time t being executed will be in accordance with the infected cases reported at time t
-
τ
τ
>
0. Keeping the above facts in mind, the dynamics of model is governed by the following system of nonlinear delay differential equations: (1)
X
′
t
=
b
-
β
X
(
t
)
Y
(
t
)
+
λ
0
X
m
t
-
λ
X
(
t
)
M
t
1
+
k
M
t
+
p
ν
Y
t
-
d
X
t
,
X
m
′
t
=
λ
X
t
M
t
1
+
k
M
t
-
λ
0
X
m
(
t
)
+
q
ν
Y
(
t
)
-
d
X
m
t
,
Y
′
t
=
β
X
t
Y
t
-
ν
Y
(
t
)
-
d
Y
t
,
M
′
t
=
m
0
+
μ
Y
t
-
τ
-
μ
0
M
t
.Here X
0
=
X
0
>
0,  X
m
0
=
X
m
0
≥
0, and Y
(
θ
)
=
Y
0
≥
0 for θ
∈
[
-
τ
,
0
] and M
(
0
)
=
M
0
≥
0.In the above model, the constantsβ and d
(
d
=
b
), respectively, represent the contact rate of the unaware susceptible with the infective and the natural death rate. All the constants in the system are assumed to be positive.Using the fact thatX
(
t
)
+
Y
(
t
)
+
X
m
(
t
)
=
1, p
+
q
=
1 system (1) is reduced to the following system: (2)
X
m
′
t
=
λ
1
-
X
m
t
-
Y
t
M
t
1
+
k
M
t
-
λ
0
+
d
X
m
t
+
q
ν
Y
t
,
Y
′
t
=
β
1
-
X
m
t
-
Y
t
Y
t
-
ν
+
d
Y
t
,
M
′
t
=
m
0
+
μ
Y
t
-
τ
-
μ
0
M
t
.Now it is sufficient to study system (2) in detail rather than system (1).For the analysis of system (2), we need the region of attraction which is given by the set Ω
=
{
(
X
m
,
Y
,
M
)
∈
R
+
3
:
0
≤
X
m
,
Y
≤
1
,
0
≤
M
≤
(
m
0
+
μ
)
/
μ
0
}, and it attracts all solutions initiating in the interior of the positive orthant.The system (2) has two equilibria.(i)
Disease-free equilibriumE
0
(
λ
m
0
/
λ
m
0
+
(
λ
0
+
d
)
(
μ
0
+
k
m
0
)
,
0
,
m
0
/
μ
0
).
(ii)
Endemic equilibriumE
*
(
X
m
*
,
Y
*
,
M
*
).Define the basic reproduction numberR
0
=
β
/
(
ν
+
d
). The existence of equilibrium E
0 is trivial; then we prove the existence of E
* in detail. When R
0
>
1, in the equilibrium E
* the values of X
m
*, M
* are (3)
X
m
*
=
1
-
ν
+
d
β
-
Y
*
,
M
*
=
m
0
u
0
+
μ
μ
0
Y
*
,with Y
* satisfying the equation (4)
A
1
Y
*
2
+
A
2
Y
*
+
A
3
=
0
,where (5)
A
1
=
μ
k
β
λ
0
+
d
+
q
ν
,
A
2
=
μ
k
(
λ
0
+
d
)
(
ν
+
d
-
β
)
+
λ
μ
ν
+
d
iiiiiii
+
β
λ
0
+
d
+
q
ν
μ
0
+
k
m
0
,
A
3
=
λ
m
0
ν
+
d
+
λ
0
+
d
ν
+
d
-
β
μ
0
+
k
m
0
.Solving (4) we get (6)
Y
*
=
-
A
2
±
A
2
2
-
4
A
1
A
3
2
A
1
.We obtain A
1
>
0, and when β
λ
<
k
(
λ
0
+
d
)
(
β
-
ν
-
d
), (7)
A
3
=
λ
m
0
ν
+
d
+
λ
0
+
d
ν
+
d
-
β
μ
0
+
k
m
0
A
3
=
β
λ
m
0
1
R
0
-
1
+
β
λ
m
0
A
3
A
3
+
β
λ
0
+
d
μ
0
+
k
m
0
1
R
0
-
1
A
3
=
β
λ
m
0
1
R
0
-
1
+
β
μ
0
λ
0
+
d
1
R
0
-
1
A
3
A
3
+
β
m
0
λ
-
k
λ
0
+
d
1
-
1
R
0
<
0
,then we get Y
*
=
-
A
2
+
A
2
2
-
4
A
1
A
3
/
2
A
1 for Y
*
>
0.Remark 1.
From the expression ofY
*, it is easy to note that d
Y
*
/
d
λ
<
0 and d
Y
*
/
d
μ
<
0, which shows that the equilibrium proportion of the infective decreases as the rate of dissemination of awareness among the susceptible and the implementation rate of awareness programs increases.


## 3. Stability Analysis

In this section we present the local stability ofE
0, E
* and explore the conditions of Hopf-bifurcation by taking delay τ as a bifurcation parameter.


### 3.1. Stability of Equilibria without Delay (τ
=
0)

Theorem 2.
The disease-free equilibriumE
0 is locally asymptotically stable if R
0
<
1.Proof.
The Jacobian matrix corresponding to the system (2) when τ
=
0 is given below:(8)
J
=
-
λ
M
1
+
k
M
-
λ
0
-
d
-
λ
M
1
+
k
M
+
ν
q
λ
1
-
Y
-
X
m
1
1
+
k
M
2
-
β
Y
β
1
-
X
m
-
ν
+
d
-
2
β
Y
0
0
μ
-
μ
0
.
The characteristic equation atE
0 is of the form (9)
η
+
μ
0
η
+
λ
m
0
k
m
0
+
μ
0
+
λ
0
+
d
iiiiiiii
×
η
+
β
λ
m
0
λ
m
0
+
λ
0
+
d
μ
0
+
k
m
0
iiiiiiiiiiiiii
+
β
1
R
0
-
1
=
0
,where η is the eigenvalue. We get (10)
η
1
=
-
μ
0
<
0
,
η
2
=
-
λ
m
0
k
m
0
+
μ
0
+
λ
0
+
d
<
0
,
η
3
=
-
β
λ
m
0
λ
m
0
+
λ
0
+
d
μ
0
+
k
m
0
+
1
R
0
-
1
,so η
3
<
0 when R
0
<
1. Thus E
0 is locally asymptotically stable if R
0
<
1.The form of characteristic equation atE
* is (11)
η
3
+
Q
1
η
2
+
Q
2
η
+
Q
3
+
Q
4
=
0
,where (12)
Q
1
=
B
1
+
B
3
+
μ
0
,
Q
2
=
μ
0
B
1
+
B
3
+
B
1
B
3
-
β
Y
*
B
2
,
Q
3
=
μ
0
B
1
B
3
-
β
μ
0
Y
*
B
2
,
Q
4
=
-
β
μ
Y
*
B
4
,
B
1
=
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
,
B
2
=
λ
M
*
1
1
+
k
M
*
-
ν
q
,
B
3
=
β
X
m
*
-
1
+
2
β
Y
*
+
ν
+
d
,
B
4
=
λ
X
m
*
+
Y
*
-
1
1
1
+
k
M
*
2
.Theorem 3.
When the endemic equilibriumE
* exists, it is locally asymptotically stable provided m
0
(
λ
0
+
d
+
q
ν
)
>
μ
μ
0.Proof.
For the characteristic equation(13)
η
3
+
Q
1
η
2
+
Q
2
η
+
Q
3
+
Q
4
=
0
,
it is easy to show(14)
Q
1
=
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
+
β
X
m
*
-
1
+
2
β
Y
*
+
ν
+
d
+
μ
0
=
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
+
β
X
m
*
-
1
+
2
β
Y
*
+
β
1
-
X
m
*
-
Y
*
+
μ
0
=
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
+
β
Y
*
+
μ
0
>
0
,
Q
2
=
μ
0
B
1
+
B
3
+
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
×
β
X
m
*
-
1
+
2
β
Y
*
+
ν
+
d
-
β
Y
*
λ
M
*
1
1
+
k
M
*
-
ν
q
=
μ
0
B
1
+
B
3
+
β
Y
*
λ
M
*
1
+
k
M
*
+
λ
0
+
d
-
β
Y
*
λ
M
*
1
+
k
M
*
-
ν
q
=
μ
0
B
1
+
B
3
+
λ
0
+
ν
q
+
d
β
Y
*
>
0
,
Q
3
+
Q
4
=
μ
0
B
1
B
3
-
β
Y
*
μ
B
4
+
μ
0
B
2
=
μ
0
β
Y
*
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
-
β
Y
*
μ
λ
X
m
*
+
Y
*
-
1
1
1
+
k
M
*
2
iiiiiiiiiiii
+
μ
0
(
λ
M
*
1
1
+
k
M
*
-
ν
q
)
=
β
Y
*
μ
λ
1
-
X
m
*
-
Y
*
1
1
+
k
M
*
2
+
μ
0
β
Y
*
λ
0
+
d
+
ν
q
μ
0
β
Y
*
>
0
.And (15)
Q
1
Q
2
-
Q
3
+
Q
4
=
B
1
+
B
3
μ
0
B
1
+
μ
0
B
3
+
B
1
B
3
-
β
Y
*
B
2
+
μ
0
2
B
1
+
B
3
-
Q
4
=
μ
0
B
1
+
B
3
2
+
μ
0
2
B
1
+
B
3
+
B
3
B
1
+
B
3
B
1
-
B
2
+
μ
B
4
,where (16)
B
1
+
B
3
B
1
-
B
2
+
μ
B
4
=
λ
M
*
1
+
k
M
*
+
λ
0
+
d
+
β
Y
*
λ
0
+
d
+
q
ν
-
μ
λ
1
-
X
m
*
-
Y
*
1
1
+
k
M
*
2
=
λ
1
+
k
M
*
M
*
λ
0
+
d
+
q
ν
-
μ
1
+
k
M
*
+
λ
0
+
d
+
β
Y
*
λ
0
+
d
+
q
ν
+
μ
λ
X
m
*
+
Y
*
1
+
k
M
*
=
λ
1
+
k
M
*
2
M
*
λ
0
+
d
+
q
ν
+
k
M
*
2
×
λ
0
+
d
+
q
ν
-
μ
+
λ
0
+
d
+
β
Y
*
λ
0
+
d
+
q
ν
+
μ
λ
X
m
*
+
Y
*
1
1
+
k
M
*
=
λ
1
+
k
M
*
2
m
0
λ
0
+
d
+
q
ν
μ
0
+
μ
Y
*
λ
0
+
d
+
q
ν
μ
0
+
k
M
*
2
λ
0
+
d
+
q
ν
-
μ
m
0
λ
0
+
d
+
q
ν
μ
0
+
λ
0
+
d
+
β
Y
*
λ
0
+
d
+
q
ν
+
μ
λ
X
m
*
+
Y
*
1
1
+
k
M
*
>
0
,so Q
1
Q
2
-
(
Q
3
+
Q
4
)
>
0. According to Hurwitz criterion, we can claim that all the eigenvalues will be either negative or having negative real part. Thus, the equilibrium E
* is locally asymptotically stable.


### 3.2.  Stability of Equilibria with Delay (τ
>
0)

Theorem 4.
The disease-free equilibriumE
0 is locally asymptotically stable if R
0
<
1.Similar to the proof of Theorem2, this proof is omitted.We linearize system (2) about E
* to study the stability of the endemic equilibrium and get (17)
d
u
d
t
=
M
1
u
t
+
M
2
u
t
-
τ
,where(18)
u
(
t
)
=
x
t
y
t
z
t
T
,
1
M
1
=
-
λ
M
*
1
+
k
M
*
-
λ
0
-
d
-
λ
M
*
1
+
k
M
*
+
ν
q
λ
1
-
Y
*
-
X
m
*
1
1
+
k
M
*
2
-
β
Y
*
β
1
-
X
m
*
-
ν
+
d
-
2
β
Y
*
0
0
0
-
μ
0
,
1
M
2
=
0
0
0
0
0
0
0
μ
0
.In the above, x, y, and z are small perturbations around E
*.Then the form of the characteristic equation of the system is(19)
η
3
+
Q
1
η
2
+
Q
2
η
+
Q
3
=
Q
4
e
-
η
τ
,where η is the eigenvalue.To show the Hopf-bifurcation, we need to show that (19) has a pair of purely imaginary roots. For this, substituting η
=
i
ω
(
ω
>
0
) into (19) and separating real and imaginary parts, we get the following transcendental equations: (20)
Q
2
ω
-
ω
3
=
Q
4
sin
ω
τ
,
Q
1
ω
2
-
Q
3
=
Q
4
cos
⁡
ω
τ
.Squaring and adding the above equations and substituting ψ
=
ω
2, we get (21)
h
(
ψ
)
=
ψ
3
+
P
1
ψ
2
+
P
2
ψ
+
P
3
=
0
,where P
1
=
Q
1
2
-
2
Q
2, P
2
=
Q
2
2
-
2
Q
1
Q
3, and P
3
=
Q
3
2
-
Q
4
2.If the coefficients inh
(
ψ
) satisfy the conditions of Routh-Hurwitz criterion, then E
* is locally asymptotically stable for all delays τ
>
0, provided it is stable in absence of delay. In the following we consider that the values of P
i
(
i
=
1
,
2
,
3
,
4
) do not satisfy the Routh-Hurwitz criterion.Lemma 5.
IfB
2
>
0
,
λ
0
+
d
+
q
ν
<
μ
λ
(
ν
+
d
)
(
1
/
β
(
1
+
k
M
*
)
2
)
, (21) will have at least one positive root ψ
0.Proof.
Forh
(
ψ
)
=
ψ
3
+
P
1
ψ
2
+
P
2
ψ
+
P
3, it is obvious that (22)
P
1
=
B
1
+
B
3
+
μ
0
2
-
2
B
1
+
B
3
μ
0
+
B
1
B
3
-
β
B
2
Y
*
P
1
=
B
1
2
+
B
3
2
+
μ
0
2
+
2
β
B
2
Y
*
>
0
,
P
2
=
B
1
+
B
3
μ
0
+
B
1
B
3
-
β
B
2
Y
*
2
-
2
B
1
+
B
3
+
μ
0
B
1
B
3
μ
0
-
μ
0
β
B
2
Y
*
P
1
=
B
1
2
μ
0
2
+
B
3
2
μ
0
2
+
B
1
B
3
2
+
B
2
B
3
2
-
2
B
1
B
2
B
3
2
+
2
B
2
B
3
μ
0
2
P
1
=
B
1
2
+
B
3
2
μ
0
2
+
B
1
B
3
-
B
2
B
3
2
+
2
μ
0
2
B
2
B
3
>
0
,
P
3
=
Q
3
+
Q
4
Q
3
-
Q
4
P
1
=
B
3
μ
0
B
1
-
μ
0
B
2
+
μ
B
4
Q
3
+
Q
4
P
1
=
B
3
Q
3
+
Q
4
i
i
i
i
i
i
i
×
1
β
1
+
k
M
*
2
μ
0
λ
0
+
d
+
q
ν
iiiiiiiiiiiiii
-
μ
λ
1
-
X
m
*
-
Y
*
1
1
+
k
M
*
2
P
1
<
B
3
Q
3
+
Q
4
1
β
1
+
k
M
*
2
λ
0
+
d
+
q
ν
iiiiiiiiiiiiiiiiiiiiiiiiiiiii
-
μ
λ
ν
+
d
1
β
1
+
k
M
*
2
<
0
.Then we have h
(
0
)
<
0 and h
(
+
∞
)
→
+
∞; thus (21) has at least one positive root ψ
0.Denoteω
0
=
ψ
0; then (21) has a pair of purely imaginary roots (
±
i
ω
0
).Now we turn to the bifurcation analysis. We use the delayτ as the bifurcation parameter. We view the solutions of (21) as a function of the bifurcation parameter τ, and let η
(
τ
)
=
γ
(
τ
)
+
i
ω
(
τ
) be the eigenvalue of (21) such that, for the initial value of the bifurcation parameter τ
0, we have γ
(
τ
0
)
=
0 and ω
(
τ
0
)
=
ω
0 (ω
0
>
0). To establish the Hopf bifurcation at τ
=
τ
0, we need to show that d
Re
η
(
τ
)
/
d
τ
|
τ
=
τ
0
>
0.Lemma 6.
One has the following transversality condition:(23)
d
Re
η
τ
d
τ
τ
=
τ
0
>
0
.Proof.
Differentiating with respect toτ from (21), we get (24)
d
η
d
τ
-
1
=
d
τ
d
η
=
3
η
2
+
2
Q
1
η
+
Q
2
-
Q
4
τ
e
-
η
τ
Q
4
η
e
-
η
τ
.So (25)
Sign
Re
d
η
d
τ
τ
=
τ
0
=
Sign
Re
d
τ
d
η
η
=
i
ω
0
=
Sign
Re
3
η
2
+
2
Q
1
η
+
Q
2
-
Q
4
τ
e
-
η
τ
Q
4
η
e
-
η
τ
η
=
i
ω
0
=
Sign
Re
3
η
2
+
2
Q
1
η
+
Q
2
Q
4
η
e
-
η
τ
η
=
i
ω
0
=
Sign
-
Re
τ
η
η
=
i
ω
0
=
Sign
Re
-
3
η
2
+
2
Q
1
η
+
Q
2
η
4
+
Q
1
η
3
+
Q
2
η
2
+
Q
3
η
η
=
i
ω
0
=
Sign
Re
-
Q
2
-
3
ω
0
2
+
2
Q
1
ω
0
i
ω
0
4
-
Q
2
ω
0
2
+
Q
3
ω
0
-
Q
1
ω
0
3
i
=
Sign
ω
0
2
3
ω
0
4
+
2
Q
1
2
-
4
Q
2
ω
0
2
+
Q
2
2
-
2
Q
1
Q
3
Q
4
2
ω
0
4
=
Sign
3
ω
0
4
+
2
P
1
ω
0
2
+
P
2
Q
4
2
>
0
.This proves the lemma. Now we have the following theorem.Theorem 7.
The endemic equilibriumE
* of the system is locally asymptotically stable when τ
<
τ
0 and becomes unstable for τ
>
τ
0 provided (26)
μ
μ
0
m
0
<
λ
0
+
d
+
q
ν
<
min
⁡
λ
M
*
1
+
k
M
*
,
μ
λ
ν
+
d
β
1
+
k
M
*
2
.When τ
=
τ
0, a Hopf-bifurcation occurs, leading a family of periodic solutions bifurcating from E
* as τ passes through the critical value τ
0, where (27)
τ
0
=
1
ω
0
arc
cos
⁡
Q
1
ω
0
2
-
Q
3
Q
4
.


## 3.1. Stability of Equilibria without Delay (τ
=
0)

Theorem 2.
The disease-free equilibriumE
0 is locally asymptotically stable if R
0
<
1.Proof.
The Jacobian matrix corresponding to the system (2) when τ
=
0 is given below:(8)
J
=
-
λ
M
1
+
k
M
-
λ
0
-
d
-
λ
M
1
+
k
M
+
ν
q
λ
1
-
Y
-
X
m
1
1
+
k
M
2
-
β
Y
β
1
-
X
m
-
ν
+
d
-
2
β
Y
0
0
μ
-
μ
0
.
The characteristic equation atE
0 is of the form (9)
η
+
μ
0
η
+
λ
m
0
k
m
0
+
μ
0
+
λ
0
+
d
iiiiiiii
×
η
+
β
λ
m
0
λ
m
0
+
λ
0
+
d
μ
0
+
k
m
0
iiiiiiiiiiiiii
+
β
1
R
0
-
1
=
0
,where η is the eigenvalue. We get (10)
η
1
=
-
μ
0
<
0
,
η
2
=
-
λ
m
0
k
m
0
+
μ
0
+
λ
0
+
d
<
0
,
η
3
=
-
β
λ
m
0
λ
m
0
+
λ
0
+
d
μ
0
+
k
m
0
+
1
R
0
-
1
,so η
3
<
0 when R
0
<
1. Thus E
0 is locally asymptotically stable if R
0
<
1.The form of characteristic equation atE
* is (11)
η
3
+
Q
1
η
2
+
Q
2
η
+
Q
3
+
Q
4
=
0
,where (12)
Q
1
=
B
1
+
B
3
+
μ
0
,
Q
2
=
μ
0
B
1
+
B
3
+
B
1
B
3
-
β
Y
*
B
2
,
Q
3
=
μ
0
B
1
B
3
-
β
μ
0
Y
*
B
2
,
Q
4
=
-
β
μ
Y
*
B
4
,
B
1
=
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
,
B
2
=
λ
M
*
1
1
+
k
M
*
-
ν
q
,
B
3
=
β
X
m
*
-
1
+
2
β
Y
*
+
ν
+
d
,
B
4
=
λ
X
m
*
+
Y
*
-
1
1
1
+
k
M
*
2
.Theorem 3.
When the endemic equilibriumE
* exists, it is locally asymptotically stable provided m
0
(
λ
0
+
d
+
q
ν
)
>
μ
μ
0.Proof.
For the characteristic equation(13)
η
3
+
Q
1
η
2
+
Q
2
η
+
Q
3
+
Q
4
=
0
,
it is easy to show(14)
Q
1
=
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
+
β
X
m
*
-
1
+
2
β
Y
*
+
ν
+
d
+
μ
0
=
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
+
β
X
m
*
-
1
+
2
β
Y
*
+
β
1
-
X
m
*
-
Y
*
+
μ
0
=
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
+
β
Y
*
+
μ
0
>
0
,
Q
2
=
μ
0
B
1
+
B
3
+
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
×
β
X
m
*
-
1
+
2
β
Y
*
+
ν
+
d
-
β
Y
*
λ
M
*
1
1
+
k
M
*
-
ν
q
=
μ
0
B
1
+
B
3
+
β
Y
*
λ
M
*
1
+
k
M
*
+
λ
0
+
d
-
β
Y
*
λ
M
*
1
+
k
M
*
-
ν
q
=
μ
0
B
1
+
B
3
+
λ
0
+
ν
q
+
d
β
Y
*
>
0
,
Q
3
+
Q
4
=
μ
0
B
1
B
3
-
β
Y
*
μ
B
4
+
μ
0
B
2
=
μ
0
β
Y
*
λ
M
*
1
1
+
k
M
*
+
λ
0
+
d
-
β
Y
*
μ
λ
X
m
*
+
Y
*
-
1
1
1
+
k
M
*
2
iiiiiiiiiiii
+
μ
0
(
λ
M
*
1
1
+
k
M
*
-
ν
q
)
=
β
Y
*
μ
λ
1
-
X
m
*
-
Y
*
1
1
+
k
M
*
2
+
μ
0
β
Y
*
λ
0
+
d
+
ν
q
μ
0
β
Y
*
>
0
.And (15)
Q
1
Q
2
-
Q
3
+
Q
4
=
B
1
+
B
3
μ
0
B
1
+
μ
0
B
3
+
B
1
B
3
-
β
Y
*
B
2
+
μ
0
2
B
1
+
B
3
-
Q
4
=
μ
0
B
1
+
B
3
2
+
μ
0
2
B
1
+
B
3
+
B
3
B
1
+
B
3
B
1
-
B
2
+
μ
B
4
,where (16)
B
1
+
B
3
B
1
-
B
2
+
μ
B
4
=
λ
M
*
1
+
k
M
*
+
λ
0
+
d
+
β
Y
*
λ
0
+
d
+
q
ν
-
μ
λ
1
-
X
m
*
-
Y
*
1
1
+
k
M
*
2
=
λ
1
+
k
M
*
M
*
λ
0
+
d
+
q
ν
-
μ
1
+
k
M
*
+
λ
0
+
d
+
β
Y
*
λ
0
+
d
+
q
ν
+
μ
λ
X
m
*
+
Y
*
1
+
k
M
*
=
λ
1
+
k
M
*
2
M
*
λ
0
+
d
+
q
ν
+
k
M
*
2
×
λ
0
+
d
+
q
ν
-
μ
+
λ
0
+
d
+
β
Y
*
λ
0
+
d
+
q
ν
+
μ
λ
X
m
*
+
Y
*
1
1
+
k
M
*
=
λ
1
+
k
M
*
2
m
0
λ
0
+
d
+
q
ν
μ
0
+
μ
Y
*
λ
0
+
d
+
q
ν
μ
0
+
k
M
*
2
λ
0
+
d
+
q
ν
-
μ
m
0
λ
0
+
d
+
q
ν
μ
0
+
λ
0
+
d
+
β
Y
*
λ
0
+
d
+
q
ν
+
μ
λ
X
m
*
+
Y
*
1
1
+
k
M
*
>
0
,so Q
1
Q
2
-
(
Q
3
+
Q
4
)
>
0. According to Hurwitz criterion, we can claim that all the eigenvalues will be either negative or having negative real part. Thus, the equilibrium E
* is locally asymptotically stable.


## 3.2.  Stability of Equilibria with Delay (τ
>
0)

Theorem 4.
The disease-free equilibriumE
0 is locally asymptotically stable if R
0
<
1.Similar to the proof of Theorem2, this proof is omitted.We linearize system (2) about E
* to study the stability of the endemic equilibrium and get (17)
d
u
d
t
=
M
1
u
t
+
M
2
u
t
-
τ
,where(18)
u
(
t
)
=
x
t
y
t
z
t
T
,
1
M
1
=
-
λ
M
*
1
+
k
M
*
-
λ
0
-
d
-
λ
M
*
1
+
k
M
*
+
ν
q
λ
1
-
Y
*
-
X
m
*
1
1
+
k
M
*
2
-
β
Y
*
β
1
-
X
m
*
-
ν
+
d
-
2
β
Y
*
0
0
0
-
μ
0
,
1
M
2
=
0
0
0
0
0
0
0
μ
0
.In the above, x, y, and z are small perturbations around E
*.Then the form of the characteristic equation of the system is(19)
η
3
+
Q
1
η
2
+
Q
2
η
+
Q
3
=
Q
4
e
-
η
τ
,where η is the eigenvalue.To show the Hopf-bifurcation, we need to show that (19) has a pair of purely imaginary roots. For this, substituting η
=
i
ω
(
ω
>
0
) into (19) and separating real and imaginary parts, we get the following transcendental equations: (20)
Q
2
ω
-
ω
3
=
Q
4
sin
ω
τ
,
Q
1
ω
2
-
Q
3
=
Q
4
cos
⁡
ω
τ
.Squaring and adding the above equations and substituting ψ
=
ω
2, we get (21)
h
(
ψ
)
=
ψ
3
+
P
1
ψ
2
+
P
2
ψ
+
P
3
=
0
,where P
1
=
Q
1
2
-
2
Q
2, P
2
=
Q
2
2
-
2
Q
1
Q
3, and P
3
=
Q
3
2
-
Q
4
2.If the coefficients inh
(
ψ
) satisfy the conditions of Routh-Hurwitz criterion, then E
* is locally asymptotically stable for all delays τ
>
0, provided it is stable in absence of delay. In the following we consider that the values of P
i
(
i
=
1
,
2
,
3
,
4
) do not satisfy the Routh-Hurwitz criterion.Lemma 5.
IfB
2
>
0
,
λ
0
+
d
+
q
ν
<
μ
λ
(
ν
+
d
)
(
1
/
β
(
1
+
k
M
*
)
2
)
, (21) will have at least one positive root ψ
0.Proof.
Forh
(
ψ
)
=
ψ
3
+
P
1
ψ
2
+
P
2
ψ
+
P
3, it is obvious that (22)
P
1
=
B
1
+
B
3
+
μ
0
2
-
2
B
1
+
B
3
μ
0
+
B
1
B
3
-
β
B
2
Y
*
P
1
=
B
1
2
+
B
3
2
+
μ
0
2
+
2
β
B
2
Y
*
>
0
,
P
2
=
B
1
+
B
3
μ
0
+
B
1
B
3
-
β
B
2
Y
*
2
-
2
B
1
+
B
3
+
μ
0
B
1
B
3
μ
0
-
μ
0
β
B
2
Y
*
P
1
=
B
1
2
μ
0
2
+
B
3
2
μ
0
2
+
B
1
B
3
2
+
B
2
B
3
2
-
2
B
1
B
2
B
3
2
+
2
B
2
B
3
μ
0
2
P
1
=
B
1
2
+
B
3
2
μ
0
2
+
B
1
B
3
-
B
2
B
3
2
+
2
μ
0
2
B
2
B
3
>
0
,
P
3
=
Q
3
+
Q
4
Q
3
-
Q
4
P
1
=
B
3
μ
0
B
1
-
μ
0
B
2
+
μ
B
4
Q
3
+
Q
4
P
1
=
B
3
Q
3
+
Q
4
i
i
i
i
i
i
i
×
1
β
1
+
k
M
*
2
μ
0
λ
0
+
d
+
q
ν
iiiiiiiiiiiiii
-
μ
λ
1
-
X
m
*
-
Y
*
1
1
+
k
M
*
2
P
1
<
B
3
Q
3
+
Q
4
1
β
1
+
k
M
*
2
λ
0
+
d
+
q
ν
iiiiiiiiiiiiiiiiiiiiiiiiiiiii
-
μ
λ
ν
+
d
1
β
1
+
k
M
*
2
<
0
.Then we have h
(
0
)
<
0 and h
(
+
∞
)
→
+
∞; thus (21) has at least one positive root ψ
0.Denoteω
0
=
ψ
0; then (21) has a pair of purely imaginary roots (
±
i
ω
0
).Now we turn to the bifurcation analysis. We use the delayτ as the bifurcation parameter. We view the solutions of (21) as a function of the bifurcation parameter τ, and let η
(
τ
)
=
γ
(
τ
)
+
i
ω
(
τ
) be the eigenvalue of (21) such that, for the initial value of the bifurcation parameter τ
0, we have γ
(
τ
0
)
=
0 and ω
(
τ
0
)
=
ω
0 (ω
0
>
0). To establish the Hopf bifurcation at τ
=
τ
0, we need to show that d
Re
η
(
τ
)
/
d
τ
|
τ
=
τ
0
>
0.Lemma 6.
One has the following transversality condition:(23)
d
Re
η
τ
d
τ
τ
=
τ
0
>
0
.Proof.
Differentiating with respect toτ from (21), we get (24)
d
η
d
τ
-
1
=
d
τ
d
η
=
3
η
2
+
2
Q
1
η
+
Q
2
-
Q
4
τ
e
-
η
τ
Q
4
η
e
-
η
τ
.So (25)
Sign
Re
d
η
d
τ
τ
=
τ
0
=
Sign
Re
d
τ
d
η
η
=
i
ω
0
=
Sign
Re
3
η
2
+
2
Q
1
η
+
Q
2
-
Q
4
τ
e
-
η
τ
Q
4
η
e
-
η
τ
η
=
i
ω
0
=
Sign
Re
3
η
2
+
2
Q
1
η
+
Q
2
Q
4
η
e
-
η
τ
η
=
i
ω
0
=
Sign
-
Re
τ
η
η
=
i
ω
0
=
Sign
Re
-
3
η
2
+
2
Q
1
η
+
Q
2
η
4
+
Q
1
η
3
+
Q
2
η
2
+
Q
3
η
η
=
i
ω
0
=
Sign
Re
-
Q
2
-
3
ω
0
2
+
2
Q
1
ω
0
i
ω
0
4
-
Q
2
ω
0
2
+
Q
3
ω
0
-
Q
1
ω
0
3
i
=
Sign
ω
0
2
3
ω
0
4
+
2
Q
1
2
-
4
Q
2
ω
0
2
+
Q
2
2
-
2
Q
1
Q
3
Q
4
2
ω
0
4
=
Sign
3
ω
0
4
+
2
P
1
ω
0
2
+
P
2
Q
4
2
>
0
.This proves the lemma. Now we have the following theorem.Theorem 7.
The endemic equilibriumE
* of the system is locally asymptotically stable when τ
<
τ
0 and becomes unstable for τ
>
τ
0 provided (26)
μ
μ
0
m
0
<
λ
0
+
d
+
q
ν
<
min
⁡
λ
M
*
1
+
k
M
*
,
μ
λ
ν
+
d
β
1
+
k
M
*
2
.When τ
=
τ
0, a Hopf-bifurcation occurs, leading a family of periodic solutions bifurcating from E
* as τ passes through the critical value τ
0, where (27)
τ
0
=
1
ω
0
arc
cos
⁡
Q
1
ω
0
2
-
Q
3
Q
4
.


## 4. Numerical Simulations and Results

To check the feasibility of our analysis in Section3 about E
0, we present some numerical computations in this section using MatLab by choosing the following set of parameter values: d
=
0.002,  β
=
0.4, μ
0
=
0.08, μ
=
0.2
,
ν
=
0.46, q
=
0.6,  λ
=
0.2, λ
0
=
0.4, m
0
=
0.04, k
=
0.5, and τ
=
10 (see Figure 1), and then we get Figure 2 with different initial values, from which we can know that E
0 is stable and all the trajectories approach it. The numerical simulations support the analysis given.Figure 1
Schematic model flow diagram.Figure 2
The stability of the disease-free equilibriumE
0 with τ
=
10.Then we choose the following set of parameter values which satisfy the condition in Theorem7: d
=
0.00004,  β
=
0.3,  μ
0
=
0.14,  μ
=
0.35,  ν
=
0.2,  q
=
0.15,  λ
=
0.019,  λ
0
=
0.001,  m
0
=
0.000005, and  k
=
0.0003. The numerical value of τ
0 computed is found to be 30.12. When τ
>
0, we give different τ as follows: τ
=
20
,
τ
=
25
,
τ
=
35
,
and
τ
=
40 (Figures 3 and 4). As shown in Figures 3 and 4 the variables approach their equilibria when τ is less than τ
0, whereas, as τ exceeds its critical value τ
0, all variables start showing oscillatory behavior. This indicates that, in the latter case, sometimes the infective will be high and sometimes will be low and it may be difficult to make the prediction regarding the size of epidemic. It is clear that τ plays a key role in the prevention and control of diseases.Figure 3
The stability of the endemic equilibriaE
* with τ
=
20 and τ
=
40, respectively.Figure 4
The stability of the endemic equilibriaE
* with τ
=
25 and τ
=
35, respectively.In the following we letτ
=
10 and make λ, μ change from 0.05 to 0.8 and research the influence of the dissemination rate λ and the implementation rate μ on the infective Y
(
t
) separately. The variation of Y
(
t
) with respect to time t for different values of λ and μ is shown in Figure 5, and we discover all of the proportions of Y
* reductions as λ and μ increase, which proves the conclusions of the remark, but λ has greater influence on Y
(
t
) than μ. In addition the reason why Y
(
t
) has a similar trend in Figure 5 is that λ and μ have a similar trend with M
(
t
). There really is an effort here to make clear that λ and μ play a key role in the prevention and control of diseases.Figure 5
Whenτ
=
10, the stability of Y
* with variational λ and μ, respectively.A comparison between the oscillations inY and M is presented in Figure 6. From this figure, we obtain that, as the infective increase, the awareness programs also start growing with some time lag due to time delay in the execution of awareness programs. As soon as the awareness programs reach a potentially high value, the infective start decreasing due to the execution of awareness programs in the infective. Awareness programs after certain time also start decreasing, which eventually results in the increase of the infective. Both of the figures show the interplay between Y and M, which proves the existence of delay.Figure 6
Comparison between variations ofY and M with time.


## 5. Discussion

In this paper, a nonlinear mathematical model with delay and awareness programs driven by media has been proposed and analyzed. It is assumed that pathogens are transmitted via direct contact between the susceptible and the infective. The model exhibits two equilibria, and the disease-free equilibrium has been shown to be stable for basic reproduction numberR
0
<
1 when τ
≥
0. For R
0
>
1, it leads to the existence and stability of an endemic equilibrium under some conditions in absence of time delay. But for τ
>
0 this equilibrium is locally asymptotically stable when the delay is suitably small; that is, τ
<
τ
0, while a loss of stability by a Hopf-bifurcation can occur as delay increases. Numerical simulations prove the stability of equilibria and show that the value of delay sustainable for the disease is 30.12. The numerical results suggest that if we want to reduce the proportion of the infective, we can increase the dissemination rate λ and implementation rate μ. At last we obtain the existence of τ from Figure 6. In short media is a key tool for influencing people behavior towards the disease to devise proper policies for controlling the epidemic.



---

*Source: 235935-2015-03-19.xml*