# Stability and Hopf Bifurcation in a Delayed HIV Infection Model with General Incidence Rate and Immune Impairment

Authors: Fuxiang Li; Wanbiao Ma; Zhichao Jiang; Dan Li
Journal: Computational and Mathematical Methods in Medicine
 (2015)

---

## Abstract

We investigate the dynamical behavior of a delayed HIV infection model with general incidence rate and immune impairment. We derive two threshold parameters, the basic reproduction numberR
0 and the immune response reproduction number R
1. By using Lyapunov functional and LaSalle invariance principle, we prove the global stability of the infection-free equilibrium and the infected equilibrium without immunity. Furthermore, the existence of Hopf bifurcations at the infected equilibrium with CTL response is also studied. By theoretical analysis and numerical simulations, the effect of the immune impairment rate on the stability of the infected equilibrium with CTL response has been studied.


---

## Body

## 1. Introduction

In recent years, mathematical models have been proved to be valuable in understanding the dynamics of viral infection (see, e.g., [1–8]). In most virus infections, cytotoxic T lymphocyte (CTL) cells play a significant role in antiviral defense by attacking virus-infected cells. In order to study the role of the population dynamics of the viral infection with CTL response, Nowak and Bangham et al. proposed a basic viral infection model describing the interactions between a replicating virus population and a specific antiviral CTL response, which takes into account four populations: uninfected cells, actively infected cells, free virus, and CTL cells (see, e.g., [1–4, 9, 10]). Now, the population dynamics of viral infection with CTL response has been paid much attention and many properties have been investigated (see, e.g., [11–16]).Furthermore, the state of latent infection cannot be ignored in many biological models. The infected cells are separated into two distinct compartments, latently infected and actively infected. These latently infected cells do not produce virus and can evade from viral cytopathic effects and host immune mechanisms (see, e.g., [17–20]). Recently, the following model with latent infection and CTL response has been proposed (see, e.g., [11]):(1)
x
˙
t
=
λ
-
β
x
t
v
t
-
μ
1
x
t
,
u
˙
t
=
β
x
t
v
t
-
σ
+
μ
2
u
t
,
y
˙
t
=
σ
u
t
-
p
y
t
z
t
-
μ
3
y
t
,
v
˙
t
=
k
y
t
-
μ
4
v
t
,
z
˙
t
=
q
y
t
z
t
-
μ
5
z
t
,where x
(
t
), u
(
t
), y
(
t
), v
(
t
), and z
(
t
) represent the numbers of uninfected cells, latently infected cells, actively infected cells, free virus, and CTLs at time t, respectively. Uninfected cells are produced at the rate λ, die at the rate μ
1, and become infected at the rate β. The constant σ is the rate of latently infected cells translating to actively infected cells and μ
3 is the death rate of actively infected cells. The constant μ
2 represents the death rate of latently infected cells. The constant p is the rate of CTL-mediated lysis and q is the rate of CTL proliferation. The constant k is the rate of production of virus by infected cells and μ
4 is the clearance rate of free virus. The removal rate of CTLs is μ
5.However, in plenty of previous papers, many models are constructed under the assumption that the presence of antigen can stimulate immunity and ignore the immune impairment (see, e.g., [8, 11, 16, 17]). In fact, some pathogens can also suppress immune response or even destroy immunity especially when the load of pathogens is too high such as HIV, HBV (see, e.g., [15, 21–25]). Regoes et al. consider an ordinary differential equation (ODE) model with an immune impairment term m
y
z (see, e.g., [12, 26, 27]), where m denotes the immune impairment rate. Time delay should be considered in models for CTL response. It is shown that time delay plays an important role to the dynamic properties in models for CTL response (see, e.g., [1, 5, 6, 8, 15]). In fact, antigenic stimulation generating CTLs may need a period of time t; that is, the CTL response at time t may depend on the numbers of CTLs and infected cells at time t
-
τ, for a time lag τ
>
0 (see, e.g., [1, 5, 13]).Motivated by the above works, in this paper, we will study a delay differential equation (DDE) model of HIV infection with immune impairment and delayed CTL response. Furthermore, we know that the actual incidence rate is probably not linear over the entire range ofx and v. Based on the works mentioned above (see, e.g., [21, 28–31]), we propose the following system with general incidence function:(2)
x
˙
t
=
λ
-
f
x
t
,
v
t
v
t
-
μ
1
x
t
,
u
˙
t
=
f
x
t
,
v
t
v
t
-
σ
+
μ
2
u
t
,
y
˙
t
=
σ
u
t
-
p
y
t
z
t
-
μ
3
y
t
,
v
˙
t
=
k
y
t
-
μ
4
v
t
,
z
˙
t
=
q
y
t
-
τ
z
t
-
τ
-
μ
5
z
t
-
m
y
t
z
t
,where the state variables x
(
t
), u
(
t
), y
(
t
), v
(
t
), and z
(
t
) and the parameters λ, σ, p, k, q, μ
1, μ
2, μ
3, μ
4, and μ
5 have the same biological meaning as in system (1). m is the immune impairment rate. Suppose all the parameters are nonnegative. We assume the incidence rate is the general incidence function f
(
x
,
v
)
v, where f
∈
C
1
(
[
0
,
+
∞
]
×
[
0
,
+
∞
]
,
R
) satisfies the following hypotheses:(H1)
f
(
x
,
v
)
v
≥
0, for all x
≥
0 and v
≥
0; f
(
x
,
v
)
=
0 if and only if x
=
0;
(H2)
∂
f
(
x
,
v
)
/
∂
x
>
0, for all x
≥
0 and v
≥
0;
(H3)
∂
f
(
x
,
v
)
/
∂
v
≤
0, for all x
≥
0 and v
≥
0;
(H4)
∂
(
f
(
x
,
v
)
v
)
/
∂
v
>
0, for all x
>
0 and v
≥
0.Clearly, the hypotheses can be satisfied by different types of the incidence rate including the mass action, the Holling type II function, the saturation incidence, Beddington-DeAngelis incidence function, Crowley-Martin incidence function, and the more generalized incidence functions (see, e.g., [4, 6, 17, 32, 33]). Further, in order to study the global stability of the equilibria of system (2) by the method of Lyapunov functionals, we assume the following hypotheses hold (see, e.g., [28]):(H5)
x
-
x
0
-
∫
x
0
x
‍
f
(
x
0
,
0
)
/
f
(
s
,
0
)
d
s
→
+
∞, as x
→
+
∞ or x
→
0
+;
(H6)
x
-
x
1
-
∫
x
1
x
‍
f
(
x
1
,
v
1
)
/
f
(
s
,
v
1
)
d
s
→
+
∞, as x
→
+
∞ or x
→
0
+;
(H7)
x
-
x
∗
-
∫
x
∗
x
‍
f
(
x
∗
,
v
∗
)
/
f
(
s
,
v
∗
)
d
s
→
+
∞, as x
→
+
∞ or x
→
0
+.The main purpose of this paper is to carry out a complete theoretical analysis on the global stability of the equilibria of system (2). The organization of this paper is as follows. In Section 2, we consider the nonnegativity and boundedness of the solutions and the existence of the equilibria of system (2). In Section 3, we consider the global stability of the infection-free equilibrium E
0 and the infected equilibrium without immunity E
1 by constructing suitable Lyapunov functionals and using LaSalle invariance principle. In Section 4, we discuss the local stability of the infected equilibrium with CTL response E
∗ and the existence of Hopf bifurcations. Finally, in Section 5, the brief conclusions are given and some numerical simulations are carried out to illustrate the main results.


## 2. Basic Results

### 2.1. The Nonnegativity and Boundedness of the Solutions

According to biological meanings, the initial condition of system (2) is given as follows:(3)
x
θ
=
φ
1
θ
,
u
θ
=
φ
2
θ
,
y
θ
=
φ
3
θ
,
v
θ
=
φ
4
θ
,
z
θ
=
φ
5
θ
,where θ
∈
[
-
τ
,
0
] and (
φ
1
,
φ
2
,
φ
3
,
φ
4
,
φ
5
)
∈
C
=
C
(
[
-
τ
,
0
]
,
R
+
5
) and C is the Banach space of the continuous functions mapping the interval [
-
τ
,
0
] into R
+
5
,
R
+
5
=
x
1
,
x
2
,
x
3
,
x
4
,
x
5
∣
x
i
≥
0
,
i
=
1,2
,
3,4
,
5.Under the initial condition (3), it easily shows that the solution of system (2) is unique and nonnegative for all t
≥
0 and ultimately bounded. It has the following result.Proposition 1.
Under the initial condition (3), the solution of system (2) is unique and nonnegative for all t
≥
0 and also ultimately bounded, when (H1)–(H7) are satisfied.Proof.
The uniqueness and nonnegativity of the solution(
x
(
t
)
,
u
(
t
)
,
y
(
t
)
,
v
(
t
)
,
z
(
t
)
) can be easily proved by using the theorems in [34, 35].
Next, fort
≥
0, define (4)
L
t
=
x
t
+
u
t
+
y
t
+
μ
3
2
k
v
t
+
p
2
q
z
t
+
τ
.By the nonnegativity of the solutions, it follows that, for t
≥
0, (5)
L
′
t
=
λ
-
μ
1
x
t
-
μ
2
u
t
-
μ
3
2
y
t
-
μ
3
μ
4
2
k
v
t
-
p
μ
5
2
q
z
t
+
τ
-
p
2
y
t
z
t
-
p
m
2
q
y
t
+
τ
z
t
+
τ
≤
λ
-
γ
L
t
,where γ
=
m
i
n
{
μ
1
,
μ
2
,
μ
3
/
2
,
μ
4
,
μ
5
}. Thus, it has that l
i
m
s
u
p
t
→
+
∞
L
(
t
)
≤
λ
/
γ, from which it has that the solution (
x
(
t
)
,
u
(
t
)
,
y
(
t
)
,
v
(
t
)
,
z
(
t
)
) is ultimately bounded.


