# Successive Vaccination and Difference in Immunity of a Delay SIR Model with a General Incidence Rate

Authors: Yongzhen Pei; Li Changguo; Qianyong Wu; Yunfei Lv
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

A delay SIR epidemic model with difference in immunity and successive vaccination is proposed to understand their effects on the disease spread. From theorems, it is obtained that the basic reproduction number governs the dynamic behavior of the system. The existence and stability of the possible equilibria are examined in terms of a certain threshold condition about the basic reproduction number. By use of new computational techniques for delay differential equations, we prove that the system is permanent. Our results indicate that the recovery rate and the vaccination rate are two factors for the dynamic behavior of the system. Numerical simulations are carried out to investigate the influence of the key parameters on the spread of the disease, to support the analytical conclusion, and to illustrate possible behavioral scenarios of the model.


---

## Body

## 1. Introduction

The current threat of some new type diseases has raised our awareness that curbing the spread of some emerging and reemerging human diseases is of public health importance such as H1N1. This emerging disease, which was first reported in Mexico, spread very quickly, due to the travel of infected persons by airplanes, trains, and buses to some other regions. It continued to spread around the world and caused about 5000 deaths. In recent years, many mathematical models have been developed for the transmission dynamics of infectious diseases such as SARS, HIV/AIDS, measles, and smallpox ([1–6], to name a few, and the references therein). These models have provided understanding of the underlying mechanisms which influence the spread of diseases and suggested some control strategies. Moreover, to our knowledge, the first effective control strategy for the elimination of infectious diseases is obtaining immunity. It has been reported that “People’s immunity to A/H1N1 flu virus is greater than previously thought after access vaccines. The WHO is working to give more nations access vaccines to fight the H1N1 flu pandemic.”Besides these above studies, many authors formulated and analyzed SIR epidemic models for the control of diseases [7–20]. In particular, some authors have studied the effects of vaccination on the spread of diseases [7–11]; others have studied the effects of treatment on the spread of diseases [12–15]. Gao et al. have proposed an epidemic model with density-dependent birth pulses and seasonal prevention [16]. Recently, some works have investigated permanent and temporary immunity [17–20]. However, in these SIR models, an unrealistic assumption is that all the rest of infected individuals acquire immunity besides death. Measles encephalitis in adults in [21, 22] shows that there is difference in immunity of infected individuals. That is, some infected individuals can acquire immunity after recovery, but some do not acquire immunity and can be infected once more. At the same time, vaccination is an important strategy for the elimination of infectious diseases [7–11].Vaccinations have many types; impulsive vaccination and successive vaccination are two main policies. Successive vaccination is that people have been vaccinated at birth to protect themselves from disease; the studies can be found in [23, 24]. Makinde in [23] studied a SIR model for the transmission dynamics of a childhood disease in the presence of a preventive vaccine and analyzed the vaccination reproductive number for disease control and eradication qualitatively. Impulsive vaccination (only at fixed time sequence we execute effectively the vaccination for the disease) is an important and effective strategy for the elimination of infectious diseases and has been studied in the literature. For example, see [10, 17, 25–27]. In above-mentioned papers, authors almost considered the vaccination of susceptible population. But, in fact, under a certain situation, the vaccine treatment also should be considered for the newborns of the susceptible, the exposed, and the removed. We find that there are few studies on the aspect of the vaccination of newborns. In this case, successive vaccination seems more reasonable than impulsive vaccination. Therefore, in this paper, a SIR model with difference in immunity and successive vaccination is considered.As far as disease transmission is concerned, the incidence rate, defined as the rate of new infection, plays a very important role in modelling infectious diseases. Bilinear incidence rateβ
S
I in [28, 29] and standard incidence rates λ
S
I
/
N in [30, 31] have often been used in epidemic models. However, it is unreasonable to consider the bilinear incidence rate (based on the law of mass action) as the number of susceptibles is large, owing to the number of susceptibles with which every infective contact within a certain time is limited. Standard incidence rate may be a good approximation if the number of available partners is large enough but it is not possible to make more contacts when the population N is small. Combine the two previous approaches by assuming that if the number of available partners N is low, the number of actual per capita partners is proportional to N, whereas if the number of available partners is large, there is a saturation effect which makes the number of actual partners constant. Considering this case, a saturation incidence rate of type f
(
I
)
S with f
(
I
)
=
k
I
/
(
1
+
α
I
) is being proposed in [32]. More general incidence rate used in the literature is the one for which f
(
I
)
=
k
I
l
/
(
1
+
α
I
h
) [33, 34], where I
l measures the infection force of the disease and f
(
I
)
=
1
/
(
1
+
α
I
h
) measures the inhibitory effect caused by behavioral changes. Note that if f
(
I
) is decreasing when I is large, this may be interpreted as the fact that susceptibles tend to reduce their social contacts if the perceived number of infectives increases over a psychologically significant value. The above saturation incidence rate depends also on the size of the infectives I termed as infectives-dependent. Particular examples of susceptibles-dependent incidence rate are f
(
S
)
=
k
S
/
(
1
+
α
S
) [35]. Very general incidence rates which are not linear in S are also used in Derrick and van den Driessche [36] (f
(
S
,
I
,
N
)
=
I
Φ
(
S
,
I
,
N
), where N
=
S
+
I), Korobeinikov and Maini [37] (f
(
S
,
I
)
=
h
1
(
I
)
h
2
(
S
)), and Moghadas and Alexander [38] (f
(
S
,
I
)
=
β
(
1
+
g
(
I
,
v
)
)
I
S).Nie et al. [19] and Ji et al. [39] respectively considered a delay SIR epidemic model with nonlinear incidence rate and density-dependent birth and death rates. Motivated by the main idea described in [6, 39], in this paper, we consider a delay SIR model with difference in immunity and successive vaccination and an abstract incidence rate. The main difference between our study and those described in [6, 39] is the difference in immunity and successive vaccination and an abstract incidence rate. An abstract incidence rate of type f
(
I
)
S is employed to model the spread of the disease which is propagated through the infective individuals, under a few biologically feasible assumptions upon f
(
I
).In view of above facts, we will formulate a mathematical model in Section2. We provide the region of biologically feasible solutions in Section 3. Then, we study the existence and stability of the steady states in the next section, analyze the permanence result in Section 5, and give some numerical simulations in Section 6. Lastly, we end the paper with a brief discussion of our results in Section 6.


## 2. Model Formulation and Invariant Region