### 2.2. The Existence of the Equilibria

Next, we consider the existence of the equilibria. The equilibrium of system (2) satisfies(6)
λ
-
f
x
,
v
v
-
μ
1
x
=
0
,
f
x
,
v
v
-
σ
+
μ
2
u
=
0
,
σ
u
-
p
y
z
-
μ
3
y
=
0
,
k
y
-
μ
4
v
=
0
,
q
y
z
-
μ
5
z
-
m
y
z
=
0
.Ifu
=
0
,
y
=
0
,
v
=
0, and z
=
0, system (2) has only one equilibrium, that is, the infection-free equilibrium E
0
=
(
x
0
,
0,0
,
0,0
), where x
0
=
λ
/
μ
1.Ifu
≠
0
,
y
≠
0
,
v
≠
0, and z
=
0, we have(7)
f
x
,
k
σ
λ
-
μ
1
x
μ
3
μ
4
σ
+
μ
2
-
μ
3
μ
4
σ
+
μ
2
k
σ
=
0
,
(8)
y
=
σ
λ
-
μ
1
x
μ
3
σ
+
μ
2
,
u
=
λ
-
μ
1
x
σ
+
μ
2
,
v
=
k
σ
λ
-
μ
1
x
μ
3
μ
4
σ
+
μ
2
.Since v
>
0, we have that x
<
λ
/
μ
1. Hence, we only need to consider the case of x
<
λ
/
μ
1.Consider the following function defined on the interval(
0
,
λ
/
μ
1
) by (9)
F
x
=
f
x
,
k
σ
λ
-
μ
1
x
μ
3
μ
4
σ
+
μ
2
-
μ
3
μ
4
σ
+
μ
2
k
σ
.Under hypotheses (H2) and (H3), we have (10)
F
′
x
=
∂
f
∂
x
+
∂
f
∂
v
-
k
σ
μ
1
μ
3
μ
4
σ
+
μ
2
>
0
.We know that the function F
(
x
) is strictly monotonically increasing with respect to x. Denote the basic reproduction number R
0 of system (2) by (11)
R
0
=
k
σ
f
λ
/
μ
1
,
0
μ
3
μ
4
σ
+
μ
2
.Clearly, we have (12)
F
0
=
-
μ
3
μ
4
σ
+
μ
2
k
σ
<
0
,
F
λ
μ
1
=
f
λ
μ
1
,
0
-
μ
3
μ
4
σ
+
μ
2
k
σ
=
μ
3
μ
4
σ
+
μ
2
k
σ
R
0
-
1
.It has that there exists a unique x
1
∈
(
0
,
λ
/
μ
1
) such that F
(
x
1
)
=
0, if R
0
>
1. Then we can compute u
1
,
y
1 and v
1 by (8). Hence, we get the unique infected equilibrium without immunity E
1
=
(
x
1
,
u
1
,
y
1
,
v
1
,
0
).Ifz
≠
0 and q
>
m, we get the following equations:(13)
f
x
,
k
μ
5
μ
4
q
-
m
k
μ
5
μ
4
q
-
m
-
λ
+
μ
1
x
=
0
,
(14)
u
=
λ
-
μ
1
x
σ
+
μ
2
,
y
=
μ
5
q
-
m
>
0
,
v
=
k
μ
5
μ
4
q
-
m
>
0
,
(15)
z
=
λ
-
μ
1
x
q
-
m
σ
-
μ
3
μ
5
σ
+
μ
2
p
μ
5
σ
+
μ
2
.Sincez
>
0, we have x
<
x
¯, where (16)
x
¯
=
λ
q
-
m
σ
-
μ
3
μ
5
σ
+
μ
2
μ
1
q
-
m
σ
.Hence, the existence of the equilibrium requires x
¯
>
0 and (13) has a solution on the interval (
0
,
x
¯
).Denote(17)
R
¯
=
λ
q
-
m
σ
μ
3
μ
5
σ
+
μ
2
.Hence, if R
¯
>
1, it has x
¯
>
0. Denote (18)
G
x
=
f
x
,
k
μ
5
μ
4
q
-
m
k
μ
5
μ
4
q
-
m
-
λ
+
μ
1
x
.Under hypothesis (H2), we know that the function G
(
x
) is strictly monotonically increasing with respect to x. Clearly, we have (19)
G
0
=
-
λ
<
0
,
G
x
¯
=
f
x
¯
,
k
μ
5
μ
4
q
-
m
k
μ
5
μ
4
q
-
m
-
λ
+
μ
1
x
¯
=
f
x
¯
,
k
μ
5
μ
4
q
-
m
k
μ
5
μ
4
q
-
m
-
μ
3
μ
5
σ
+
μ
2
q
-
m
σ
=
μ
3
μ
5
σ
+
μ
2
q
-
m
σ
R
1
-
1
,where (20)
R
1
=
k
σ
f
x
¯
,
k
μ
5
/
μ
4
q
-
m
μ
3
μ
4
σ
+
μ
2
.Hence, we have that there exists x
∗
∈
(
0
,
x
¯
) such that G
(
x
∗
)
=
0, if R
¯
>
1 and R
1
>
1. Then we can compute u
∗, y
∗, v
∗, and z
∗ by (14) and (15).Denote the immune response reproduction number of system (2) as R
1. Therefore, we have that there exists a unique infected equilibrium with CTL response E
∗
=
(
x
∗
,
u
∗
,
y
∗
,
v
∗
,
z
∗
), if R
¯
>
1 and R
1
>
1. This proves the following theorem.Theorem 2.
Suppose that hypotheses (H1)–(H4) are satisfied; the following conclusions hold.(i)
System (2) always has an infection-free equilibrium E
0.
(ii)
System (2) has an infected equilibrium without immunity E
1 if R
0
>
1.
(iii)
System (2) has an infected equilibrium with immunity E
∗ if R
¯
>
1 and R
1
>
1.From hypotheses (H1)–(H3), it is clear that R
1
<
R
0. In order to study the global stability of the infected equilibrium E
1 in the next section, we give the following remark.Remark 3.
Suppose thatR
¯
>
1 is satisfied; then the following results hold:(i)
IfR
1
>
1, then (
q
-
m
)
y
1
/
μ
5
>
1.
(ii)
IfR
1
≤
1, then (
q
-
m
)
y
1
/
μ
5
≤
1.Let us give the proof of Remark3. Firstly, for Case (i), since R
1
>
1, then (21)
F
x
¯
=
f
x
¯
,
k
σ
λ
-
μ
1
x
¯
μ
3
μ
4
σ
+
μ
2
-
μ
3
μ
4
σ
+
μ
2
k
σ
=
μ
3
μ
4
σ
+
μ
2
k
σ
R
1
-
1
>
0
.Since the function F
(
x
) is strictly monotonically increasing with respect to x and F
(
x
1
)
=
0, we have x
1
<
x
¯. Therefore (22)
λ
-
μ
1
x
1
>
λ
-
μ
1
x
¯
=
μ
3
μ
5
σ
+
μ
2
σ
q
-
m
.Then (23)
q
-
m
y
1
μ
5
=
q
-
m
μ
5
·
σ
λ
-
μ
1
x
1
μ
3
σ
+
μ
2
>
1
.Secondly, for Case (ii), sinceR
1
≤
1, then (24)
F
x
¯
=
f
x
¯
,
k
σ
λ
-
μ
1
x
¯
μ
3
μ
4
σ
+
μ
2
-
μ
3
μ
4
σ
+
μ
2
k
σ
=
μ
3
μ
4
σ
+
μ
2
k
σ
R
1
-
1
≤
0
.We have x
1
≥
x
¯. Therefore (25)
λ
-
μ
1
x
1
≤
λ
-
μ
1
x
¯
=
μ
3
μ
5
σ
+
μ
2
σ
q
-
m
.Then (26)
q
-
m
y
1
μ
5
=
q
-
m
μ
5
·
σ
λ
-
μ
1
x
1
μ
3
σ
+
μ
2
≤
1
.


## 2.1. The Nonnegativity and Boundedness of the Solutions

According to biological meanings, the initial condition of system (2) is given as follows:(3)
x
θ
=
φ
1
θ
,
u
θ
=
φ
2
θ
,
y
θ
=
φ
3
θ
,
v
θ
=
φ
4
θ
,
z
θ
=
φ
5
θ
,where θ
∈
[
-
τ
,
0
] and (
φ
1
,
φ
2
,
φ
3
,
φ
4
,
φ
5
)
∈
C
=
C
(
[
-
τ
,
0
]
,
R
+
5
) and C is the Banach space of the continuous functions mapping the interval [
-
τ
,
0
] into R
+
5
,
R
+
5
=
x
1
,
x
2
,
x
3
,
x
4
,
x
5
∣
x
i
≥
0
,
i
=
1,2
,
3,4
,
5.Under the initial condition (3), it easily shows that the solution of system (2) is unique and nonnegative for all t
≥
0 and ultimately bounded. It has the following result.Proposition 1.
Under the initial condition (3), the solution of system (2) is unique and nonnegative for all t
≥
0 and also ultimately bounded, when (H1)–(H7) are satisfied.Proof.
The uniqueness and nonnegativity of the solution(
x
(
t
)
,
u
(
t
)
,
y
(
t
)
,
v
(
t
)
,
z
(
t
)
) can be easily proved by using the theorems in [34, 35].
Next, fort
≥
0, define (4)
L
t
=
x
t
+
u
t
+
y
t
+
μ
3
2
k
v
t
+
p
2
q
z
t
+
τ
.By the nonnegativity of the solutions, it follows that, for t
≥
0, (5)
L
′
t
=
λ
-
μ
1
x
t
-
μ
2
u
t
-
μ
3
2
y
t
-
μ
3
μ
4
2
k
v
t
-
p
μ
5
2
q
z
t
+
τ
-
p
2
y
t
z
t
-
p
m
2
q
y
t
+
τ
z
t
+
τ
≤
λ
-
γ
L
t
,where γ
=
m
i
n
{
μ
1
,
μ
2
,
μ
3
/
2
,
μ
4
,
μ
5
}. Thus, it has that l
i
m
s
u
p
t
→
+
∞
L
(
t
)
≤
λ
/
γ, from which it has that the solution (
x
(
t
)
,
u
(
t
)
,
y
(
t
)
,
v
(
t
)
,
z
(
t
)
) is ultimately bounded.


## 2.2. The Existence of the Equilibria

Next, we consider the existence of the equilibria. The equilibrium of system (2) satisfies(6)
λ
-
f
x
,
v
v
-
μ
1
x
=
0
,
f
x
,
v
v
-
σ
+
μ
2
u
=
0
,
σ
u
-
p
y
z
-
μ
3
y
=
0
,
k
y
-
μ
4
v
=
0
,
q
y
z
-
μ
5
z
-
m
y
z
=
0
.Ifu
=
0
,
y
=
0
,
v
=
0, and z
=
0, system (2) has only one equilibrium, that is, the infection-free equilibrium E
0
=
(
x
0
,
0,0
,
0,0
), where x
0
=
λ
/
μ
1.Ifu
≠
0
,
y
≠
0
,
v
≠
0, and z
=
0, we have(7)
f
x
,
k
σ
λ
-
μ
1
x
μ
3
μ
4
σ
+
μ
2
-
μ
3
μ
4
σ
+
μ
2
k
σ
=
0
,
(8)
y
=
σ
λ
-
μ
1
x
μ
3
σ
+
μ
2
,
u
=
λ
-
μ
1
x
σ
+
μ
2
,
v
=
k
σ
λ
-
μ
1
x
μ
3
μ
4
σ
+
μ
2
.Since v
>
0, we have that x
<
λ
/
μ
1. Hence, we only need to consider the case of x
<
λ
/
μ
1.Consider the following function defined on the interval(
0
,
λ
/
μ
1
) by (9)
F
x
=
f
x
,
k
σ
λ
-
μ
1
x
μ
3
μ
4
σ
+
μ
2
-
μ
3
μ
4
σ
+
μ
2
k
σ
.Under hypotheses (H2) and (H3), we have (10)
F
′
x
=
∂
f
∂
x
+
∂
f
∂
v
-
k
σ
μ
1
μ
3
μ
4
σ
+
μ
2
>
0
.We know that the function F
(
x
) is strictly monotonically increasing with respect to x. Denote the basic reproduction number R
0 of system (2) by (11)
R
0
=
k
σ
f
λ
/
μ
1
,
0
μ
3
μ
4
σ
+
μ
2
.Clearly, we have (12)
F
0
=
-
μ
3
μ
4
σ
+
μ
2
k
σ
<
0
,
F
λ
μ
1
=
f
λ
μ
1
,
0
-
μ
3
μ
4
σ
+
μ
2
k
σ
=
μ
3
μ
4
σ
+
μ
2
k
σ
R
0
-
1
.It has that there exists a unique x
1
∈
(
0
,
λ
/
μ
1
) such that F
(
x
1
)
=
0, if R
0
>
1. Then we can compute u
1
,
y
1 and v
1 by (8). Hence, we get the unique infected equilibrium without immunity E
1
=
(
x
1
,
u
1
,
y
1
,
v
1
,
0
).Ifz
≠
0 and q
>
m, we get the following equations:(13)
f
x
,
k
μ
5
μ
4
q
-
m
k
μ
5
μ
4
q
-
m
-
λ
+
μ
1
x
=
0
,
(14)
u
=
λ
-
μ
1
x
σ
+
μ
2
,
y
=
μ
5
q
-
m
>
0
,
v
=
k
μ
5
μ
4
q
-
m
>
0
,
(15)
z
=
λ
-
μ
1
x
q
-
m
σ
-
μ
3
μ
5
σ
+
μ
2
p
μ
5
σ
+
μ
2
.Sincez
>
0, we have x
<
x
¯, where (16)
x
¯
=
λ
q
-
m
σ
-
μ
3
μ
5
σ
+
μ
2
μ
1
q
-
m
σ
.Hence, the existence of the equilibrium requires x
¯
>
0 and (13) has a solution on the interval (
0
,
x
¯
).Denote(17)
R
¯
=
λ
q
-
m
σ
μ
3
μ
5
σ
+
μ
2
.Hence, if R
¯
>
1, it has x
¯
>
0. Denote (18)
G
x
=
f
x
,
k
μ
5
μ
4
q
-
m
k
μ
5
μ
4
q
-
m
-
λ
+
μ
1
x
.Under hypothesis (H2), we know that the function G
(
x
) is strictly monotonically increasing with respect to x. Clearly, we have (19)
G
0
=
-
λ
<
0
,
G
x
¯
=
f
x
¯
,
k
μ
5
μ
4
q
-
m
k
μ
5
μ
4
q
-
m
-
λ
+
μ
1
x
¯
=
f
x
¯
,
k
μ
5
μ
4
q
-
m
k
μ
5
μ
4
q
-
m
-
μ
3
μ
5
σ
+
μ
2
q
-
m
σ
=
μ
3
μ
5
σ
+
μ
2
q
-
m
σ
R
1
-
1
,where (20)
R
1
=
k
σ
f
x
¯
,
k
μ
5
/
μ
4
q
-
m
μ
3
μ
4
σ
+
μ
2
.Hence, we have that there exists x
∗
∈
(
0
,
x
¯
) such that G
(
x
∗
)
=
0, if R
¯
>
1 and R
1
>
1. Then we can compute u
∗, y
∗, v
∗, and z
∗ by (14) and (15).Denote the immune response reproduction number of system (2) as R
1. Therefore, we have that there exists a unique infected equilibrium with CTL response E
∗
=
(
x
∗
,
u
∗
,
y
∗
,
v
∗
,
z
∗
), if R
¯
>
1 and R
1
>
1. This proves the following theorem.Theorem 2.
Suppose that hypotheses (H1)–(H4) are satisfied; the following conclusions hold.(i)
System (2) always has an infection-free equilibrium E
0.
(ii)
System (2) has an infected equilibrium without immunity E
1 if R
0
>
1.
(iii)
System (2) has an infected equilibrium with immunity E
∗ if R
¯
>
1 and R
1
>
1.From hypotheses (H1)–(H3), it is clear that R
1
<
R
0. In order to study the global stability of the infected equilibrium E
1 in the next section, we give the following remark.Remark 3.
Suppose thatR
¯
>
1 is satisfied; then the following results hold:(i)
IfR
1
>
1, then (
q
-
m
)
y
1
/
μ
5
>
1.
(ii)
IfR
1
≤
1, then (
q
-
m
)
y
1
/
μ
5
≤
1.Let us give the proof of Remark3. Firstly, for Case (i), since R
1
>
1, then (21)
F
x
¯
=
f
x
¯
,
k
σ
λ
-
μ
1
x
¯
μ
3
μ
4
σ
+
μ
2
-
μ
3
μ
4
σ
+
μ
2
k
σ
=
μ
3
μ
4
σ
+
μ
2
k
σ
R
1
-
1
>
0
.Since the function F
(
x
) is strictly monotonically increasing with respect to x and F
(
x
1
)
=
0, we have x
1
<
x
¯. Therefore (22)
λ
-
μ
1
x
1
>
λ
-
μ
1
x
¯
=
μ
3
μ
5
σ
+
μ
2
σ
q
-
m
.Then (23)
q
-
m
y
1
μ
5
=
q
-
m
μ
5
·
σ
λ
-
μ
1
x
1
μ
3
σ
+
μ
2
>
1
.Secondly, for Case (ii), sinceR
1
≤
1, then (24)
F
x
¯
=
f
x
¯
,
k
σ
λ
-
μ
1
x
¯
μ
3
μ
4
σ
+
μ
2
-
μ
3
μ
4
σ
+
μ
2
k
σ
=
μ
3
μ
4
σ
+
μ
2
k
σ
R
1
-
1
≤
0
.We have x
1
≥
x
¯. Therefore (25)
λ
-
μ
1
x
1
≤
λ
-
μ
1
x
¯
=
μ
3
μ
5
σ
+
μ
2
σ
q
-
m
.Then (26)
q
-
m
y
1
μ
5
=
q
-
m
μ
5
·
σ
λ
-
μ
1
x
1
μ
3
σ
+
μ
2
≤
1
.


## 3. The Global Stability of the Equilibria