In this section, we will present a delay SIR epidemic model with a general nonlinear incidence rate. The total populationN
(
t
) is divided into three subclasses, namely, the susceptibles S
(
t
), the infectives I
(
t
), and the recovered individuals R
(
t
). Based on the SIR model in [12, 39], we considered following system:
(1)
d
S
d
t
=
[
b
-
a
r
N
K
]
N
-
β
e
-
d
1
τ
S
f
(
I
(
t
-
τ
)
)
-
[
d
+
(
1
-
a
)
r
N
K
]
S
-
θ
S
+
μ
1
I
,
d
I
d
t
=
β
e
-
d
1
τ
f
(
I
(
t
-
τ
)
)
S
-
[
d
+
(
1
-
a
)
r
N
K
]
I
-
μ
1
I
-
e
-
d
1
ω
δ
I
(
t
-
ω
)
,
d
R
d
t
=
e
-
d
1
ω
δ
I
(
t
-
ω
)
-
[
d
+
(
1
-
a
)
r
N
K
]
R
+
θ
S
,
h
h
h
h
h
h
h
h
h
h
h
N
(
0
)
=
S
0
>
0
,
I
(
0
)
=
φ
(
θ
)
≥
0
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
∀
θ
∈
[
-
τ
¯
,
0
]
,
R
(
0
)
=
R
0
≥
0
,
where τ
¯
=
max
⁡
{
τ
,
ω
} and φ
∈
C
(
[
-
τ
¯
,
0
]
,
R
). We give the following useful assumptions.(1)
There are no disease induced deaths, and all the newborns are susceptible.
(2)
f
(
I
) is the nonlinear incidence rate satisfying the following assumptions:
(2)
f
(
0
)
=
0
,
f
′
(
I
)
>
0
,
f
′′
(
I
)
<
0
,
lim
⁡
t
→
∞
f
(
I
)
=
c
<
+
∞
.
(3)
The force of infection at any timet is dominated by β
e
-
d
1
τ
S
(
t
)
f
(
I
(
t
-
τ
)
), where τ is incubation period and 0
<
e
-
d
1
τ
≤
1 represents the survival probability of individuals in the population after time τ [20]. It is also assumed that d
1
≤
d in [
-
τ
,
0
], where d is the death rate and d
1 is the death rate in the time interval [
-
τ
,
0
].
(4)
The parametersa are a convex combination constant, r
=
b
-
d
>
0 is the intrinsic growth rate (b is the birth rate), and K
>
0 is the carrying capacity of the population. The term (
b
-
(
a
r
N
(
t
)
/
K
)
) has a density-dependent per capita birth rate and the term (
d
+
(
(
1
-
a
)
r
N
(
t
)
/
K
)
) has a density-dependent per capita death rate [39].
(5)
For0
<
a
<
1, the birth and death rates are consistent with the limited resources associated with density dependence. The birth rate is density independent when a
=
0 and the death rate is density independent when a
=
1. Thus, the spread of the disease (animals such as rodents, etc.) is assumed to be governed by the following system of logistic equations with time delay.
(6)
The total population is assumed to be large enough to be adequately described by a deterministic model and is divided into compartments based on the disease status [40].
(7)
The successive vaccination rateθ is positive. The positive constant μ
1 is the recovery rate of the infectious individuals from compartment I to S. The parameters β are the effective per capita contact rate constant of infected individuals. The parameters δ are the recovery rate of infected individuals.
(8)
Models are formulated as functional differential and/or integral equations when time delay is included [40]. Ours follows the former with the assumption that the I-equation satisfies a certain integral condition [41].Models with multiple delays are not common, but few authors have in the past considered these-Beretta et al. [42]-to name but a few. Since N
(
t
)
=
S
(
t
)
+
I
(
t
)
+
R
(
t
), thus the governing equation (1) can be rewritten as
(3)
d
N
d
t
=
r
[
1
-
N
K
]
N
,
d
I
d
t
=
β
e
-
d
1
τ
(
N
-
I
-
R
)
f
(
I
(
t
-
τ
)
)
-
[
d
+
(
1
-
a
)
r
N
K
]
I
-
μ
1
I
-
e
-
d
1
ω
δ
I
(
t
-
ω
)
,
d
R
d
t
=
e
-
d
1
ω
δ
I
(
t
-
ω
)
-
[
d
+
(
1
-
a
)
r
N
K
]
R
+
θ
(
N
-
I
-
R
)
.
Let τ
¯
=
max
⁡
(
ω
,
τ
). Then (3) satisfies the following initial conditions
(4)
N
(
0
)
=
S
0
>
0
,
I
(
0
)
=
φ
(
θ
)
≥
0
,
hhhhhhh
∀
θ
∈
[
-
τ
¯
,
0
]
,
R
(
0
)
=
R
0
≥
0
.In this paper, we will consider two different delaysτ, ω which are important parameters on the dynamic behavior. So, the present study is continuation of the previous work τ
=
ω by Naresh et al. [43].Lemma 1.
All solutions of the model system (3) starting in R
+
3 are bounded and eventually enter the compact attracting set
(5)
Φ
=
{
(
S
,
I
,
R
)
∈
R
+
3
:
S
(
t
)
+
I
(
t
)
+
R
(
t
)
=
N
(
t
)
≤
K
}
.Lemma 2.
Let the initial data beN
(
0
)
=
S
0
>
0, I
(
0
)
=
I
0
(
u
)
≥
0, for all u
∈
[
-
τ
¯
,
0
], with I
0
(
0
)
>
0, R
(
0
)
=
R
0
≥
0. Then, the solution (
S
(
t
)
,
I
(
t
)
,
R
(
t
)
) of the model remains positive for all time t
>
0.Lemma 3 (see [44]).
For the characteristic equation in the formp
(
λ
)
+
q
(
λ
)
e
-
r
λ
=
0, where p and q are polynomials with real coefficients and r
>
0 is the delay, suppose(a)
p
(
λ
)
≠
0, R
(
λ
)
>
0;
(b)
|
q
(
i
y
)
|
<
|
p
(
i
y
)
|; 0
≤
y
<
∞;
(c)
lim
⁡
|
λ
|
→
∞
,
R
(
λ
)
≥
0
|
q
(
λ
)
/
p
(
λ
)
|
=
0.
ThenR
(
λ
)
<
0 for every root λ and all r
>
0.


## 3. Equilibrium and Stability Analysis

In this section, we focus on the existence and local stability of equilibria. Let the right-hand side of equalities in model (3) be zero. Then, there are two equilibria; namely,(i)
E
0
=
(
K
,
0
,
p
), p
=
K
θ
/
[
d
+
(
1
-
a
)
r
+
θ
]
), disease-free equilibrium;
(ii)
E
*
=
(
N
*
,
I
*
,
R
*
), endemic equilibrium,where the values ofN
*, I
*, and R
* are given in Section 3.2.


### 3.1. Community Matrix

Firstly, after computing the Jacobian or community matrix of model (3) at point (
N
,
I
,
R
), the characteristic equation is given by(6)
|
r
-
2
r
K
N
-
λ
0
0
β
e
-
d
1
τ
f
(
I
)
-
(
1
-
a
)
r
K
I
m
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
ω
-
n
-
λ
β
e
-
d
1
τ
f
(
I
)
θ
-
(
1
-
a
)
r
K
R
δ
e
-
(
d
1
+
λ
)
ω
-
θ
-
[
d
+
(
1
-
a
)
r
K
N
+
θ
]
-
λ
|
=
0
,where m
=
β
(
N
-
I
-
R
)
f
′
(
I
), n
=
β
e
-
d
1
τ
f
(
I
)
+
d
+
(
(
1
-
a
)
r
/
K
)
N
+
μ
1.Now, we analyze the equilibria stability of system (3). Computing the Jacobian of system (3) evaluated at E
0, one gets the following matrix(7)
J
(
E
0
)
=
(
-
r
-
λ
0
0
0
a
22
-
λ
0
(
d
+
θ
)
θ
d
+
(
1
-
a
)
r
+
θ
δ
e
-
(
d
1
+
λ
)
τ
-
θ
-
[
d
+
(
1
-
a
)
r
+
θ
]
-
λ
)
,where
(8)
a
22
=
β
(
K
-
p
)
f
′
(
0
)
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
ω
-
[
d
+
(
1
-
a
)
r
+
μ
1
]
.
Denote
(9)
A
=
β
(
K
-
p
)
f
′
(
0
)
,
C
=
d
+
(
1
-
a
)
r
+
μ
1
;
then
(10)
a
22
=
A
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
ω
-
C
.
Denote
(11)
h
(
λ
)
=
A
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
ω
-
C
-
λ
.The eigenvalues of the system (3) about the steady state E
0 are λ
1
=
-
r, h
(
λ
)
=
0 and λ
3
=
-
[
d
+
(
1
-
a
)
r
+
θ
]. All the parameters of the model are assumed to be nonnegative. Therefore, λ
1 and λ
3 are negative. Next, we discuss the roots of h
(
λ
)
=
0 in five cases.Case 1.
Forτ
=
ω
≠
0, from the second equation of the system (3), we can get the following.Proposition 4.
Forτ
=
ω
>
0, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(12)
(
A
-
δ
)
e
-
d
1
τ
<
C
.Proof.
From the above analysis,λ
2 satisfies the following characteristic equation:
(13)
g
(
λ
)
=
(
A
-
δ
)
e
-
(
d
1
+
λ
)
τ
-
C
-
λ
=
0
.
(1)
Clearly,λ
=
0 is not a root of g
(
λ
)
=
0.
(2)
From the fact thatg
(
0
)
<
0, g
′
(
λ
)
<
0 for λ
>
0, it is obtained that g
(
λ
)
=
0 has no positive real root.
(3)
It is sufficient to show thatg
(
λ
)
=
0 does not admit a purely imaginary root. In fact, if λ
=
i
v
(
v
>
0
) is a root of (g
(
λ
)
=
0), then by separating the real part, one gets
(14)
(
A
-
δ
)
e
-
d
1
τ
cos
⁡
(
v
τ
)
=
C
.Together with the condition of Proposition 4, we have
(15)
cos
⁡
(
v
τ
)
>
1
.
This is impossible.(4)
It is easy to show thatg
(
λ
)
=
0 has no imaginary root whose real part is positive. Otherwise, there is an imaginary root λ
=
u
+
i
v with u
>
0. Without any loss of generality, we consider v
>
0. Then, we take the real and imaginary parts of g
(
λ
)
=
0; namely,
(16)
(
A
-
δ
)
e
-
(
d
1
+
u
)
τ
cos
⁡
(
v
τ
)
=
C
+
u
.Combined with (9), we have
(17)
C
>
(
A
-
δ
)
e
-
d
1
τ
>
(
A
-
δ
)
e
-
(
d
1
+
u
)
τ
cos
⁡
(
v
τ
)
=
C
+
u
.
This is a contradiction which implies that all eigenvalues roots of g
(
λ
) have negative real parts. Therefore, the disease-free equilibrium of the system (3) is locally asymptotically stable when (9) holds. The proof is completed.Case 2.
Forτ
≠
0
,
ω
=
0, by the same way as in Case 1, one gets the following.Proposition 5.
For allτ
≠
0
,
ω
=
0, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(18)
A
e
-
d
1
τ
-
δ
<
C
.Case 3.
Forω
≠
0
,
τ
=
0, one gets the following.Proposition 6.
Forω
≠
0
,
τ
=
0, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(19)
δ
e
-
d
1
ω
<
|
A
-
C
|
.Proof.
By the fact thath
(
λ
)
=
0 is equivalent to p
(
λ
)
+
q
(
λ
)
e
-
λ
ω
=
0 with p
(
λ
)
=
A
-
C
-
λ, q
(
λ
)
=
-
δ
e
-
d
1
ω.(i)
Supposeλ
=
u
+
i
v
(
u
>
0
). Then, p
(
λ
)
=
A
-
C
-
u
-
i
v
≠
0.
(ii)
By|
q
(
i
v
)
|
=
δ
e
-
d
1
ω, |
p
(
i
v
)
|
=
|
A
-
C
-
i
v
|
=
(
A
-
c
)
2
+
v
2, together with the condition of Proposition 6, we know that |
q
(
i
v
)
|
<
|
p
(
i
v
)
|.
(iii)
Supposeλ
=
u
+
i
v
,
(
u
>
0
). Then,
(20)
lim
⁡
u
2
+
v
2
→
+
∞
|
q
(
λ
)
p
(
λ
)
|
=
δ
e
-
d
1
ω
,
lim
⁡
u
2
+
v
2
→
+
∞
1
(
A
-
C
-
u
)
2
+
v
2
=
0
.Then, using Lemma 3, we have R
(
λ
)
<
0 for all ω.Case 4.
Forω
>
τ
>
0, let ɛ
=
ω
-
τ. Then, ω
=
τ
+
ɛ. For fixed τ,
(21)
h
(
λ
)
=
A
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
(
τ
+
ɛ
)
-
C
-
λ
.
Let λ
=
u
+
i
v
(
u
>
0
). Then, we take the real and imaginary parts of h
(
λ
)
=
0; namely,
(22)
A
e
-
(
d
1
+
u
)
τ
cos
⁡
(
v
τ
)
-
δ
e
-
(
d
1
+
u
)
(
τ
+
ɛ
)
cos
⁡
(
v
(
τ
+
ɛ
)
)
=
C
+
u
,
-
A
e
-
(
d
1
+
u
)
τ
sin
(
v
τ
)
-
δ
e
-
(
d
1
+
u
)
(
τ
+
ɛ
)
sin
(
v
(
τ
+
ɛ
)
)
=
v
.
Sum of squares of the above equalities is
(23)
A
2
e
-
2
(
d
1
+
u
)
τ
+
δ
2
e
-
2
(
d
1
+
u
)
(
τ
+
ɛ
)
-
2
A
δ
e
-
(
d
1
+
u
)
(
2
τ
+
ɛ
)
cos
⁡
(
v
ɛ
)
-
(
C
+
u
)
2
-
v
2
=
0
.
Then, we have
(24)
∂
u
∂
ɛ
|
ɛ
=
0
=
δ
(
u
+
d
1
)
(
A
-
δ
)
(
A
-
δ
)
2
τ
+
(
C
+
u
)
e
2
τ
(
u
+
δ
)
.
Obviously, ∂
u
/
∂
ɛ
|
ɛ
=
0
<
0 when A
-
δ
<
0. Combined with (
A
-
δ
)
e
-
d
1
τ
<
A
e
-
d
1
τ
-
δ
e
-
d
1
ω, we have the following.Proposition 7.
Forω
>
τ
>
0, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(25)
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
<
C
,
A
-
δ
<
0
,
0
<
ω
-
τ
≪
1
.Case 5.
For0
<
τ
<
ω, in the same way as in Case 4, we have the following.Proposition 8.
For0
<
τ
<
ω, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(26)
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
<
C
,
A
-
δ
>
0
,
0
<
τ
-
ω
≪
1
.From what has been discussed above, we get the following.Theorem 9.
The disease-free equilibrium of the system (3) is locally asymptotically stable if one of the following conditions holds.(a)
τ
=
ω
≠
0, (
A
-
δ
)
e
-
d
1
τ
<
C.
(b)
τ
≠
0, ω
=
0, A
e
-
d
1
τ
-
δ
<
C.
(c)
τ
=
0, ω
≠
0, δ
e
-
d
1
ω
<
|
A
-
C
|.
(d)
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
<
C, A
-
δ
<
0, 0
<
ω
-
τ
≪
1.
(e)
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
<
C, A
-
δ
>
0, 0
<
τ
-
ω
≪
1.


### 3.2. Existence of Endemic Equilibrium