In this section, we study the global stability of the equilibria of system (2). Firstly, we analyze the global stability of the infection-free equilibrium E
0.Theorem 4.
Suppose that hypotheses (H1)–(H7) are satisfied. If R
0
≤
1, then the infection-free equilibrium E
0 is globally asymptotically stable for any time delay τ
≥
0. If R
0
>
1, then the infection-free equilibrium E
0 is unstable for any time delay τ
≥
0.Proof.
Let(
x
(
t
), u
(
t
), y
(
t
), v
t, z
(
t
)
) be a positive solution of system (2) with the initial condition (3) for t
≥
0. Motivated by the works in [14, 28, 31, 36, 37], we consider the following Lyapunov functional:(27)
V
1
=
x
-
x
0
-
∫
x
0
x
f
x
0
,
0
f
s
,
0
d
s
+
u
+
σ
+
μ
2
σ
y
+
μ
3
σ
+
μ
2
k
σ
v
+
σ
+
μ
2
σ
p
q
-
m
z
+
σ
+
μ
2
σ
p
q
-
m
∫
t
-
τ
t
q
y
θ
z
θ
d
θ
,where λ
=
μ
1
x
0. By (H1)–(H5), it is obvious that V
1 is positive definite with respect to E
0. For t
≥
0, the time derivative of V
1 along the solutions of system (2) is (28)
V
˙
1
=
1
-
f
x
0
,
0
f
x
,
0
x
˙
+
u
˙
+
σ
+
μ
2
σ
y
˙
+
μ
3
σ
+
μ
2
k
σ
v
˙
+
σ
+
μ
2
σ
p
q
-
m
z
˙
+
σ
+
μ
2
σ
p
q
-
m
q
y
t
z
t
-
q
y
t
-
τ
z
t
-
τ
=
μ
1
1
-
f
x
0
,
0
f
x
,
0
x
0
-
x
+
f
x
0
,
0
f
x
,
0
f
x
,
v
v
-
μ
3
μ
4
σ
+
μ
2
k
σ
v
-
σ
+
μ
2
σ
p
q
-
m
μ
5
z
=
μ
1
1
-
f
x
0
,
0
f
x
,
0
x
0
-
x
-
μ
3
μ
4
σ
+
μ
2
k
σ
1
-
f
x
,
v
f
x
,
0
R
0
v
-
σ
+
μ
2
σ
p
q
-
m
μ
5
z
.
Since hypotheses(
H
1
)
–
(
H
3
) and R
0
≤
1, we have(29)
μ
1
1
-
f
x
0
,
0
f
x
,
0
x
0
-
x
≤
0
,
1
-
f
x
,
v
f
x
,
0
R
0
≥
0
.Therefore, V
˙
1
≤
0 if R
0
≤
1. Then it follows from stability theorems in [34, 35] that the infection-free equilibrium E
0 is stable for any time delay τ
≥
0 if R
0
≤
1.
Furthermore, note that, for eacht
≥
0, V
˙
1
=
0 implies that x
(
t
)
=
x
0, z
(
t
)
=
0. Let M be the largest invariant set in the set (30)
Γ
1
=
φ
1
,
φ
2
,
φ
3
,
φ
4
,
φ
5
∈
C
∣
V
˙
1
=
0
⊂
φ
1
,
φ
2
,
φ
3
,
φ
4
,
φ
5
∈
C
∣
φ
1
0
=
x
0
,
φ
5
0
=
0
.We have from the first four equations of system (2) and the invariance of M that M
=
{
E
0
}. Since any solution of system (2) is bounded, it follows from LaSalle invariance principle (see, e.g., [34, 35]) that the infection-free equilibrium E
0 is also globally attractive for any time delay τ
≥
0 if R
0
≤
1.
The characteristic equation of system (2) at the infection-free equilibrium E
0 is(31)
s
+
μ
1
s
+
μ
5
s
3
+
σ
+
μ
2
+
μ
3
+
μ
4
s
2
+
μ
2
+
σ
μ
3
+
μ
4
+
μ
3
μ
4
s
+
μ
3
μ
4
σ
+
μ
2
1
-
R
0
=
0
.Clearly, if R
0
>
1, (31) has at least a positive real root. Thus, the infection-free equilibrium E
0 is unstable.Next we study the global stability of the infected equilibrium without immunityE
1.Theorem 5.
Suppose that hypotheses (H1)–(H7) and R
¯
>
1 are satisfied. If R
0
>
1
≥
R
1, then the infected equilibrium without immunity E
1 is globally asymptotically stable for any time delay τ
≥
0. If R
1
>
1, then the infected equilibrium without immunity E
1 is unstable for any time delay τ
≥
0.Proof.
Let(
x
(
t
), u
(
t
), y
(
t
), v
t, z
(
t
)
) be a positive solution of system (2) with the initial condition (3) for t
≥
0. Consider the following Lyapunov functional: (32)
V
2
=
x
-
x
1
-
∫
x
1
x
f
x
1
,
v
1
f
s
,
v
1
d
s
+
u
-
u
1
-
u
1
ln
⁡
u
u
1
+
σ
+
μ
2
σ
y
-
y
1
-
y
1
ln
⁡
y
y
1
+
μ
3
σ
+
μ
2
k
σ
v
-
v
1
-
v
1
ln
⁡
v
v
1
+
σ
+
μ
2
σ
p
q
-
m
z
+
σ
+
μ
2
σ
p
q
q
-
m
∫
t
-
τ
t
y
θ
z
θ
d
θ
.
Letψ
(
x
)
=
x
-
x
1
-
∫
x
1
x
‍
f
(
x
1
,
v
1
)
/
f
(
s
,
v
1
)
d
s. Then, ψ
(
x
) has the global minimum at x
=
x
1 and ψ
(
x
1
)
=
0. Furthermore, ψ
(
x
)
>
0 for x
>
0. Hence, V
2 is positive definite with respect to E
1. For t
≥
0, the time derivative of V
2 along the solutions of system (2) is(33)
V
˙
2
=
1
-
f
x
1
,
v
1
f
x
,
v
1
x
˙
+
1
-
u
1
u
u
˙
+
σ
+
μ
2
σ
1
-
y
1
y
y
˙
+
μ
3
σ
+
μ
2
k
σ
1
-
v
1
v
v
˙
+
σ
+
μ
2
σ
p
q
-
m
z
˙
+
σ
+
μ
2
σ
·
q
q
-
m
p
y
t
z
t
-
p
y
t
-
τ
z
t
-
τ
=
1
-
f
x
1
,
v
1
f
x
,
v
1
λ
-
f
x
,
v
v
-
μ
1
x
+
1
-
u
1
u
f
x
,
v
v
-
σ
+
μ
2
u
+
σ
+
μ
2
σ
1
-
y
1
y
σ
u
-
p
y
z
-
μ
3
y
+
μ
3
σ
+
μ
2
k
σ
1
-
v
1
v
k
y
-
μ
4
v
+
σ
+
μ
2
σ
·
p
q
-
m
q
y
t
-
τ
z
t
-
τ
-
μ
5
z
-
m
y
z
+
σ
+
μ
2
σ
·
q
q
-
m
p
y
t
z
t
-
p
y
t
-
τ
z
t
-
τ
.
Note thatλ
=
f
(
x
1
,
v
1
)
v
1
+
μ
1
x
1, f
(
x
1
,
v
1
)
v
1
=
(
σ
+
μ
2
)
u
1, and μ
3
y
1
=
σ
u
1; we have (34)
V
˙
2
=
1
-
f
x
1
,
v
1
f
x
,
v
1
f
x
1
,
v
1
v
1
+
μ
1
x
1
-
f
x
,
v
v
-
μ
1
x
+
1
-
u
1
u
f
x
,
v
v
-
f
x
1
,
v
1
v
1
u
1
u
+
σ
+
μ
2
σ
1
-
y
1
y
-
p
y
z
+
f
x
1
,
v
1
v
1
1
-
y
1
y
u
u
1
-
y
y
1
+
μ
3
σ
+
μ
2
k
σ
1
-
v
1
v
k
y
-
μ
4
v
+
σ
+
μ
2
σ
·
p
q
-
m
q
y
t
-
τ
z
t
-
τ
-
μ
5
z
-
m
y
z
+
σ
+
μ
2
σ
q
q
-
m
p
y
t
z
t
-
p
y
t
-
τ
z
t
-
τ
=
μ
1
1
-
f
x
1
,
v
1
f
x
,
v
1
x
1
-
x
+
σ
+
μ
2
σ
p
z
y
1
-
μ
5
q
-
m
+
f
x
1
,
v
1
v
1
5
-
y
1
y
u
u
1
-
f
x
,
v
v
f
x
1
,
v
1
v
1
u
1
u
-
y
y
1
v
1
v
-
f
x
1
,
v
1
f
x
,
v
1
-
f
x
,
v
1
f
x
,
v
+
f
x
1
,
v
1
v
1
-
1
-
v
v
1
+
f
x
,
v
v
f
x
,
v
1
v
1
+
f
x
,
v
1
f
x
,
v
.
Since the arithmetic mean is greater than or equal to the geometric mean, it has(35)
5
-
y
1
y
u
u
1
-
f
x
,
v
v
f
x
1
,
v
1
v
1
u
1
u
-
y
y
1
v
1
v
-
f
x
1
,
v
1
f
x
,
v
1
-
f
x
,
v
1
f
x
,
v
≤
0
.From hypotheses (H3)-(H4), we have (36)
-
1
-
v
v
1
+
f
x
,
v
v
f
x
,
v
1
v
1
+
f
x
,
v
1
f
x
,
v
=
f
x
,
v
-
f
x
,
v
1
f
x
,
v
1
f
x
,
v
v
-
f
x
,
v
1
v
1
f
x
,
v
v
1
≤
0
.
Note Remark3, we have y
1
≤
μ
5
/
(
q
-
m
). Therefore, V
˙
2
≤
0 if R
1
≤
1. Then it follows from stability theorems in [34, 35] that the infected equilibrium without immunity E
1 is stable for any time delay τ
≥
0 if R
1
≤
1.
Furthermore, note that, for eacht
≥
0, V
˙
2
=
0 implies that x
(
t
)
=
x
1, u
(
t
)
=
u
1, y
(
t
)
=
y
1, and v
(
t
)
=
v
1. Let M be the largest invariant set in the set (37)
Γ
2
=
φ
1
,
φ
2
,
φ
3
,
φ
4
,
φ
5
∈
C
∣
V
˙
2
=
0
⊂
φ
1
,
φ
2
,
φ
3
,
φ
4
,
φ
5
∈
C
∣
φ
1
0
=
x
1
,
φ
2
0
=
u
1
,
φ
3
0
=
y
1
,
φ
4
0
=
v
1
.We have from system (2) and the invariance of M that M
=
{
E
1
}. Since any solution of system (2) is bounded, it follows from LaSalle invariance principle (see, e.g., [34, 35]) that the infected equilibrium without immunity E
1 is also globally attractive for any time delay τ
≥
0 if R
1
≤
1.
The characteristic equation of system (2) at E
1 takes the form(38)
s
+
μ
5
+
m
y
1
-
q
y
1
e
-
s
τ
ψ
0
s
=
0
,where ψ
0
(
s
) is a polynomial with respect to s. Let (39)
ψ
1
s
=
s
+
μ
5
+
m
y
1
-
q
y
1
e
-
s
τ
.Thus we have l
i
m
s
→
+
∞
ψ
1
(
s
)
>
0 and ψ
1
(
0
)
=
μ
5
-
(
q
-
m
)
y
1. From Remark 3, we have that (
q
-
m
)
y
1
/
μ
5
>
1 if R
1
>
1. Thus, ψ
1
(
0
)
<
0 if R
1
>
1. Hence, if R
1
>
1, then ψ
1
(
s
)
=
0 has at least a positive real root; that is, (38) has at least a positive real root. Therefore, the infected equilibrium without immunity E
1 is unstable.


## 4. The Local Stability of the Infected Equilibrium and Hopf Bifurcation

The characteristic equation of system (2) at the infected equilibrium with CTL response E
∗ is given by(40)
s
5
+
A
1
s
4
+
A
2
s
3
+
A
3
s
2
+
A
4
s
+
A
5
+
e
-
s
τ
B
1
s
4
+
B
2
s
3
+
B
3
s
2
+
B
4
s
+
B
5
=
0
,where(41)
A
1
=
A
+
D
+
E
+
μ
4
+
q
y
∗
,
A
2
=
E
+
μ
4
q
y
∗
+
μ
4
E
+
D
E
+
μ
4
+
q
y
∗
+
A
D
+
E
+
μ
4
+
q
y
∗
-
p
y
∗
m
z
∗
,
A
3
=
μ
4
E
q
y
∗
+
D
μ
4
+
E
q
y
∗
-
σ
k
B
+
A
μ
4
+
E
q
y
∗
+
μ
4
E
+
D
E
+
μ
4
+
q
y
∗
-
p
y
∗
m
z
∗
A
+
D
+
μ
4
,
A
4
=
-
B
k
σ
q
y
∗
+
A
μ
4
E
q
y
∗
+
D
μ
4
+
E
q
y
∗
-
B
k
σ
+
C
k
σ
B
+
H
-
p
y
∗
m
z
∗
D
μ
4
+
A
μ
4
+
A
D
,
A
5
=
C
σ
k
B
+
H
q
y
∗
-
p
y
∗
m
z
∗
A
D
μ
4
-
A
B
σ
k
q
y
∗
,
B
1
=
-
q
y
∗
,
B
2
=
-
q
y
∗
A
+
D
+
μ
3
+
μ
4
,
B
3
=
-
q
y
∗
μ
3
μ
4
+
D
μ
3
+
μ
4
+
A
D
+
μ
3
+
μ
4
,
B
4
=
-
q
y
∗
D
μ
3
μ
4
+
A
μ
3
μ
4
+
D
μ
3
+
D
μ
4
-
σ
k
B
+
H
,
B
5
=
-
q
y
∗
A
D
μ
3
μ
4
-
A
σ
k
B
+
H
+
C
σ
k
B
+
H
,
a
=
∂
f
x
∗
,
v
∗
∂
x
>
0
,
d
=
∂
f
x
∗
,
v
∗
∂
v
≤
0
,
A
=
a
v
∗
+
μ
1
,
B
=
d
v
∗
,
C
=
a
v
∗
,
D
=
σ
+
μ
2
,
E
=
μ
3
+
p
z
∗
,
F
=
q
-
m
z
∗
,
H
=
f
x
∗
,
v
∗
,
M
=
p
y
∗
F
.When τ
=
0, (40) becomes(42)
s
5
+
α
1
s
4
+
α
2
s
3
+
α
3
s
2
+
α
4
s
+
α
5
=
0
,where (43)
α
1
=
A
1
+
B
1
=
A
+
D
+
E
+
μ
4
>
0
,
α
2
=
A
2
+
B
2
=
M
+
A
D
+
E
+
μ
4
+
E
μ
4
+
D
μ
4
+
D
E
>
0
,
α
3
=
A
3
+
B
3
=
A
+
D
+
μ
4
M
+
A
E
μ
4
+
D
μ
4
+
D
E
-
B
k
σ
>
0
,
α
4
=
A
4
+
B
4
=
D
μ
4
+
A
μ
4
+
A
D
M
+
G
>
0
,
α
5
=
A
5
+
B
5
=
A
D
μ
4
M
>
0
,
G
=
C
D
E
μ
4
-
μ
1
B
k
σ
>
0
.Denote (44)
Δ
1
=
α
1
,
Δ
2
=
α
1
α
2
-
α
3
,
Δ
3
=
α
3
Δ
2
+
α
1
α
5
-
α
1
2
α
4
,
Δ
4
=
α
4
Δ
3
+
α
1
α
5
-
α
2
α
5
Δ
2
-
α
5
2
,
Δ
5
=
α
5
Δ
4
.Since μ
4
(
σ
+
μ
2
)
(
μ
3
+
p
z
∗
)
=
σ
k
f
(
x
∗
,
v
∗
) and (
H
4
), we have σ
k
H
=
μ
4
D
E, μ
4
D
E
+
B
k
σ
>
0, and A
D
E
μ
4
-
G
>
0
.Thus,(45)
Δ
1
=
α
1
>
0
,
Δ
2
=
E
M
+
A
2
D
+
E
+
μ
4
+
D
+
E
+
μ
4
A
D
+
E
+
μ
4
+
D
E
+
μ
4
+
E
μ
4
+
B
k
σ
>
0
,
Δ
3
=
M
A
A
2
E
+
A
E
2
+
E
2
μ
4
+
D
E
μ
4
+
E
M
+
B
k
σ
+
D
A
2
E
+
D
2
E
+
A
D
E
+
A
E
2
+
D
E
2
+
D
E
μ
4
+
E
M
+
B
k
σ
+
μ
4
A
2
E
+
D
2
E
+
D
E
μ
4
+
A
D
E
+
A
E
2
+
D
E
2
+
D
E
μ
4
+
E
2
μ
4
+
E
μ
4
2
+
E
M
+
B
k
σ
+
μ
1
D
E
μ
4
+
B
k
σ
A
+
D
+
E
+
μ
4
2
+
A
D
E
A
2
D
+
A
2
E
+
A
D
2
+
A
D
E
+
D
2
E
+
D
E
μ
4
+
A
D
E
+
A
E
2
+
D
E
2
+
E
M
+
B
k
σ
+
A
D
μ
4
A
2
D
+
A
2
E
+
A
2
μ
4
+
A
D
2
+
A
D
E
+
A
D
μ
4
+
D
2
E
+
D
2
μ
4
+
A
D
E
+
A
E
2
+
A
E
μ
4
+
D
E
2
+
D
E
μ
4
+
A
D
μ
4
+
A
μ
4
2
+
D
E
μ
4
+
D
μ
4
2
+
E
M
+
B
k
σ
+
A
E
μ
4
-
B
k
σ
Δ
2
>
0
,
Δ
4
=
N
1
M
3
+
N
2
M
2
+
N
3
M
+
N
4
,where (46)
N
1
=
E
A
2
D
+
μ
4
+
D
+
μ
4
A
μ
4
+
D
μ
4
+
A
D
>
0
,
N
2
=
A
2
D
A
2
E
+
A
E
2
+
A
2
μ
4
A
2
E
+
A
E
2
+
D
2
μ
4
D
2
E
+
D
E
2
+
D
E
μ
4
+
B
k
σ
+
A
D
μ
4
A
2
E
+
D
2
E
+
D
E
2
+
A
D
2
A
2
E
+
D
2
E
+
A
D
E
+
A
E
2
+
D
E
2
+
D
E
μ
4
+
E
2
μ
4
+
D
E
μ
4
+
B
k
σ
+
D
μ
4
2
D
2
E
+
D
E
2
+
D
E
μ
4
+
E
2
μ
4
+
D
E
μ
4
+
E
μ
4
2
+
B
k
σ
+
A
μ
4
2
+
a
d
μ
4
A
2
E
+
D
2
E
+
A
E
2
+
A
D
E
+
D
E
2
+
D
E
μ
4
+
E
2
μ
4
+
A
E
μ
4
+
D
E
μ
4
+
E
μ
4
2
+
B
k
σ
+
G
E
A
+
D
+
μ
4
+
A
2
E
2
D
2
+
D
μ
4
+
μ
4
2
-
B
k
σ
E
A
D
+
A
μ
4
+
D
μ
4
+
A
A
D
+
A
μ
4
+
D
μ
4
D
E
μ
4
+
B
k
σ
>
0
,
N
3
=
A
2
E
2
+
A
3
E
+
D
3
E
+
A
D
E
2
+
D
2
E
2
+
A
E
2
μ
4
+
E
2
μ
4
2
+
E
μ
4
3
+
A
+
D
+
μ
4
B
k
σ
-
E
B
k
σ
-
B
k
σ
A
D
+
A
μ
4
+
D
μ
4
Δ
2
-
E
M
+
A
2
E
μ
4
2
+
A
2
D
2
E
A
D
E
+
D
E
μ
4
+
D
2
E
+
A
D
E
+
A
E
2
+
A
E
μ
4
+
D
E
μ
4
+
D
E
2
+
A
D
μ
4
+
A
E
μ
4
+
E
μ
4
2
+
D
E
μ
4
+
A
2
D
+
A
2
E
+
A
2
μ
4
+
B
k
σ
+
A
2
D
E
μ
4
A
D
E
+
D
E
μ
4
+
A
E
2
+
A
E
μ
4
+
A
2
E
+
B
k
σ
+
A
3
E
μ
4
4
+
A
2
D
2
E
A
D
2
+
E
2
μ
4
-
G
2
D
2
μ
4
2
+
2
A
D
2
μ
4
+
2
A
D
μ
4
2
+
D
3
μ
4
+
D
μ
4
3
,
N
4
=
G
A
E
μ
4
-
B
k
σ
D
+
E
+
μ
4
·
A
D
+
A
E
+
A
μ
4
+
E
μ
4
+
D
μ
4
+
D
E
+
A
2
D
+
E
+
μ
4
+
B
k
σ
+
G
A
D
E
A
D
2
+
A
D
E
+
D
E
μ
4
+
D
2
E
+
A
D
E
+
A
E
2
+
D
E
2
+
A
2
D
+
A
2
E
+
B
k
σ
+
G
A
D
μ
4
A
D
2
+
A
D
E
+
D
E
μ
4
+
D
2
E
+
A
D
E
+
A
E
2
+
D
E
2
+
A
2
D
+
A
2
E
+
A
D
μ
4
+
D
2
μ
4
+
A
E
μ
4
+
A
D
μ
4
+
A
μ
4
2
+
D
μ
4
2
+
D
E
μ
4
+
A
2
μ
4
+
B
k
σ
>
0
.Assume further that(H8)
E
≥
m
a
x
{
μ
4
,
D
}; that is, μ
3
+
p
z
∗
≥
μ
4 and μ
3
+
p
z
∗
≥
σ
+
μ
2
.We have(47)
N
3
≥
G
A
2
E
2
+
A
3
E
-
E
B
k
σ
-
B
k
σ
A
μ
4
+
D
μ
4
·
Δ
2
-
E
M
-
B
k
σ
A
D
Δ
2
-
A
E
μ
4
-
D
E
μ
4
-
E
μ
4
2
-
B
k
σ
A
D
E
μ
4
-
G
A
+
D
+
μ
4
+
A
2
E
μ
4
2
+
A
2
D
2
E
A
D
E
+
D
E
μ
4
+
D
2
E
+
A
D
E
+
A
E
2
+
A
E
μ
4
+
D
E
μ
4
+
D
E
2
+
A
D
μ
4
+
A
E
μ
4
+
E
μ
4
2
+
D
E
μ
4
+
A
2
D
+
A
2
E
+
A
2
μ
4
+
B
k
σ
+
A
2
D
E
μ
4
A
D
E
+
D
E
μ
4
+
A
E
2
+
A
E
μ
4
+
A
2
E
+
B
k
σ
+
A
3
E
μ
4
4
+
A
2
D
2
E
A
D
2
+
E
2
μ
4
>
0
.Therefore, Δ
4
>
0, Δ
5
>
0
. By Routh-Hurwitz criterion, all the roots of (42) have negative real parts. Hence we have the following result.Proposition 6.
Whenτ
=
0, if R
¯
>
1, R
1
>
1, and (H8) hold, then the infected equilibrium with CTL response E
∗ is locally asymptotically stable.In fact, whenτ
=
0, we can show that if R
¯
>
1 and R
1
>
1 hold, the infected equilibrium with CTL response E
∗ is globally asymptotically stable by constructing suitable Lyapunov function.Proposition 7.
Suppose that hypotheses (H1)–(H7) and R
¯
>
1 are satisfied. If R
1
>
1, then the infected equilibrium with CTL response E
∗ is globally asymptotically stable when τ
=
0.Proof.
By the following Lyapunov function,(48)
V
3
=
x
-
x
∗
-
∫
x
∗
x
f
x
∗
,
v
∗
f
s
,
v
∗
d
s
+
u
-
u
∗
-
u
∗
ln
⁡
u
u
∗
+
σ
+
μ
2
σ
y
-
y
∗
-
y
∗
ln
⁡
y
y
∗
+
f
x
∗
,
v
∗
μ
4
v
-
v
∗
-
v
∗
ln
⁡
v
v
∗
+
σ
+
μ
2
σ
p
q
-
m
z
-
z
∗
-
z
∗
ln
⁡
z
z
∗
,
V
3 is positive definite with respect to E
∗. For t
≥
0, the time derivative of V
3 along the solutions of system (2) is (49)
V
˙
3
=
1
-
f
x
∗
,
v
∗
f
x
,
v
∗
x
˙
+
1
-
u
∗
u
u
˙
+
σ
+
μ
2
σ
1
-
y
∗
y
y
˙
+
f
x
∗
,
v
∗
μ
4
1
-
v
∗
v
v
˙
+
σ
+
μ
2
σ
p
q
-
m
1
-
z
∗
z
z
˙
=
1
-
f
x
∗
,
v
∗
f
x
,
v
∗
λ
-
f
x
,
v
v
-
μ
1
x
+
1
-
u
∗
u
f
x
,
v
v
-
σ
+
μ
2
u
+
σ
+
μ
2
σ
1
-
y
∗
y
σ
u
-
p
y
z
-
μ
3
y
+
f
x
∗
,
v
∗
μ
4
1
-
v
∗
v
k
y
-
μ
4
v
+
σ
+
μ
2
σ
p
q
-
m
1
-
z
∗
z
q
y
z
-
μ
5
z
-
m
y
z
.Note that λ
=
f
(
x
∗
,
v
∗
)
v
∗
+
μ
1
x
∗, f
(
x
∗
,
v
∗
)
v
∗
=
(
σ
+
μ
2
)
u
∗, and μ
3
y
∗
=
σ
u
∗
-
p
y
∗
z
∗; we have (50)
V
˙
3
=
1
-
f
x
∗
,
v
∗
f
x
,
v
∗
f
x
∗
,
v
∗
v
∗
+
μ
1
x
∗
-
f
x
,
v
v
-
μ
1
x
+
1
-
u
∗
u
f
x
,
v
v
-
f
x
∗
,
v
∗
v
∗
u
∗
u
+
σ
+
μ
2
σ
1
-
y
∗
y
p
y
z
∗
-
p
y
z
+
f
x
∗
,
v
∗
v
∗
σ
u
∗
1
-
y
∗
y
σ
u
-
σ
u
∗
y
y
∗
+
f
x
∗
,
v
∗
μ
4
1
-
v
∗
v
k
y
-
μ
4
v
+
σ
+
μ
2
σ
·
p
q
-
m
1
-
z
∗
z
q
y
z
-
μ
5
z
-
m
y
z
=
μ
1
x
∗
-
x
1
-
f
x
∗
,
v
∗
f
x
,
v
∗
+
f
x
∗
,
v
∗
·
v
∗
-
1
+
f
x
,
v
f
x
,
v
∗
v
v
∗
-
v
v
∗
+
f
x
,
v
∗
f
x
,
v
+
f
x
∗
,
v
∗
·
v
∗
5
-
f
x
∗
,
v
∗
f
x
,
v
∗
-
f
x
,
v
f
x
∗
,
v
∗
v
v
∗
u
∗
u
-
y
∗
y
u
u
∗
-
f
x
,
v
∗
f
x
,
v
-
y
y
∗
v
∗
v
.Since the arithmetic mean is greater than or equal to the geometric mean, it has (51)
5
-
f
x
∗
,
v
∗
f
x
,
v
∗
-
f
x
,
v
f
x
∗
,
v
∗
v
v
∗
u
∗
u
-
y
∗
y
u
u
∗
-
f
x
,
v
∗
f
x
,
v
-
y
y
∗
v
∗
v
≤
0
.From hypotheses (H3)-
(
H
4
), we have (52)
-
1
+
f
x
,
v
f
x
,
v
∗
v
v
∗
-
v
v
∗
+
f
x
,
v
∗
f
x
,
v
=
f
x
,
v
-
f
x
,
v
∗
f
x
,
v
∗
f
x
,
v
v
-
f
x
,
v
∗
v
∗
f
x
,
v
v
∗
≤
0
.Therefore, V
˙
3
≤
0 if R
1
>
1. Then it follows from stability theorems in [34, 35] that the infected equilibrium CTL response E
∗ is stable for τ
=
0 if R
1
>
1. Similarly, by LaSalle invariance principle, we can show that the infected equilibrium CTL response E
∗ is also globally attractive for τ
=
0 if R
1
>
1.Next, we consider the case whenτ
>
0. Since α
5
>
0, s
=
0 is not a root of (40). We suppose (40) has a purely imaginary root s
=
i
ω
(
ω
>
0
) for some τ
>
0. Substituting s
=
i
ω into (40) and separating the real and imaginary parts, we have(53)
ω
5
-
A
2
ω
3
+
A
4
ω
=
B
1
ω
4
-
B
3
ω
2
+
B
5
sin
⁡
ω
τ
+
B
2
ω
3
-
B
4
ω
cos
⁡
ω
τ
,
A
1
ω
4
-
A
3
ω
2
+
A
5
=
-
B
1
ω
4
-
B
3
ω
2
+
B
5
cos
⁡
ω
τ
+
B
2
ω
3
-
B
4
ω
sin
⁡
ω
τ
.Squaring and adding the two equations of (53), it follows that(54)
ω
10
+
C
1
ω
8
+
C
2
ω
6
+
C
3
ω
4
+
C
4
ω
2
+
C
5
=
0
,where (55)
C
1
=
A
1
2
-
2
A
2
-
B
1
2
,
C
2
=
A
2
2
+
2
A
4
-
2
A
1
A
3
+
2
B
1
B
3
-
B
2
2
,
C
3
=
A
3
2
-
2
A
2
A
4
-
B
3
2
+
2
B
2
B
4
+
2
A
1
A
5
-
2
B
1
B
5
,
C
4
=
A
4
2
-
B
4
2
-
2
A
3
A
5
+
2
B
3
B
5
,
C
5
=
A
5
2
-
B
5
2
.Letting ν
=
ω
2, (54) can be written as(56)
h
ν
=
ν
5
+
C
1
ν
4
+
C
2
ν
3
+
C
3
ν
2
+
C
4
ν
+
C
5
=
0
.Then we have(57)
h
′
ν
=
5
ν
4
+
4
C
1
ν
3
+
3
C
2
ν
2
+
2
C
3
ν
+
C
4
.Denote (58)
p
1
=
-
6
25
C
1
2
+
3
5
C
2
,
q
1
=
8
125
C
1
3
-
6
25
C
1
C
2
+
2
5
C
3
,
r
1
=
-
3
625
C
1
4
+
3
125
C
1
2
C
2
-
2
25
C
1
C
3
+
1
5
C
4
,
Θ
0
=
p
1
2
-
4
r
1
,
p
2
=
-
1
3
p
1
2
-
4
r
1
,
q
2
=
-
2
27
p
1
3
+
8
3
p
1
r
1
-
q
1
2
,
Θ
1
=
1
27
p
2
3
+
1
4
q
2
2
,
s
∗
=
-
q
2
2
+
Θ
1
3
+
-
q
2
2
-
Θ
1
3
+
1
3
p
1
,
Θ
2
=
-
s
∗
-
p
1
+
2
q
1
s
∗
-
p
1
,
Θ
3
=
-
s
∗
-
p
1
-
2
q
1
s
∗
-
p
1
.By a similar argument as that in [38], we have the following results.Lemma 8.
For the polynomial equation (56), the following results hold.(i)
Equation (56) has at least one positive root, if one of the following conditions (
a
)
–
(
d
) holds:
(a)
C
5
<
0
.
(b)
C
5
≥
0
,
q
1
=
0
,
Θ
0
≥
0, and p
1
<
0 or r
1
≤
0 and there exists ν
∗
∈
{
ν
1
,
ν
2
,
ν
3
,
ν
4
} such that ν
∗
>
0 and h
(
ν
∗
)
≤
0, where ν
i
=
y
i
-
1
/
5
C
1
(
i
=
1,2
,
3,4
), and (59)
y
1
=
-
p
1
+
Θ
0
2
,
y
2
=
-
-
p
1
+
Θ
0
2
,
y
3
=
-
p
1
-
Θ
0
2
,
y
4
=
-
-
p
1
-
Θ
0
2
.
(c)
C
5
≥
0
,
q
1
≠
0
,
s
∗
>
p
1
,
Θ
2
≥
0, or Θ
3
≥
0 and there exists ν
∗
∈
{
ν
1
∗
,
ν
2
∗
,
ν
3
∗
,
ν
4
∗
} such that ν
∗
>
0 and h
(
ν
∗
)
≤
0, where ν
i
=
y
i
-
1
/
5
C
1
(
i
=
1,2
,
3,4
), and (60)
y
1
=
-
s
∗
-
p
1
+
Θ
2
2
,
y
2
=
-
s
∗
-
p
1
-
Θ
2
2
,
y
3
=
s
∗
-
p
1
+
Θ
3
2
,
y
4
=
s
∗
-
p
1
-
Θ
3
2
.
(d)
C
5
≥
0
,
q
≠
0
,
s
∗
<
p
1
,
q
1
2
/
4
(
p
1
-
s
∗
)
2
+
1
/
2
s
∗
=
0
,
ν
¯
>
0, and h
(
ν
¯
)
≤
0, where ν
¯
=
q
1
/
2
(
p
1
-
s
∗
)
-
1
/
5
C
1.
(ii)
If the conditions(a)–(d) of (
i
) are all not satisfied, then (56) has no positive real root.Suppose thath
(
ν
)
=
0 has positive real roots. Without loss of generality, we may assume that (56) has k
¯
(
1
≤
k
¯
≤
5
) positive real roots, denoted, respectively, as ν
1
,
ν
2
,
…
,
ν
k
¯. Then, (54) has positive real roots ω
k
¯
=
v
k
¯. From (40), we get (61)
cos
⁡
ω
τ
=
ω
5
-
A
2
ω
3
+
A
4
ω
B
2
ω
3
-
B
4
ω
-
A
1
ω
4
-
A
3
ω
2
+
A
5
B
1
ω
4
-
B
3
ω
2
+
B
5
B
2
ω
3
-
B
4
ω
2
+
B
1
ω
4
-
B
3
ω
2
+
B
5
2
≡
L
ω
.Therefore, let(62)
τ
k
j
=
1
ω
k
arccos
⁡
L
ω
k
+
2
π
j
,where k
=
1,2
,
…
,
k
¯
,
j
=
0,1
,
…. Then ±
i
ω
k are a pair of purely imaginary roots of (54) with τ
=
τ
k
(
j
).Define(63)
τ
0
=
τ
k
0
0
=
min
k
∈
1,2
,
…
,
k
¯
⁡
τ
k
0
,
ω
0
=
ω
k
0
.Let s
(
τ
)
=
ξ
(
τ
)
+
i
ω
(
τ
) be a root of (40) satisfying ξ
(
τ
k
(
j
)
)
=
0 and ω
(
τ
k
(
j
)
)
=
ω
k
. Differentiating the two sides of (40) with respect to τ and noticing that s is a function of τ, it follows that (64)
d
s
d
τ
-
1
=
-
5
s
4
+
4
A
1
s
3
+
3
A
2
s
2
+
2
A
3
s
+
A
4
s
s
5
+
A
1
s
4
+
A
2
s
3
+
A
3
s
2
+
A
4
s
+
A
5
+
4
B
1
s
3
+
3
B
2
s
2
+
2
B
3
s
+
B
4
s
B
1
s
4
+
B
2
s
3
+
B
3
s
2
+
B
4
s
+
B
5
-
τ
s
.Thus, we get (65)
d
Res
τ
d
τ
τ
=
τ
k
j
-
1
=
-
5
ω
k
4
-
3
A
2
ω
k
2
+
A
4
-
ω
k
6
+
A
2
ω
k
4
-
A
4
ω
k
2
-
4
A
1
ω
k
3
-
2
A
3
ω
k
A
1
ω
k
5
-
A
3
ω
k
3
+
A
5
ω
k
-
ω
k
6
+
A
2
ω
k
4
-
A
4
ω
k
2
2
+
A
1
ω
k
5
-
A
3
ω
k
3
+
A
5
ω
k
2
+
-
3
B
2
ω
k
2
+
B
4
B
2
ω
k
4
-
B
4
ω
k
2
+
-
4
B
1
ω
k
3
+
2
B
3
ω
k
B
1
ω
k
5
-
B
3
ω
k
3
+
B
5
ω
k
B
2
ω
k
4
-
B
4
ω
k
2
2
+
B
1
ω
k
5
-
B
3
ω
k
3
+
B
5
ω
k
2
.From (40), we attain (66)
ω
5
-
A
2
ω
3
+
A
4
ω
2
+
A
1
ω
4
-
A
3
ω
2
+
A
5
2
=
B
2
ω
3
-
B
4
ω
2
+
B
1
ω
4
-
B
3
ω
2
+
B
5
2
.Then(67)
d
Res
τ
d
τ
τ
=
τ
k
j
-
1
=
5
ν
k
4
+
4
C
1
ν
k
3
+
3
C
2
ν
k
2
+
2
C
3
ν
k
+
C
4
B
1
ω
k
4
-
B
3
ω
k
2
+
B
5
2
+
B
2
ω
k
2
-
B
4
2
ω
k
2
=
h
′
ν
k
B
1
ω
k
4
-
B
3
ω
k
2
+
B
5
2
+
B
2
ω
k
2
-
B
4
2
ω
k
2
.Therefore, it follows that (68)
sign
⁡
d
Res
τ
d
τ
τ
=
τ
k
j
=
sign
⁡
d
Res
τ
d
τ
τ
=
τ
k
j
-
1
=
sign
⁡
h
′
ν
k
.Since ν
k
>
0, we can know that R
e
d
s
k
(
τ
)
/
d
τ
|
τ
=
τ
k
(
j
) and h
′
ν
k have the same sign.From the above analysis, we have the following results.Theorem 9.
Letτ
k
(
j
)
,
τ
0, and ω
0 be defined by (62) and (63). If R
¯
>
1 and R
1
>
1 are satisfied, then the following results hold:(i)
If the conditions(
a
)
–
(
d
) of Lemma 8 are all not satisfied, then the infected equilibrium with CTL response E
∗ is locally asymptotically stable for all time delay τ
>
0
.
(ii)
If one of the conditions(
a
)
–
(
d
) of Lemma 8 is satisfied, then the infected equilibrium with CTL response E
∗ is locally asymptotically stable for τ
∈
[
0
,
τ
0
) and unstable for τ
>
τ
0.
(iii)
If all the conditions as stated in (ii) hold andh
′
(
ν
k
)
≠
0, then system (2) undergoes a Hopf bifurcation at E
∗ when τ
=
τ
k
j
(
j
=
0,1
,
2
,
…
).