Thus, by Theorem9, we may define the basic reproduction number as
(27)
R
0
=
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
C
.This thresholdR
0 defines the average number of secondary infections generated by a typical infectious individual in a completely susceptible population in a steady demographic state.In Theorem9, we have already shown that the system (3) has an infection-free steady state which is locally asymptotically stable under condition R
0
<
1. The disease-free equilibrium is unstable when R
0
>
1, and the system (3) has a nontrivial endemic equilibrium E
*
=
(
N
*
,
I
*
,
R
*
) when R
0
>
1. From (3),
(28)
N
*
=
K
>
0
,
R
*
=
δ
e
-
d
1
ω
-
θ
d
+
(
1
-
a
)
r
+
θ
I
*
+
K
θ
d
+
(
1
-
a
)
r
+
θ
≐
q
I
*
+
p
,
where q
=
(
e
-
d
1
ω
δ
-
θ
)
/
(
d
+
(
1
-
a
)
r
+
θ
). Substituting these values of N
* and R
* in the second equation of (3), we get the following equation for I:
(29)
G
(
I
)
=
β
e
-
d
1
τ
(
K
-
(
1
+
q
)
I
-
p
)
f
(
I
)
-
[
C
+
δ
e
-
d
1
ω
]
I
.
Obviously, I
=
0 is one of the roots of (29) as f
(
0
)
=
0. Therefore, to exclude that root, choose
(30)
H
(
I
)
=
β
e
-
d
1
τ
(
K
-
(
1
+
q
)
I
-
p
)
f
(
I
)
I
-
[
C
+
δ
e
-
d
1
ω
]
.It can easily be seen that the functionH
(
I
) is negative for large positive I; that is,
(31)
H
(
K
)
=
-
β
e
-
d
1
τ
(
K
q
+
p
)
f
(
k
)
K
-
[
C
+
δ
e
-
d
1
ω
]
<
0
.Next, we determine the sign of its derivative(32)
H
′
(
I
)
=
β
e
-
d
1
τ
(
K
-
p
)
f
′
(
I
)
I
-
f
(
I
)
I
2
-
β
e
-
d
1
τ
(
1
+
q
)
f
′
(
I
)
.
It can easily be seen that K
>
p. In addition, from the properties of the function f
(
I
), in particular from f
(
0
)
=
0 and f
′′
(
0
)
<
0, it follows that f
(
I
)
-
f
′
(
I
)
I
>
0, and consequently H
′
(
I
)
<
0 for all I
>
0. Therefore, for a positive root of H
(
I
)
=
0 to exist, H
(
I
) has to satisfy H
(
0
)
>
0; that is,
(33)
H
(
0
)
=
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
-
C
=
(
A
e
-
d
1
τ
-
e
-
d
1
ω
δ
C
-
1
)
C
=
(
R
0
-
1
)
C
.
Hence, one needs the requirement that R
0
>
1 to ensure the existence of the endemic equilibrium. From the above analysis, we have the following theorem.Theorem 10.
The system (3) has a nontrivial endemic equilibrium E
*
=
(
N
*
,
I
*
,
R
*
) when R
0
>
1.


### 3.3. Local Stability of the Endemic Equilibrium

In this section, we analyze the local stability of the endemic equilibriumE
* for τ
=
ω. Its characteristic equation is given by(34)
|
-
r
-
λ
0
0
β
e
-
d
1
τ
f
(
I
*
)
-
(
1
-
a
)
r
I
*
K
m
e
-
(
d
1
+
λ
)
τ
-
n
-
λ
β
e
-
d
1
τ
f
(
I
*
)
θ
-
(
1
-
a
)
r
R
*
K
δ
e
-
(
d
1
+
λ
)
τ
-
θ
-
[
d
+
(
1
-
a
)
r
+
θ
]
-
λ
|
=
0
,where m
=
β
(
K
-
I
*
-
R
*
)
f
′
(
I
*
)
-
δ, n
=
β
e
-
d
1
τ
f
(
I
*
)
+
d
+
(
1
-
a
)
r
+
μ
1.The Jacobin matrix leads to the characteristic equation(35)
(
λ
+
r
)
[
λ
2
+
m
1
λ
+
m
0
+
(
n
1
λ
+
n
0
)
e
-
λ
τ
]
=
0
,
where
(36)
m
1
=
β
e
-
d
1
τ
f
(
I
*
)
+
2
d
+
2
(
1
-
a
)
r
+
θ
+
μ
1
>
0
,
m
0
=
[
β
e
-
d
1
τ
f
(
I
*
)
+
d
+
(
1
-
a
)
r
+
μ
1
]
×
(
d
+
(
1
-
a
)
r
+
θ
)
+
(
d
+
(
1
-
a
)
r
+
μ
1
)
θ
>
0
,
n
1
=
-
[
β
f
′
(
I
*
)
(
K
-
I
*
-
R
*
)
-
δ
]
e
-
d
1
τ
,
n
0
=
[
β
f
′
(
I
*
)
(
K
-
I
*
-
R
*
)
-
δ
]
[
d
+
(
1
-
a
)
r
+
θ
]
×
e
-
d
1
τ
+
δ
f
(
I
*
)
e
-
2
d
1
τ
.
Since all the model parameters are assumed to be nonnegative, it follows that one eigenvalue is negative; that is, λ
1
=
-
r. Thus, the stability of E
* depends on the roots of the quasi-polynomial
(37)
λ
2
+
m
1
λ
+
m
0
+
(
n
1
λ
+
n
0
)
e
-
λ
τ
=
0
.
We note that m
1
>
0 and m
0
>
0, whereas n
1 and n
0 may be positive or negative. For τ
=
0, we state the following results that follow directly from (39). The endemic steady state is locally asymptotically stable if the following conditions hold:
(
H
1
)
β
f
(
I
*
)
+
2
d
+
2
(
1
-
a
)
r
+
θ
+
μ
1
+
δ
β
f
′
(
I
*
)
>
β
(
K
-
I
*
-
R
*
)
f
′
(
I
*
)
,
(
H
2
)
[
β
f
(
I
*
)
+
d
+
(
1
-
a
)
r
+
μ
1
+
δ
-
β
(
K
-
I
*
-
R
*
)
f
′
(
I
*
)
]
×
[
d
+
(
1
-
a
)
r
+
θ
]
+
f
(
I
*
)
δ
>
0
.
The main purpose of this paper is to study the stability behavior of E
* in the case τ
≠
0. Obviously, i
η
(
η
>
0
) is the root of (29) if and only if η satisfies
(38)
-
η
2
+
m
1
i
η
+
m
0
=
-
(
n
1
i
η
+
n
0
)
(
cos
⁡
η
τ
-
i
sin
η
τ
)
.
Separating the real and imaginary parts, we have
(39)
-
η
2
+
m
0
=
-
n
0
cos
⁡
η
τ
-
n
1
η
sin
η
τ
,
(40)
m
1
η
=
-
n
1
η
cos
⁡
η
τ
+
n
0
sin
η
τ
.
Eliminating τ by squaring and adding (39) and (40), we obtain a polynomial in η as
(41)
η
4
+
(
m
1
2
-
n
1
2
-
2
m
0
)
η
2
+
m
0
2
-
n
0
2
=
0
.
Suppose that the conditions
(
H
3
)
m
1
2
>
n
1
2
+
2
m
0
,
m
0
2
>
n
0
2
hold for all τ
≥
0. Then, the infected steady state of the system (3) is locally asymptotically stable.Theorem 11.
Forτ
=
ω, if R
0
>
1, then the endemic equilibrium of the system (3) is locally asymptotically stable, when conditions (
H
1
)–(
H
3
) hold.Corollary 12.
Forτ
=
ω, if μ
1
<
μ
1
* or θ
<
θ
*, then the endemic equilibrium of the system (3) is locally asymptotically stable, when conditions (
H
1
)–(
H
3
) hold.


## 3.1. Community Matrix

Firstly, after computing the Jacobian or community matrix of model (3) at point (
N
,
I
,
R
), the characteristic equation is given by(6)
|
r
-
2
r
K
N
-
λ
0
0
β
e
-
d
1
τ
f
(
I
)
-
(
1
-
a
)
r
K
I
m
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
ω
-
n
-
λ
β
e
-
d
1
τ
f
(
I
)
θ
-
(
1
-
a
)
r
K
R
δ
e
-
(
d
1
+
λ
)
ω
-
θ
-
[
d
+
(
1
-
a
)
r
K
N
+
θ
]
-
λ
|
=
0
,where m
=
β
(
N
-
I
-
R
)
f
′
(
I
), n
=
β
e
-
d
1
τ
f
(
I
)
+
d
+
(
(
1
-
a
)
r
/
K
)
N
+
μ
1.Now, we analyze the equilibria stability of system (3). Computing the Jacobian of system (3) evaluated at E
0, one gets the following matrix(7)
J
(
E
0
)
=
(
-
r
-
λ
0
0
0
a
22
-
λ
0
(
d
+
θ
)
θ
d
+
(
1
-
a
)
r
+
θ
δ
e
-
(
d
1
+
λ
)
τ
-
θ
-
[
d
+
(
1
-
a
)
r
+
θ
]
-
λ
)
,where
(8)
a
22
=
β
(
K
-
p
)
f
′
(
0
)
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
ω
-
[
d
+
(
1
-
a
)
r
+
μ
1
]
.
Denote
(9)
A
=
β
(
K
-
p
)
f
′
(
0
)
,
C
=
d
+
(
1
-
a
)
r
+
μ
1
;
then
(10)
a
22
=
A
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
ω
-
C
.
Denote
(11)
h
(
λ
)
=
A
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
ω
-
C
-
λ
.The eigenvalues of the system (3) about the steady state E
0 are λ
1
=
-
r, h
(
λ
)
=
0 and λ
3
=
-
[
d
+
(
1
-
a
)
r
+
θ
]. All the parameters of the model are assumed to be nonnegative. Therefore, λ
1 and λ
3 are negative. Next, we discuss the roots of h
(
λ
)
=
0 in five cases.Case 1.
Forτ
=
ω
≠
0, from the second equation of the system (3), we can get the following.Proposition 4.
Forτ
=
ω
>
0, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(12)
(
A
-
δ
)
e
-
d
1
τ
<
C
.Proof.
From the above analysis,λ
2 satisfies the following characteristic equation:
(13)
g
(
λ
)
=
(
A
-
δ
)
e
-
(
d
1
+
λ
)
τ
-
C
-
λ
=
0
.
(1)
Clearly,λ
=
0 is not a root of g
(
λ
)
=
0.
(2)
From the fact thatg
(
0
)
<
0, g
′
(
λ
)
<
0 for λ
>
0, it is obtained that g
(
λ
)
=
0 has no positive real root.
(3)
It is sufficient to show thatg
(
λ
)
=
0 does not admit a purely imaginary root. In fact, if λ
=
i
v
(
v
>
0
) is a root of (g
(
λ
)
=
0), then by separating the real part, one gets
(14)
(
A
-
δ
)
e
-
d
1
τ
cos
⁡
(
v
τ
)
=
C
.Together with the condition of Proposition 4, we have
(15)
cos
⁡
(
v
τ
)
>
1
.
This is impossible.(4)
It is easy to show thatg
(
λ
)
=
0 has no imaginary root whose real part is positive. Otherwise, there is an imaginary root λ
=
u
+
i
v with u
>
0. Without any loss of generality, we consider v
>
0. Then, we take the real and imaginary parts of g
(
λ
)
=
0; namely,
(16)
(
A
-
δ
)
e
-
(
d
1
+
u
)
τ
cos
⁡
(
v
τ
)
=
C
+
u
.Combined with (9), we have
(17)
C
>
(
A
-
δ
)
e
-
d
1
τ
>
(
A
-
δ
)
e
-
(
d
1
+
u
)
τ
cos
⁡
(
v
τ
)
=
C
+
u
.
This is a contradiction which implies that all eigenvalues roots of g
(
λ
) have negative real parts. Therefore, the disease-free equilibrium of the system (3) is locally asymptotically stable when (9) holds. The proof is completed.Case 2.
Forτ
≠
0
,
ω
=
0, by the same way as in Case 1, one gets the following.Proposition 5.
For allτ
≠
0
,
ω
=
0, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(18)
A
e
-
d
1
τ
-
δ
<
C
.Case 3.
Forω
≠
0
,
τ
=
0, one gets the following.Proposition 6.
Forω
≠
0
,
τ
=
0, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(19)
δ
e
-
d
1
ω
<
|
A
-
C
|
.Proof.
By the fact thath
(
λ
)
=
0 is equivalent to p
(
λ
)
+
q
(
λ
)
e
-
λ
ω
=
0 with p
(
λ
)
=
A
-
C
-
λ, q
(
λ
)
=
-
δ
e
-
d
1
ω.(i)
Supposeλ
=
u
+
i
v
(
u
>
0
). Then, p
(
λ
)
=
A
-
C
-
u
-
i
v
≠
0.
(ii)
By|
q
(
i
v
)
|
=
δ
e
-
d
1
ω, |
p
(
i
v
)
|
=
|
A
-
C
-
i
v
|
=
(
A
-
c
)
2
+
v
2, together with the condition of Proposition 6, we know that |
q
(
i
v
)
|
<
|
p
(
i
v
)
|.
(iii)
Supposeλ
=
u
+
i
v
,
(
u
>
0
). Then,
(20)
lim
⁡
u
2
+
v
2
→
+
∞
|
q
(
λ
)
p
(
λ
)
|
=
δ
e
-
d
1
ω
,
lim
⁡
u
2
+
v
2
→
+
∞
1
(
A
-
C
-
u
)
2
+
v
2
=
0
.Then, using Lemma 3, we have R
(
λ
)
<
0 for all ω.Case 4.
Forω
>
τ
>
0, let ɛ
=
ω
-
τ. Then, ω
=
τ
+
ɛ. For fixed τ,
(21)
h
(
λ
)
=
A
e
-
(
d
1
+
λ
)
τ
-
δ
e
-
(
d
1
+
λ
)
(
τ
+
ɛ
)
-
C
-
λ
.
Let λ
=
u
+
i
v
(
u
>
0
). Then, we take the real and imaginary parts of h
(
λ
)
=
0; namely,
(22)
A
e
-
(
d
1
+
u
)
τ
cos
⁡
(
v
τ
)
-
δ
e
-
(
d
1
+
u
)
(
τ
+
ɛ
)
cos
⁡
(
v
(
τ
+
ɛ
)
)
=
C
+
u
,
-
A
e
-
(
d
1
+
u
)
τ
sin
(
v
τ
)
-
δ
e
-
(
d
1
+
u
)
(
τ
+
ɛ
)
sin
(
v
(
τ
+
ɛ
)
)
=
v
.
Sum of squares of the above equalities is
(23)
A
2
e
-
2
(
d
1
+
u
)
τ
+
δ
2
e
-
2
(
d
1
+
u
)
(
τ
+
ɛ
)
-
2
A
δ
e
-
(
d
1
+
u
)
(
2
τ
+
ɛ
)
cos
⁡
(
v
ɛ
)
-
(
C
+
u
)
2
-
v
2
=
0
.
Then, we have
(24)
∂
u
∂
ɛ
|
ɛ
=
0
=
δ
(
u
+
d
1
)
(
A
-
δ
)
(
A
-
δ
)
2
τ
+
(
C
+
u
)
e
2
τ
(
u
+
δ
)
.
Obviously, ∂
u
/
∂
ɛ
|
ɛ
=
0
<
0 when A
-
δ
<
0. Combined with (
A
-
δ
)
e
-
d
1
τ
<
A
e
-
d
1
τ
-
δ
e
-
d
1
ω, we have the following.Proposition 7.
Forω
>
τ
>
0, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(25)
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
<
C
,
A
-
δ
<
0
,
0
<
ω
-
τ
≪
1
.Case 5.
For0
<
τ
<
ω, in the same way as in Case 4, we have the following.Proposition 8.
For0
<
τ
<
ω, R
(
λ
)
<
0 for every root λ of h
(
λ
)
=
0 when
(26)
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
<
C
,
A
-
δ
>
0
,
0
<
τ
-
ω
≪
1
.From what has been discussed above, we get the following.Theorem 9.
The disease-free equilibrium of the system (3) is locally asymptotically stable if one of the following conditions holds.(a)
τ
=
ω
≠
0, (
A
-
δ
)
e
-
d
1
τ
<
C.
(b)
τ
≠
0, ω
=
0, A
e
-
d
1
τ
-
δ
<
C.
(c)
τ
=
0, ω
≠
0, δ
e
-
d
1
ω
<
|
A
-
C
|.
(d)
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
<
C, A
-
δ
<
0, 0
<
ω
-
τ
≪
1.
(e)
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
<
C, A
-
δ
>
0, 0
<
τ
-
ω
≪
1.