## 5. Conclusion and Numerical Simulations

In this paper, we proposed a class of delayed HIV infection model (2) with general incidence rate and immune impairment. This general incidence only satisfies some general hypotheses and includes many types of special incidence functions as special cases. First, we discussed the nonnegativity and boundedness of the solutions and the existence of equilibria of system (2). Then, by constructing suitable Lyapunov functionals and using Lyapunov-LaSalle invariance principle and Hopf bifurcation theorem, we proved the following results.IfR
0
≤
1, the infection-free equilibrium E
0 is globally asymptotically stable for any time delay τ
≥
0; that is, any solution (
x
(
t
)
,
u
(
t
)
,
y
(
t
)
,
v
(
t
)
,
z
(
t
)
)
→
E
0
=
(
x
0
,
0,0
,
0,0
). In biology, this means that the virus can be finally cleared from the body and the disease dies out. At the same time, as the time t increases, the numbers of latently infected cells, actively infected cells, and CTLs trends to zero and the number of uninfected cells trends to a constant x
0.IfR
0
>
1
≥
R
1 and R
¯
>
1, the infected equilibrium without immunity E
1 is globally asymptotically stable for any time delay τ
≥
0; that is, any solution (
x
(
t
)
,
u
(
t
)
,
y
(
t
)
,
v
(
t
)
,
z
(
t
)
)
→
E
1
=
(
x
1
,
u
1
,
y
1
,
v
1
,
0
). In biology, this indicates that the HIV infection will finally become chronic with no persistent CTL response.IfR
1
>
1 and R
¯
>
1, there exists a unique infected equilibrium with CTL response E
∗. The result of Theorem 9 implies that the time delay τ can destabilize the stability of the infected equilibrium with CTL response E
∗ and leads to the occurrence of Hopf bifurcations.If the time delayτ
∈
[
0
,
τ
0
), the infected equilibrium with CTL response E
∗ is locally asymptotically stable. In biology, this implies that the HIV infection may become chronic and the CTL immune response may be persistent. When the time delay τ passes through the critical value τ
0, the infected equilibrium with CTL response E
∗ will become unstable and a Hopf bifurcation occurs under some conditions. In biology, this suggests that as the time delay τ increases, the numbers of the uninfected cells, latently infected cells, actively infected cells, free virus, and CTLs will first attend constant values and then become oscillated.We now give numerical simulations to illustrate the main results in Sections3 and 4.Let us choosef
(
x
,
v
)
v
=
β
x
v. Then we have that R
0
=
k
λ
β
σ
/
(
μ
1
μ
3
μ
4
(
σ
+
μ
2
)
) and R
1
=
R
0
-
k
β
μ
5
/
(
μ
1
μ
4
(
q
-
m
)
). Based on the numerical simulations in [12, 19, 21, 22, 27], let us take the following data: (69)
λ
=
270
,
β
=
0.001
,
k
=
6
,
m
=
0.001
,
p
=
0.04
,
q
=
0.025
,
σ
=
0.001
,
μ
1
=
0.02
,
μ
2
=
0.1
,
μ
3
=
0.8
,
μ
4
=
1.2
,
μ
5
=
0.05
.Direct calculations show that R
0
=
0.8354
<
1 and R
1
=
0.3146
<
1; system (2) has the infection-free equilibrium E
0
=
(
13500,0
,
0,0
,
0
). By Theorem 4, the infection-free equilibrium E
0 is globally asymptotically stable for any time delay τ
≥
0. Figure 1 gives the phase trajectories of system (2) with suitable initial condition.Figure 1
Phase trajectories of system (2) with R
0
≤
1.Next, let us choose the following data:(70)
λ
=
270
,
β
=
0.001
,
k
=
6
,
m
=
0.01
,
p
=
0.04
,
q
=
0.025
,
σ
=
0.002
,
μ
1
=
0.02
,
μ
2
=
0.1
,
μ
3
=
0.8
,
μ
4
=
1.2
,
μ
5
=
0.05
.Direct computations show that R
0
=
1.6544
>
1, R
¯
=
1.9853
>
1, and R
1
=
0.8211
<
1; system (2) has the infected equilibrium without immunity E
1
=
(
8160,1047.0588,2.6176,13.0882,0
). Therefore, by Theorem 5, the infected equilibrium without immunity E
1 is globally asymptotically stable for any time delay τ
≥
0. Figure 2 gives the phase trajectories of system (2) with suitable initial condition.Figure 2
Phase trajectories of system (2) with R
0
>
1
≥
R
1.Furthermore, let us choose the following data:(71)
λ
=
270
,
β
=
0.001
,
k
=
6
,
m
=
0.012
,
p
=
0.04
,
q
=
0.025
,
σ
=
0.004
,
μ
1
=
0.02
,
μ
2
=
0.1
,
μ
3
=
0.8
,
μ
4
=
1.2
,
μ
5
=
0.05
.Then we have that R
0
=
3.2452
>
1, R
¯
=
3.375
>
1, and R
1
=
2.2837
>
1 and (56) has no positive root. System (2) has the infected equilibrium with CTL response E
∗
=
(
6882.3529,1272.6244,3.84615,19.2308,13.0882
). From Theorem 9(i), the infected equilibrium with CTL response E
∗ is locally asymptotically stable for any time delay τ
>
0. Figure 3 gives the phase trajectories of system (2) with suitable initial condition.Figure 3
Phase trajectories of system (2) with R
1
>
1
,
τ
=
28. The initial condition is (7970, 1000, 2.6, 0.2, 17).Finally, let us choose the following data:(72)
λ
=
270
,
β
=
0.001
,
k
=
6
,
m
=
0.01
,
p
=
0.04
,
q
=
0.025
,
σ
=
0.004
,
μ
1
=
0.02
,
μ
2
=
0.1
,
μ
3
=
0.8
,
μ
4
=
1.2
,
μ
5
=
0.05
.Then we have that R
0
=
3.2452
>
1, R
¯
=
3.8942
>
1, and R
1
=
2.4119
>
1, (56) has two positive roots, and h
′
(
ν
k
)
≠
0. By simple computations, we have ω
0
≈
0.0394 and τ
0
≈
27.2546. From Theorem 9(ii), the infected equilibrium with CTL response E
∗
=
(
7363.6364,1180.0699,3.3333,16.6667,15.4021
) is asymptotically stable if 0
<
τ
<
τ
0 and unstable if τ
>
τ
0. Figure 4 gives the phase trajectories of system (2) with τ
<
τ
0 and suitable initial condition. Figure 5 gives the phase trajectories of system (2) with τ
>
τ
0 and suitable initial condition and shows the occurrence of the Hopf bifurcations.Figure 4
Phase trajectories of system (2) with R
1
>
1
,
τ
=
25
<
τ
0. The initial condition is (7970, 1000, 2.6, 0.2, 17).Figure 5
Phase trajectories of system (2) with R
1
>
1
,
τ
=
28
>
τ
0. The initial condition is (7970, 1000, 2.6, 0.2, 17).Sincev
∗
=
k
μ
5
/
(
μ
4
(
q
-
m
)
) and z
∗
=
μ
1
μ
3
μ
4
(
q
-
m
)
(
R
1
-
1
)
/
(
p
μ
1
μ
4
(
q
-
m
)
+
p
β
k
μ
5
), it is easy to see that the number of free viruses is increased and the number of CTLs is decreased with respect to the immune impairment rate m. For example, if we choose m
=
0.015, then v
1
∗
=
25 and z
1
∗
=
8.8462. If we choose m
=
0.016, then v
2
∗
=
27.7778 and z
2
∗
=
7.1691. If we choose m
=
0.017, then v
3
∗
=
31.25 and z
3
∗
=
5.3283. Figure 6 shows that, for any initial conditions, when m
=
0.015,0.016,0.017, the numbers of free viruses trend to v
1
∗, v
2
∗, and v
3
∗, respectively. Figure 7 shows that, for any initial conditions, when m
=
0.015,0.016,0.017, the numbers of CTLs trend to z
1
∗, z
2
∗, and z
3
∗, respectively. In Figures 6 and 7, all the data are chosen as in Figure 4 except the immune impairment rate m.Figure 6
The curves of the free virus of system (2) with τ
=
20, m
=
0.015, 0.016, and 0.017. The initial conditions are chosen as (
7000,1000,8
,
35,10
), (
7000,1000,8
,
20,10
), and (
7000,1000,8
,
8,10
), respectively.Figure 7
The curves of the CTLs of system (2) with τ
=
20, m
=
0.015, 0.016, and 0.017. The initial conditions are chosen as (
7000,1000,8
,
2,8
), (
7000,1000,8
,
2,13
), and (
7000,1000,8
,
2,18
), respectively.As immune impairment ratem increases, the CTL response gradually becomes weak and the individuals eventually develop AIDS. Thus, in order to control the HIV infection, we should decrease the value of m. Numerical simulations show the similar known results (see, e.g., [22]).



---

*Source: 206205-2015-08-04.xml*