## 3.2. Existence of Endemic Equilibrium

Thus, by Theorem9, we may define the basic reproduction number as
(27)
R
0
=
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
C
.This thresholdR
0 defines the average number of secondary infections generated by a typical infectious individual in a completely susceptible population in a steady demographic state.In Theorem9, we have already shown that the system (3) has an infection-free steady state which is locally asymptotically stable under condition R
0
<
1. The disease-free equilibrium is unstable when R
0
>
1, and the system (3) has a nontrivial endemic equilibrium E
*
=
(
N
*
,
I
*
,
R
*
) when R
0
>
1. From (3),
(28)
N
*
=
K
>
0
,
R
*
=
δ
e
-
d
1
ω
-
θ
d
+
(
1
-
a
)
r
+
θ
I
*
+
K
θ
d
+
(
1
-
a
)
r
+
θ
≐
q
I
*
+
p
,
where q
=
(
e
-
d
1
ω
δ
-
θ
)
/
(
d
+
(
1
-
a
)
r
+
θ
). Substituting these values of N
* and R
* in the second equation of (3), we get the following equation for I:
(29)
G
(
I
)
=
β
e
-
d
1
τ
(
K
-
(
1
+
q
)
I
-
p
)
f
(
I
)
-
[
C
+
δ
e
-
d
1
ω
]
I
.
Obviously, I
=
0 is one of the roots of (29) as f
(
0
)
=
0. Therefore, to exclude that root, choose
(30)
H
(
I
)
=
β
e
-
d
1
τ
(
K
-
(
1
+
q
)
I
-
p
)
f
(
I
)
I
-
[
C
+
δ
e
-
d
1
ω
]
.It can easily be seen that the functionH
(
I
) is negative for large positive I; that is,
(31)
H
(
K
)
=
-
β
e
-
d
1
τ
(
K
q
+
p
)
f
(
k
)
K
-
[
C
+
δ
e
-
d
1
ω
]
<
0
.Next, we determine the sign of its derivative(32)
H
′
(
I
)
=
β
e
-
d
1
τ
(
K
-
p
)
f
′
(
I
)
I
-
f
(
I
)
I
2
-
β
e
-
d
1
τ
(
1
+
q
)
f
′
(
I
)
.
It can easily be seen that K
>
p. In addition, from the properties of the function f
(
I
), in particular from f
(
0
)
=
0 and f
′′
(
0
)
<
0, it follows that f
(
I
)
-
f
′
(
I
)
I
>
0, and consequently H
′
(
I
)
<
0 for all I
>
0. Therefore, for a positive root of H
(
I
)
=
0 to exist, H
(
I
) has to satisfy H
(
0
)
>
0; that is,
(33)
H
(
0
)
=
A
e
-
d
1
τ
-
δ
e
-
d
1
ω
-
C
=
(
A
e
-
d
1
τ
-
e
-
d
1
ω
δ
C
-
1
)
C
=
(
R
0
-
1
)
C
.
Hence, one needs the requirement that R
0
>
1 to ensure the existence of the endemic equilibrium. From the above analysis, we have the following theorem.Theorem 10.
The system (3) has a nontrivial endemic equilibrium E
*
=
(
N
*
,
I
*
,
R
*
) when R
0
>
1.


## 3.3. Local Stability of the Endemic Equilibrium

In this section, we analyze the local stability of the endemic equilibriumE
* for τ
=
ω. Its characteristic equation is given by(34)
|
-
r
-
λ
0
0
β
e
-
d
1
τ
f
(
I
*
)
-
(
1
-
a
)
r
I
*
K
m
e
-
(
d
1
+
λ
)
τ
-
n
-
λ
β
e
-
d
1
τ
f
(
I
*
)
θ
-
(
1
-
a
)
r
R
*
K
δ
e
-
(
d
1
+
λ
)
τ
-
θ
-
[
d
+
(
1
-
a
)
r
+
θ
]
-
λ
|
=
0
,where m
=
β
(
K
-
I
*
-
R
*
)
f
′
(
I
*
)
-
δ, n
=
β
e
-
d
1
τ
f
(
I
*
)
+
d
+
(
1
-
a
)
r
+
μ
1.The Jacobin matrix leads to the characteristic equation(35)
(
λ
+
r
)
[
λ
2
+
m
1
λ
+
m
0
+
(
n
1
λ
+
n
0
)
e
-
λ
τ
]
=
0
,
where
(36)
m
1
=
β
e
-
d
1
τ
f
(
I
*
)
+
2
d
+
2
(
1
-
a
)
r
+
θ
+
μ
1
>
0
,
m
0
=
[
β
e
-
d
1
τ
f
(
I
*
)
+
d
+
(
1
-
a
)
r
+
μ
1
]
×
(
d
+
(
1
-
a
)
r
+
θ
)
+
(
d
+
(
1
-
a
)
r
+
μ
1
)
θ
>
0
,
n
1
=
-
[
β
f
′
(
I
*
)
(
K
-
I
*
-
R
*
)
-
δ
]
e
-
d
1
τ
,
n
0
=
[
β
f
′
(
I
*
)
(
K
-
I
*
-
R
*
)
-
δ
]
[
d
+
(
1
-
a
)
r
+
θ
]
×
e
-
d
1
τ
+
δ
f
(
I
*
)
e
-
2
d
1
τ
.
Since all the model parameters are assumed to be nonnegative, it follows that one eigenvalue is negative; that is, λ
1
=
-
r. Thus, the stability of E
* depends on the roots of the quasi-polynomial
(37)
λ
2
+
m
1
λ
+
m
0
+
(
n
1
λ
+
n
0
)
e
-
λ
τ
=
0
.
We note that m
1
>
0 and m
0
>
0, whereas n
1 and n
0 may be positive or negative. For τ
=
0, we state the following results that follow directly from (39). The endemic steady state is locally asymptotically stable if the following conditions hold:
(
H
1
)
β
f
(
I
*
)
+
2
d
+
2
(
1
-
a
)
r
+
θ
+
μ
1
+
δ
β
f
′
(
I
*
)
>
β
(
K
-
I
*
-
R
*
)
f
′
(
I
*
)
,
(
H
2
)
[
β
f
(
I
*
)
+
d
+
(
1
-
a
)
r
+
μ
1
+
δ
-
β
(
K
-
I
*
-
R
*
)
f
′
(
I
*
)
]
×
[
d
+
(
1
-
a
)
r
+
θ
]
+
f
(
I
*
)
δ
>
0
.
The main purpose of this paper is to study the stability behavior of E
* in the case τ
≠
0. Obviously, i
η
(
η
>
0
) is the root of (29) if and only if η satisfies
(38)
-
η
2
+
m
1
i
η
+
m
0
=
-
(
n
1
i
η
+
n
0
)
(
cos
⁡
η
τ
-
i
sin
η
τ
)
.
Separating the real and imaginary parts, we have
(39)
-
η
2
+
m
0
=
-
n
0
cos
⁡
η
τ
-
n
1
η
sin
η
τ
,
(40)
m
1
η
=
-
n
1
η
cos
⁡
η
τ
+
n
0
sin
η
τ
.
Eliminating τ by squaring and adding (39) and (40), we obtain a polynomial in η as
(41)
η
4
+
(
m
1
2
-
n
1
2
-
2
m
0
)
η
2
+
m
0
2
-
n
0
2
=
0
.
Suppose that the conditions
(
H
3
)
m
1
2
>
n
1
2
+
2
m
0
,
m
0
2
>
n
0
2
hold for all τ
≥
0. Then, the infected steady state of the system (3) is locally asymptotically stable.Theorem 11.
Forτ
=
ω, if R
0
>
1, then the endemic equilibrium of the system (3) is locally asymptotically stable, when conditions (
H
1
)–(
H
3
) hold.Corollary 12.
Forτ
=
ω, if μ
1
<
μ
1
* or θ
<
θ
*, then the endemic equilibrium of the system (3) is locally asymptotically stable, when conditions (
H
1
)–(
H
3
) hold.


## 4. Permanence

In this section, we investigate a permanence result [5]. The following is our main result of this paper. We will give the following result by using some techniques given in [8, 11]. The proof of the permanence with nonlinear incidence is a daunting task. Consequently, for simplicity and mathematical convenience, let us choose a linear incidence rate f
(
I
)
=
I. The result holds with the nonlinear incidence, as shown numerically, but the algebraic proof is long and tedious, and the conditions to impose on some of the parameters may be very restrictive. Now, let us firstly give the following theorem.Theorem 13.
IfR
0
>
1 holds, then the system (1) with τ
=
ω is permanent; that is, there are positive constants c
i
(
i
=
1,2
,
3
) such that
(42)
c
1
<
lim
⁡
t
→
∞
inf
⁡
S
(
t
)
≤
lim
⁡
t
→
∞
sup
⁡
S
(
t
)
≤
K
,
c
2
<
lim
⁡
t
→
∞
inf
⁡
I
(
t
)
≤
lim
⁡
t
→
∞
sup
⁡
I
(
t
)
≤
K
,
c
3
<
lim
⁡
t
→
∞
inf
⁡
R
(
t
)
≤
lim
⁡
t
→
∞
sup
⁡
R
(
t
)
≤
K
hold for any solution of (1) with (
ϕ
1
(
θ
)
,
ϕ
2
(
θ
)
,
ϕ
3
(
θ
)
) in the interior of Φ for all θ
∈
[
-
τ
,
0
]. In fact, c
i
(
i
=
1,2
,
3
) can be chosen explicitly as
(43)
c
1
=
(
b
-
a
r
)
K
β
K
e
-
d
1
τ
+
d
+
(
1
-
a
)
r
+
θ
,
c
2
=
I
*
e
-
(
d
+
(
1
-
a
)
r
+
μ
1
+
e
-
d
1
τ
δ
)
,
c
3
=
δ
e
-
d
1
τ
c
2
+
θ
c
1
d
+
(
1
-
a
)
r
.Proof.
Note that0
<
N
(
t
)
<
K for all t
≥
0 and that lim
⁡
t
→
∞
N
(
t
)
=
K. It is easy to see that lim
⁡
t
→
∞
inf
⁡
S
(
t
)
≥
c
1. In fact, let ϵ
<
K be arbitrary. Choose T
1
>
τ so large that N
(
t
)
>
K
-
ϵ for t
>
T
1. We have the following inequality:
(44)
S
˙
(
t
)
>
-
[
β
K
e
-
d
1
τ
+
d
+
(
1
-
a
)
r
+
θ
]
S
(
t
)
+
(
b
-
a
r
)
(
K
-
ϵ
)
,
for all t
≥
T
1, which implies that
(45)
liminf
⁡
t
→
∞
S
(
t
)
≥
(
b
-
a
r
)
(
K
-
ϵ
)
β
K
e
-
d
1
τ
+
d
+
(
1
-
a
)
r
+
θ
.
Note that ϵ may be arbitrarily small so that lim
⁡
t
→
∞
inf
⁡
S
(
t
)
≥
c
1.
Next, we will showlim
⁡
t
→
∞
inf
⁡
I
(
t
)
≥
c
2. For any ξ
:
0
<
ξ
<
1, we see the inequality S
*
<
[
(
b
-
a
r
)
K
+
μ
1
I
*
]
/
(
β
e
-
d
1
τ
ξ
I
*
+
d
+
(
1
-
a
)
r
+
θ
). There exist sufficiently large ρ
≥
1 and sufficiently small ϵ such that S
*
<
{
[
(
b
-
a
r
)
(
K
-
ϵ
)
+
μ
1
I
*
]/(
β
e
-
d
1
τ
ξ
I
*
+
d
+
(
1
-
a
)
r
+
θ
)
}
(
1
-
e
-
(
β
e
-
d
1
τ
ξ
I
*
+
d
+
(
1
-
a
)
r
+
θ
)
ρ
τ
)
≡
S
Δ. We show that I
(
t
0
)
>
q
I
* for some t
0
≥
ρ
τ. In fact, if not, it follows from the first equation of (1) that, for all t
≥
ρ
τ
+
τ
≥
T
1
+
τ,
(46)
S
˙
(
t
)
≥
(
b
-
a
r
)
(
K
-
ϵ
)
+
μ
1
I
*
-
[
β
e
-
d
1
τ
ξ
I
*
+
d
+
(
1
-
a
)
r
+
θ
]
S
(
t
)
.
Hence, for t
≥
ρ
τ
+
τ,
(47)
S
(
t
)
≥
e
-
(
β
e
-
d
1
τ
ξ
I
*
+
d
+
(
1
-
a
)
r
+
θ
)
(
t
-
ρ
τ
-
τ
)
×
[
S
(
ρ
τ
+
τ
)
+
(
b
-
a
r
)
(
K
-
ϵ
)
+
μ
1
I
*
]
×
∫
ρ
τ
+
τ
t
‍
e
-
(
β
e
-
d
1
τ
ξ
I
*
+
d
+
(
1
-
a
)
r
+
θ
)
(
t
-
ρ
τ
-
τ
)
d
θ
>
[
(
b
-
a
r
)
(
K
-
ϵ
)
+
μ
1
I
*
]
β
e
-
d
1
τ
ξ
I
*
+
d
+
(
1
-
a
)
r
+
θ
×
(
1
-
e
-
(
β
e
-
d
1
τ
ξ
I
*
+
d
+
(
1
-
a
)
r
+
θ
)
(
t
-
ρ
τ
-
τ
)
)
,
which gives us, for t
≥
2
ρ
τ
+
τ,
(48)
S
(
t
)
>
S
Δ
>
S
*
.
For t
≥
0, we define a positive differentiable function V
(
t
) as follows:
(49)
V
(
t
)
=
I
(
t
)
+
[
β
(
K
-
p
)
-
δ
]
e
-
d
1
τ
R
0
∫
t
-
τ
t
I
(
s
)
d
s
.
We obtain the inequality, for t
≥
2
ρ
τ
+
τ,
(50)
V
˙
(
t
)
=
[
β
e
-
d
1
τ
(
S
(
t
)
-
S
*
)
I
(
t
-
τ
)
+
(
1
-
a
)
r
]
×
(
1
-
N
(
t
)
K
)
I
(
t
)
>
β
e
-
d
1
τ
(
S
(
t
)
-
S
*
)
I
(
t
-
τ
)
>
β
e
-
d
1
τ
(
S
Δ
-
S
*
)
I
(
t
-
τ
)
.
Leti
_
=
min
⁡
θ
∈
[
-
τ
,
0
]
I
(
2
ρ
τ
+
2
τ
+
θ
). Now, let us show that I
(
t
)
≥
i
_ for all t
≥
2
ρ
τ
+
τ. In fact, if there exists T
2
≥
0 such that I
(
t
)
≥
i
_ for 2
ρ
τ
+
τ
≤
t
≤
2
ρ
τ
+
2
τ
+
T
2, I
(
2
ρ
τ
+
2
τ
+
T
2
)
=
i
_ and I
˙
(
2
ρ
τ
+
2
τ
+
T
2
)
≤
0. Direct calculation using the second equation of (1) and (48) gives
(51)
I
˙
(
2
ρ
τ
+
2
τ
+
T
2
)
>
[
(
d
+
(
1
-
a
)
r
+
μ
1
+
e
-
d
1
τ
δ
)
β
e
-
d
1
τ
(
S
(
2
ρ
τ
+
2
τ
+
T
2
)
h
h
h
h
h
h
h
-
(
d
+
(
1
-
a
)
r
+
μ
1
+
e
-
d
1
τ
δ
)
]
i
_
>
(
d
+
(
1
-
a
)
r
+
μ
1
+
e
-
d
1
τ
δ
)
[
S
Δ
S
*
-
1
]
i
_
>
0
.
This contradicts the definition of T
2. Thus, we have shown that I
(
t
)
≥
i
_ for all t
≥
2
ρ
τ
+
τ. Hence, for all t
≥
2
ρ
τ
+
2
τ,
(52)
V
˙
(
t
)
>
β
e
-
d
1
τ
(
S
Δ
-
S
*
)
i
_
,
which implies that V
(
t
)
→
+
∞ as t
→
+
∞. This contradicts the boundedness of V
(
t
). Consequently, I
(
t
0
)
>
ξ
I
* for some t
0
≥
ρ
τ.
In the rest, we now need to consider two cases:(i)
I
(
t
)
≥
ξ
I
* for all large t;
(ii)
I
(
t
) oscillates about ξ
I
* for all large t.
We now need to show thatI
(
t
)
≥
ξ
c
2 for large t. Obviously, it suffices to show that it holds only for case (ii). We suppose that for any large T there exists t
1, t
2
>
T such that I
(
t
1
)
=
I
(
t
2
)
=
ξ
I
* and I
(
t
)
<
ξ
I
* for t
1
<
t
<
t
2. If t
2
-
t
1
≤
τ, the second equation (1) gives us I
˙
(
t
)
>
-
(
d
+
(
1
-
a
)
r
+
μ
1
)
I
(
t
), which implies that I
(
t
)
>
I
(
t
1
)
e
-
(
d
+
(
1
-
a
)
r
+
μ
1
)
(
t
-
t
1
) on (
t
1
,
t
2
). Thus, I
(
t
)
>
ξ
c
2. On the other hand, if t
2
-
t
1
>
τ, applying the same manner gives I
(
t
)
≥
ξ
c
2 on [
t
1
,
t
1
+
τ
], and hence the remaining work is to show I
(
t
)
≥
ξ
c
2 on [
t
1
+
τ
,
t
2
]. In fact, assuming that there exists T
3
>
0 such that I
(
t
)
≥
ξ
c
2 on [
t
1
,
t
1
+
τ
+
T
3
], I
(
t
1
+
τ
+
T
3
)
=
ξ
c
2, and I
˙
(
t
1
+
τ
+
T
3
)
≤
0, it follows from (1) that
(53)
I
˙
(
t
1
+
τ
+
T
3
)
≥
[
(
d
+
(
1
-
a
)
r
+
μ
1
+
e
-
d
1
τ
δ
)
β
e
-
d
1
τ
S
(
t
1
+
τ
+
T
3
)
h
h
h
-
(
d
+
(
1
-
a
)
r
+
μ
1
+
e
-
d
1
τ
δ
)
]
ξ
c
2
>
(
d
+
(
1
-
a
)
r
+
μ
1
+
e
-
d
1
τ
δ
)
[
S
Δ
S
*
-
1
]
ξ
c
2
>
0
.
This contradicts the definition of T
3. Hence, I
(
t
)
≥
ξ
c
2 on [
t
1
,
t
2
]. Consequently, I
(
t
)
≥
ξ
c
2 for large t in the case (ii). Therefore, lim
⁡
t
→
∞
inf
⁡
⁡
I
(
t
)
≥
ξ
c
2. Note that q may be so close to 1 that lim
⁡
t
→
∞
inf
⁡
I
(
t
)
≥
c
2.
Finally, let us show thatlim
⁡
t
→
∞
inf
⁡
R
(
t
)
≥
(
δ
e
-
d
1
τ
c
2
+
θ
c
1
)
/
(
d
+
(
1
-
a
)
r
). The third equation gives us
(54)
R
˙
(
t
)
≥
[
δ
e
-
d
1
τ
I
+
θ
S
-
[
d
+
(
1
-
a
)
r
]
R
≥
[
δ
e
-
d
1
τ
ξ
c
2
+
θ
ξ
c
1
-
[
d
+
(
1
-
a
)
r
]
R
for large t. Hence, lim
⁡
t
→
∞
inf
⁡
R
(
t
)
≥
(
δ
e
-
d
1
τ
ξ
c
2
+
θ
ξ
c
1
)
/
(
d
+
(
1
-
a
)
r
). In a similar manner, we could show lim
⁡
t
→
∞
inf
⁡
R
(
t
)
≥
c
3. This proves the theorem.Corollary 14.
Ifμ
1
<
μ
1
* and θ
<
θ
*, then the system (1) with τ
=
ω is permanent.


## 5. Numerical Analysis

Since it is important to visualize the dynamical behavior of the model, the model system (3) is integrated numerically with the help of MATLAB 7.0 using the following set of parameters.(1) Letr
=
0.5, k
=
8, d
=
0.04, d
1
=
0.04, β
=
1, a
=
0.3, δ
=
0.2, μ
1
=
0.8, θ
=
0.01, and τ
=
5. It is easy to compute that E
0
=
(
8,0
,
0.2
) and R
0
=
0.94
<
1. In Figure 1, the infective population and recovered population, respectively, are plotted against the total population. We see from the figure that for any initial start the solution curves tend to the equilibrium E
0. Hence, we infer that the system (3) may be stable about the disease-free equilibrium point E
0, which satisfies Theorem 9.Figure 1
The disease-free equilibriumE
0 is locally asymptotically stable. Variation of infective population I
(
t
) and recovered population R
(
t
) with total population N
(
t
).(2) Letr
=
0.5, k
=
8, d
=
0.04, d
1
=
0.04, β
=
1, a
=
0.3, δ
=
0.2, μ
1
=
0.2, θ
=
0.02, and τ
=
5. We get E
*
=
(
8,2.23,1.11
) and this set of parameter values satisfies the local asymptotic stability conditions of E
*. It is easy to verify that R
0
=
1.83
>
1 and all other conditions of Theorem 11 are satisfied. So, we can obtain from Figure 2 that the system (3) is stable at the endemic equilibrium point E
*.Figure 2
The endemic equilibriumE
* is locally asymptotically stable. Variation of infective population I
(
t
) and recovered population R
(
t
) with total population N
(
t
).(3) The results of numerical simulation are displayed graphically in Figures3 and 4. In Figure 3, the variation of the infective population and recovered population is shown with time for different values of the removal rate constant from groups I to S, μ
1. It is found that both the infective population and the recovered population decrease as μ
1 increases. Figure 4 depicts the variation of infective population and recovered population, respectively, with time for the different successive vaccination rate, θ. As θ increases, the infective population decreases whereas the recovered population increases.Figure 3
Variation of infective populationI
(
t
) and recovered population R
(
t
) with time for different values of μ
1.Figure 4
Variation of infective populationI
(
t
) and recovered population R
(
t
) with time for different values of θ.


## 6. Discussion

In this paper, we will consider two different delays which are important parameters on the dynamic behavior. So, the present study is continuation of the previous work by [43]. Furthermore, from biological epidemic point of view, we investigate successive vaccination and difference in immunity in our system. From mathematical point of view, we study the stability of disease-free equilibrium and the existence of endemic equilibrium for different delay and consider the permanence of the system in the new paper.In Theorems9, 10, 11, and 13 corresponding to their corollaries, when the effect of the successive vaccination rate and the transfer rate from the infectious group to the susceptible group after treatment is strong, that is, θ
>
θ
* and μ
1
>
μ
1
*, the basic reproduction number R
0 being unity is a strict threshold for the control of the disease; the disease will be extinct or otherwise will tend to break out and persist. The other results are displayed graphically from our numerical simulation. We show the variation of the infective population and recovered population with time for different values of μ
1. It is found that both the infective population and the recovered population decrease as μ
1 increases. The infective population decreases whereas the recovered population increases as the successive vaccination rate increases, θ, respectively.



---

*Source: 678723-2014-06-16.xml*