# Stability of an HIV/AIDS Treatment Model with Different Stages

Authors: Hai-Feng Huo; Rui Chen
Journal: Discrete Dynamics in Nature and Society
 (2015)

---

## Abstract

An HIV/AIDS treatment model with different stages is proposed in this paper. The stage of the HIV infection is divided into two stages, that is, HIV-positive in the asymptomatic stage of HIV infection and HIV-positive individuals in the pre-AIDS stage. The fact that some individuals with HIV-positive individuals after the treatment can be transformed into the compartment of HIV-positive individuals in the asymptomatic stage of HIV infection, the compartment of HIV-positive individuals in the pre-AIDS stage, or the compartment of individuals with full-blown AIDS is also considered. Mathematical analyses establish the idea that the global dynamics of the HIV/AIDS model are determined by the basic reproduction numberR0. The disease-free equilibrium is globally asymptotically stable if R0<1. The endemic equilibrium is globally asymptotically stable if R0>1 for a special case. Numerical simulations are also conducted to support the analytic results.


---

## Body

## 1. Introduction

Human immunodeficiency virus (HIV) is one of the major life threatening viruses that are spreading worldwide. AIDS is caused by the human immunodeficiency virus (HIV), which has developed into a global pandemic since the first patient was identified in 1981, making it one of the most destructive epidemics in history [1, 2]. The HIV prevalence among pregnant women attending antenatal clinics rose from 0.8% in 1990 to 30.2% in 2005 and declined slightly to 29.1% in 2006 [3, 4].Because of the fact that there is no vaccine, there are many obstacles in the AIDS treatment. Recently, the most prevalent treatment strategy for HIV infected patients is highly active antiretroviral therapies (HAART), which can prolong the lifespans and improve their life quality of infected individuals [5]. Current antiretroviral HIV/AIDS therapy involves the simultaneous administration of two or more antiviral drugs, typically chosen from two major classes: reverse transcriptase inhibitors (RTIs) and protease inhibitors (PIs) [6]. With the help of the antiretroviral therapy, patients can prolong the lifespan during which they are free of HIV-related symptoms and therefore improve their life quality. At least 3.0 millions of people have been saved in the United States, which highlights the significant advances made in HIV disease treatment [7]. That is the reason why the early treatment of AIDS is very meaningful and of great urgency.It is well known that mathematical modeling is a very important approach to understand the dynamics of any epidemic and to further develop various control and prevention policies [8, 9]. For the sake of the importance in biological meanings, dynamical properties of HIV/AIDS models and general theory on such dynamical systems have been studied by many authors in recent years [10–13] and references cited therein. They propose mathematical models that described the interactions of the HIV virus, CD4 cells, and CTLs within host, and these mathematical models can provide some insights into the dynamics of HIV viral load in vivo and may play a significant role in the development of a better understanding of HIV/AIDs and drug therapies.Huo and Chen [14] introduce an HIV/AIDS epidemic model with antiretroviral treatment and establish the idea that the global dynamics of the HIV/AIDS model are determined by the basic reproduction number R0. They show that the early treatment of AIDS is necessary and meaningful. But they do not consider the compartments of HIV-positive individuals in the asymptomatic stage of HIV infection and the compartment of HIV-positive individuals in the pre-AIDS stage not receiving antiretroviral treatment. Yusuf and Benyah [15] present a deterministic model for controlling the spread of the disease using change in sexual habits and antiretroviral (ARV) therapy as control measures and formulate a fixed time optimal control problem subject to the model dynamics with the goal of finding the optimal combination of the two control measures that will minimize the cost of the control efforts as well as the incidence of the disease. But they neglect the fact that some individuals with HIV-positive after the treatment can be transformed into the compartment of HIV-positive individuals in the asymptomatic stage of HIV infection, the compartment of HIV-positive individuals in the pre-AIDS stage, or the compartment of individuals with full-blown AIDS.Motivated by the above, in this paper, we will divide the stage of the HIV infection into two stages, that is, HIV-positive individuals in the asymptomatic stage of HIV infection and HIV-positive individuals in the pre-AIDS stage. We also consider the fact that some individuals with HIV-positive after the treatment can be transformed into the compartment of HIV-positive individuals in the asymptomatic stage of HIV infection, the compartment of HIV-positive individuals in the pre-AIDS stage, or the compartment of individuals with full-blown AIDS.The organization of this paper is as follows. In the next section, we propose an HIV/AIDS model. In Section3, the existence and the global stability of equilibria are investigated. Some numerical simulations are given in Section 4. Some discussions are given in last section.


## 2. The Model

### 2.1. System Description

The total population is divided into six compartments:S(t), H(t), P(t), A(t), T(t), and R(t). S(t) represents the number of susceptible individuals; H(t) represents the number of HIV-positive individuals in the asymptomatic stage of HIV infection; P(t) represents the number of HIV-positive individuals in the pre-AIDS stage not receiving treatment; A(t) represents the number of individuals with full-blown AIDS not receiving treatment; T(t) represents the number of individuals who are receiving ARV treatment; R(t) represents the number of individuals who change their sexual habits and maintain the habits for the rest of their lives. The total number of population at time t is given by
(1)N(t)=S(t)+H(t)+P(t)+A(t)+T(t)+R(t).
The model structure is shown in Figure 1. The transfer diagram leads to the following system of ordinary differential equations:
(2)S˙=Λ-β1SH-β2SP-dS-μ1S,H˙=β1SH+β2SP-k1H-k3H-dH+(α1+α3)T,P˙=k1H-k4P-k2P-dP+(α2+α4)T,A˙=k2P+α5T-(d+δ1)A,T˙=k3H+k4P-(α1+α3)T-(α2+α4)T-α5T-(d+δ2)T,R˙=μ1S-dR.
For simplicity, denote k1+k3+d=a, k2+k4+d=b, α1+α2+α3+α4+α5+d+δ2=m, α1+α3=h, α2+α4=g, and d+δ1=n, and then system 2 becomes
(3)S˙=Λ-β1SH-β2SP-dS-μ1S,H˙=β1SH+β2SP-aH+hT,P˙=k1H-bP+gT,A˙=k2P+α5T-nA,T˙=k3H+k4P-mT,R˙=μ1S-dR,
where Λ is recruitment rate, d is the natural death rate, β1 is the contact rate between the susceptible and the infected individuals. β2 is the contact rate between the susceptible and the pre-AIDS stages. k1 is the rate at which individuals leave the infection class and become individuals in the pre-AIDS stage. k2 is the rate at which individuals leave the pre-AIDS class and become individuals with full-blown AIDS. k3 is the rate at which an individual in the asymptomatic stage of HIV infection receives treatment. k4 is the rate at which an individual in the pre-AIDS stage receives treatment. δ1 and δ2 are the disease-induced death rates for individuals in compartments A(t) and T(t). μ1 is the rate at which susceptible individuals change their sexual habits per unit time. α1 is the rate at which a treated individual leaves compartment T(t). α2 is the rate at which treatment of an individual with HIV fails. α3 is the rate at which an individual in the pre-AIDS stage is treated successfully. α4 is the rate at which tratment of an individual in the pre-AIDS stage fails. α5 is the rate at which treatment of an individual in the pre-AIDS stage fails. All the parameters can be found in Table 1.Table 1
Description and estimation of parameters.
Parameter
Description
Estimated value
Date source
Λ
Recruitment rate of the population
0.55 year−1
Estimate
β
1
Transmission coefficient for contact with the H class
0.3 year−1
Estimate
β
2
Transmission coefficient for contact with the P class
0.2 year−1
Estimate
d
Natural death rate
0.0196 year−1
[15]
k
1
Progression rate from the H class into the P class
0.498 year−1
Estimate
k
2
Progression rate from the P class into the A class
0.08 year−1
Estimate
k
3
Proportion of the H class receiving treatment
Variable
Estimate
k
4
Proportion of the P class receiving treatment
Variable
Estimate
α
1
The proportion of the H class successful treatment
0.2 year−1
Estimate
α
2
The proportion of the P class treatment failure
0.05 year−1
Estimate
α
3
The proportion of the P class successful treatment
0.15 year−1
Estimate
α
4
The proportion of the P class treatment failure
0.03 year−1
Estimate
α
5
The proportion of the P class treatment failure and seriousness
0.01 year−1
Estimate
δ
1
The disease-related death rate of the AIDS
0.0909 year−1
[15]
δ
2
The disease-related death rate of being treated
0.0667 year−1
[15]
μ
1
The rate of susceptible individuals who changed their habits
0.03 year−1
[15]Figure 1
Transfer diagram of the model2.


### 2.2. Basic Properties

#### 2.2.1. Invariant Region

It is necessary to prove that all solutions of system2 with positive initial data will remain positive for all times t>0. This will be established by the following lemma.Lemma 1.
The feasible regionΩ defined by
(4)Ω=St,Ht,Pt,At,Tt,RtΛdhhhh∈R+6∣S(t)+H(t)hhhh+P(t)+A(t)+T(t)+R(t)≤Λd,
with initial condition S(0)≥0, H(0)≥0, P(0)≥0, A(0)≥0, T(0)≥0, and R(0)≥0, is positive invariant for system 2.Proof.
Adding the equations of system2, we obtain
(5)dNdt=Λ-dN-δ1A-δ2T≤Λ-dN.
It follows that
(6)0≤N(t)≤Λd+N(0)e-dt,
where N(0) represents the initial values of the total population. Thus lim⁡t→∞sup⁡N(t)≤Λ/d. It implies that the region Ω={(S(t),H(t),P(t),A(t),T(t),R(t))∈R+6∣S(t)+H(t)+P(t)+A(t)+T(t)+R(t)≤Λ/d} is a positively invariant set for system 2.


#### 2.2.2. Positivity of Solutions

Lemma 2.
IfS(0)≥0, H(0)≥0, P(0)≥0, A(0)≥0, R(0)≥0, and T(0)≥0, the solutions S(t), H(t), P(t), A(t), T(t), and R(t) of system 3 are positive for all t>0.Proof.
Under the given initial conditions, we need to prove that the solutions of system3 are positive [16]; if not, we assume that there is a contradiction: there exists a first time t1 such that
(7)St1=0,S′t1<0,Ht≥0,Pt≥0,At≥0,Tt≥0,Rt≥0,0<t<t1,
there exists a t2 such that
(8)Ht2=0,H′t2<0,St≥0,Pt≥0,At≥0,Tt≥0,Rt≥0,0<t<t2,
there exists a t3 such that
(9)Pt3=0,P′t3<0,St≥0,Ht≥0,At≥0,Tt≥0,Rt≥0,0<t<t3,
there exists a t4 such that
(10)At4=0,A′t4<0,St≥0,Ht≥0,Pt≥0,Tt≥0,Rt≥0,0<t<t4.
there exists a t5 such that
(11)Tt5=0,T′t5<0,St≥0,Ht≥0,Pt≥0,At≥0,Rt≥0,0<t<t5,
and there exists a t6 such that
(12)Rt6=0,R′t6<0,St≥0,Ht≥0,Pt≥0,At≥0,Tt≥0,0<t<t6.
In the first case, we have
(13)S′(t1)=Λ>0,
which is a contradiction, so that we can derive that S(t)≥0, t≥0. In the second case, we have
(14)H′(t2)=(α1+α3)T(t2)+β2S(t2)P(t2)≥0,
which is a contradiction, so that we can derive that H(t)≥0, t≥0. Similarly, it can be shown that P(t)≥0, A(t)≥0, T(t)≥0, and R(t)≥0 for all t≥0.
Thus, the solutionsS(t), H(t), P(t), A(t), T(t), and R(t) of system 3 remain positive for all t>0.


## 2.1. System Description

The total population is divided into six compartments:S(t), H(t), P(t), A(t), T(t), and R(t). S(t) represents the number of susceptible individuals; H(t) represents the number of HIV-positive individuals in the asymptomatic stage of HIV infection; P(t) represents the number of HIV-positive individuals in the pre-AIDS stage not receiving treatment; A(t) represents the number of individuals with full-blown AIDS not receiving treatment; T(t) represents the number of individuals who are receiving ARV treatment; R(t) represents the number of individuals who change their sexual habits and maintain the habits for the rest of their lives. The total number of population at time t is given by
(1)N(t)=S(t)+H(t)+P(t)+A(t)+T(t)+R(t).
The model structure is shown in Figure 1. The transfer diagram leads to the following system of ordinary differential equations:
(2)S˙=Λ-β1SH-β2SP-dS-μ1S,H˙=β1SH+β2SP-k1H-k3H-dH+(α1+α3)T,P˙=k1H-k4P-k2P-dP+(α2+α4)T,A˙=k2P+α5T-(d+δ1)A,T˙=k3H+k4P-(α1+α3)T-(α2+α4)T-α5T-(d+δ2)T,R˙=μ1S-dR.
For simplicity, denote k1+k3+d=a, k2+k4+d=b, α1+α2+α3+α4+α5+d+δ2=m, α1+α3=h, α2+α4=g, and d+δ1=n, and then system 2 becomes
(3)S˙=Λ-β1SH-β2SP-dS-μ1S,H˙=β1SH+β2SP-aH+hT,P˙=k1H-bP+gT,A˙=k2P+α5T-nA,T˙=k3H+k4P-mT,R˙=μ1S-dR,
where Λ is recruitment rate, d is the natural death rate, β1 is the contact rate between the susceptible and the infected individuals. β2 is the contact rate between the susceptible and the pre-AIDS stages. k1 is the rate at which individuals leave the infection class and become individuals in the pre-AIDS stage. k2 is the rate at which individuals leave the pre-AIDS class and become individuals with full-blown AIDS. k3 is the rate at which an individual in the asymptomatic stage of HIV infection receives treatment. k4 is the rate at which an individual in the pre-AIDS stage receives treatment. δ1 and δ2 are the disease-induced death rates for individuals in compartments A(t) and T(t). μ1 is the rate at which susceptible individuals change their sexual habits per unit time. α1 is the rate at which a treated individual leaves compartment T(t). α2 is the rate at which treatment of an individual with HIV fails. α3 is the rate at which an individual in the pre-AIDS stage is treated successfully. α4 is the rate at which tratment of an individual in the pre-AIDS stage fails. α5 is the rate at which treatment of an individual in the pre-AIDS stage fails. All the parameters can be found in Table 1.Table 1
Description and estimation of parameters.
Parameter
Description
Estimated value
Date source
Λ
Recruitment rate of the population
0.55 year−1
Estimate
β
1
Transmission coefficient for contact with the H class
0.3 year−1
Estimate
β
2
Transmission coefficient for contact with the P class
0.2 year−1
Estimate
d
Natural death rate
0.0196 year−1
[15]
k
1
Progression rate from the H class into the P class
0.498 year−1
Estimate
k
2
Progression rate from the P class into the A class
0.08 year−1
Estimate
k
3
Proportion of the H class receiving treatment
Variable
Estimate
k
4
Proportion of the P class receiving treatment
Variable
Estimate
α
1
The proportion of the H class successful treatment
0.2 year−1
Estimate
α
2
The proportion of the P class treatment failure
0.05 year−1
Estimate
α
3
The proportion of the P class successful treatment
0.15 year−1
Estimate
α
4
The proportion of the P class treatment failure
0.03 year−1
Estimate
α
5
The proportion of the P class treatment failure and seriousness
0.01 year−1
Estimate
δ
1
The disease-related death rate of the AIDS
0.0909 year−1
[15]
δ
2
The disease-related death rate of being treated
0.0667 year−1
[15]
μ
1
The rate of susceptible individuals who changed their habits
0.03 year−1
[15]Figure 1
Transfer diagram of the model2.


## 2.2. Basic Properties

### 2.2.1. Invariant Region

It is necessary to prove that all solutions of system2 with positive initial data will remain positive for all times t>0. This will be established by the following lemma.Lemma 1.
The feasible regionΩ defined by
(4)Ω=St,Ht,Pt,At,Tt,RtΛdhhhh∈R+6∣S(t)+H(t)hhhh+P(t)+A(t)+T(t)+R(t)≤Λd,
with initial condition S(0)≥0, H(0)≥0, P(0)≥0, A(0)≥0, T(0)≥0, and R(0)≥0, is positive invariant for system 2.Proof.
Adding the equations of system2, we obtain
(5)dNdt=Λ-dN-δ1A-δ2T≤Λ-dN.
It follows that
(6)0≤N(t)≤Λd+N(0)e-dt,
where N(0) represents the initial values of the total population. Thus lim⁡t→∞sup⁡N(t)≤Λ/d. It implies that the region Ω={(S(t),H(t),P(t),A(t),T(t),R(t))∈R+6∣S(t)+H(t)+P(t)+A(t)+T(t)+R(t)≤Λ/d} is a positively invariant set for system 2.


### 2.2.2. Positivity of Solutions

Lemma 2.
IfS(0)≥0, H(0)≥0, P(0)≥0, A(0)≥0, R(0)≥0, and T(0)≥0, the solutions S(t), H(t), P(t), A(t), T(t), and R(t) of system 3 are positive for all t>0.Proof.
Under the given initial conditions, we need to prove that the solutions of system3 are positive [16]; if not, we assume that there is a contradiction: there exists a first time t1 such that
(7)St1=0,S′t1<0,Ht≥0,Pt≥0,At≥0,Tt≥0,Rt≥0,0<t<t1,
there exists a t2 such that
(8)Ht2=0,H′t2<0,St≥0,Pt≥0,At≥0,Tt≥0,Rt≥0,0<t<t2,
there exists a t3 such that
(9)Pt3=0,P′t3<0,St≥0,Ht≥0,At≥0,Tt≥0,Rt≥0,0<t<t3,
there exists a t4 such that
(10)At4=0,A′t4<0,St≥0,Ht≥0,Pt≥0,Tt≥0,Rt≥0,0<t<t4.
there exists a t5 such that
(11)Tt5=0,T′t5<0,St≥0,Ht≥0,Pt≥0,At≥0,Rt≥0,0<t<t5,
and there exists a t6 such that
(12)Rt6=0,R′t6<0,St≥0,Ht≥0,Pt≥0,At≥0,Tt≥0,0<t<t6.
In the first case, we have
(13)S′(t1)=Λ>0,
which is a contradiction, so that we can derive that S(t)≥0, t≥0. In the second case, we have
(14)H′(t2)=(α1+α3)T(t2)+β2S(t2)P(t2)≥0,
which is a contradiction, so that we can derive that H(t)≥0, t≥0. Similarly, it can be shown that P(t)≥0, A(t)≥0, T(t)≥0, and R(t)≥0 for all t≥0.
Thus, the solutionsS(t), H(t), P(t), A(t), T(t), and R(t) of system 3 remain positive for all t>0.


## 2.2.1. Invariant Region

It is necessary to prove that all solutions of system2 with positive initial data will remain positive for all times t>0. This will be established by the following lemma.Lemma 1.
The feasible regionΩ defined by
(4)Ω=St,Ht,Pt,At,Tt,RtΛdhhhh∈R+6∣S(t)+H(t)hhhh+P(t)+A(t)+T(t)+R(t)≤Λd,
with initial condition S(0)≥0, H(0)≥0, P(0)≥0, A(0)≥0, T(0)≥0, and R(0)≥0, is positive invariant for system 2.Proof.
Adding the equations of system2, we obtain
(5)dNdt=Λ-dN-δ1A-δ2T≤Λ-dN.
It follows that
(6)0≤N(t)≤Λd+N(0)e-dt,
where N(0) represents the initial values of the total population. Thus lim⁡t→∞sup⁡N(t)≤Λ/d. It implies that the region Ω={(S(t),H(t),P(t),A(t),T(t),R(t))∈R+6∣S(t)+H(t)+P(t)+A(t)+T(t)+R(t)≤Λ/d} is a positively invariant set for system 2.


## 2.2.2. Positivity of Solutions

Lemma 2.
IfS(0)≥0, H(0)≥0, P(0)≥0, A(0)≥0, R(0)≥0, and T(0)≥0, the solutions S(t), H(t), P(t), A(t), T(t), and R(t) of system 3 are positive for all t>0.Proof.
Under the given initial conditions, we need to prove that the solutions of system3 are positive [16]; if not, we assume that there is a contradiction: there exists a first time t1 such that
(7)St1=0,S′t1<0,Ht≥0,Pt≥0,At≥0,Tt≥0,Rt≥0,0<t<t1,
there exists a t2 such that
(8)Ht2=0,H′t2<0,St≥0,Pt≥0,At≥0,Tt≥0,Rt≥0,0<t<t2,
there exists a t3 such that
(9)Pt3=0,P′t3<0,St≥0,Ht≥0,At≥0,Tt≥0,Rt≥0,0<t<t3,
there exists a t4 such that
(10)At4=0,A′t4<0,St≥0,Ht≥0,Pt≥0,Tt≥0,Rt≥0,0<t<t4.
there exists a t5 such that
(11)Tt5=0,T′t5<0,St≥0,Ht≥0,Pt≥0,At≥0,Rt≥0,0<t<t5,
and there exists a t6 such that
(12)Rt6=0,R′t6<0,St≥0,Ht≥0,Pt≥0,At≥0,Tt≥0,0<t<t6.
In the first case, we have
(13)S′(t1)=Λ>0,
which is a contradiction, so that we can derive that S(t)≥0, t≥0. In the second case, we have
(14)H′(t2)=(α1+α3)T(t2)+β2S(t2)P(t2)≥0,
which is a contradiction, so that we can derive that H(t)≥0, t≥0. Similarly, it can be shown that P(t)≥0, A(t)≥0, T(t)≥0, and R(t)≥0 for all t≥0.
Thus, the solutionsS(t), H(t), P(t), A(t), T(t), and R(t) of system 3 remain positive for all t>0.


## 3. Analysis of the Model

There exist one disease-free equilibriumE0 and one endemic equilibrium E* for system 3.


### 3.1. Disease-Free Equilibrium and the Reproduction Number

The model has a disease-free equilibrium given by(15)E0=Λμ1+d,0,0,0,0,Λμ1dμ1+d.
In the following, the basic reproduction number of system 3 will be obtained by the next generation matrix method formulated in [17, 18].Lettingx=(H,P,A,T,S,R)T, then system 3 can be written as
(16)dxdt=F(x)-V(x),
where
(17)F(x)=β1SH+β2SP00000,V(x)=aH-hTbP-gT-k1HnA-k2P-α5TmT-k3H-k4Pβ1SH+β2SP+dS+μ1S-ΛdR-μ1S.
The Jacobian matrices of F(x) and V(x) at the free equilibrium E0 are, respectively,
(18)DF(E0)=F4×400000000,DV(E0)=V4×404×104×1β1Λd+μ1β2Λd+μ100d+μ1-μ10000d,
where
(19)F=β1Λd+μ1β2Λd+μ100000000000000,V=a00-h-k1b0-g0-k2n-α5-k3-k40m.
The model reproduction number, denoted by R0, is thus given by
(20)R0=ρ(FV-1)=β1Λ(bm-k4g)+β2Λ(k1m+k3g)(d+μ1)(abm-k1k4h-k3hb-k4ga).


### 3.2. Global Stability ofE0

Theorem 3.
The disease-free equilibriumE0 is globally asymptotically stable for R0<1 and unstable otherwise.Proof.
According to [19, 20], we introduce the Lyapunov candidate function
(21)V(t)=H+ΘP+cT,Θ>0,c>0.
Differentiating V(t) with respect to time yields
(22)V˙=H˙+ΘP˙+cT˙.
Substituting the model equations into 22, we get
(23)V˙=β1SH+β2SP-aH+hT+Θ(k1H-bP+gT)+ck3H+k4P-mT.
With S0=Λ/(d+μ1) and R¯0=Λμ1/d(d+μ1), we have
(24)V˙≤β1ΛHd+μ1+β2ΛPd+μ1-aH+hT+Θk1H-bP+gT+ck3H+k4P-mT=β1Λd+μ1-a+k1Θ+ck3H+β2Λd+μ1-bΘ+ck4P+h+Θg-cmT.
In order to get the value of c and Θ, we can use
(25)β2Λd+μ1-bΘ+ck4=0,h+Θg-cm=0.
From 25, we can calculate
(26)c=β2Λg+hb(d+μ1)(bm-k4g)(d+μ1),Θ=β2Λm+k4h(d+μ1)(bm-k4g)(d+μ1).
Taking c and Θ into 24, we can get
(27)V˙≤β1Λd+μ1-a+k1β2Λm+k1k4h(d+μ1)(bm-k4g)(d+μ1)hhh+β2Λgk3+hbk3d+μ1bm-k4gd+μ1H=β1Λd+μ1(bm-k4g)+β2Λd+μ1(mk1+gk3)hhhhh-abm+ak4g+k1k4h+k3hbβ2Λd+μ1hhh·bm-k4g-1β1Λd+μ1H.
Since
(28)R0=ρ(FV-1)=β1Λ(bm-k4g)+β2Λk1m+k3gd+μ1abm-k1k4h-k3hb-k4ga,
we can obtain
(29)R0-1=β1Λd+μ1(bm-k4g)+β2Λd+μ1(mk1+gk3)hhhh-abm+ak4g+k1k4h+k3hbβ1Λd+μ1·abm-k1k4h-k3hb-k4ga-1,
when R0≤1,
(30)V˙≤abm-k1k4h-k3hb-k4gaR0-1bm-k4gH≤0.
WhileV˙=0, only H=T=P=0. Substituting H=T=P=0 into system 3, we can have S→Λ/(d+μ1), A→0, and R→Λμ1/d(d+μ1) as t→∞. According to LaSalle’s invariance principle [21], the disease-free equilibrium is globally asymptotically stable.


### 3.3. Endemic Equilibrium

#### 3.3.1. Existence of the Endemic Equilibrium

IfR0>1, system 3 has a unique endemic equilibrium E*(S*,H*,P*,A*,T*,R*), where
(31)S*=ΛR0(μ1+d),H*=(bm-gk4)P*gk3+mk1,P*=(R0-1)(μ1+d)(gk3+mk1)β1(bm-gk4)+β2(gk3+mk1),A*=k2P*+α5T*n,T*=(bk3+k1k4)P*gk3+mk1,R*=μ1ΛdR0(μ1+d).


#### 3.3.2. Global Stability of the Endemic Equilibrium

Theorem 4.
IfR0>1, k3=0, the endemic equilibrium E* is globally asymptotically stable.Proof.
To study the globally asymptotic stability of the endemic equilibrium, motivated by [22, 23], we use a Lyapunov function V as follows:
(32)V=(S-S*ln⁡S)+B(H-H*ln⁡H)+C(P-P*ln⁡P)+DT-T*ln⁡T.
Applying the replacements
(33)SS*=x,HH*=y,PP*=z,TT*=μ,
the derivative of V is given by
(34)V˙=dSdt-S*dSSdt+BdHdt-H*dHHdt+CdPdt-P*dPPdt+DdTdt-T*dTTdt=dSdt1-S*S+B1-H*HdHdt+C1-P*PdPdt+D1-T*TdTdt=1-1xΛ-β1SH-β2SP-dS-μ1S+B1-1yβ1SH+β2SP-aH+hT+C1-1zk1H-bP+gT+D1-1μk4P-mT.
Since E(t) satisfies the equation of the model with system 3,
(35)Λ-β1S*H*-β2S*P*-dS*-μ1S*=0,β1S*H*+β2S*P*-aH*+hT*=0,k1H*-bP*+gT*=0,k2P*+α5T*-nA*=0,k3H*+k4P*-mT*=0,μ1S*-dR*=0,V˙=1-1xβ1S*H*+β2S*P*+dS*+μ1S*hhhhhhhhhh-β1S*H*xy-β2S*P*xz-dS*x-μ1S*x+B1-1yβ1S*H*xy+β2S*P*xzhhhhhhhhhhhhh-aH*y+hT*μ+C1-1zk1H*y-bP*z+gT*μ+D1-1μk4P*z-mT*μ=-1-x2x(d+μ1)S*+β1S*H*1-xy+β2S*P*1-xz·1-1x+B1-1yβ1S*H*xy+β2S*P*xzβ1S*H*+β2S*P*+hT*H*hhhhhhhhhhhhh-β1S*H*+β2S*P*+hT*H*H*yhhhhhhhhhhhhhβ1S*H*+β2S*P*+hT*H*+hT*μ+C1-1zk1H*y-k1H*+gT*P*P*z+gT*μ+D1-1μk4P*z-k4P*T*T*μ=-1-x2xd+μ1S*+β1S*H*1-xy-1x+y+β2S*P*1-xz-1x+z+B1-1y·β1S*H*xy-y+β2S*P*xz-yhhhh+hT*μ-y+C1-1zk1H*y-z+gT*μ-z+Dk4P*1-1μz-μ=-1-x2xd+μ1S*+β1S*H*1-xy-1x+y+β2S*P*1-xz-1x+z+Bβ1S*H*xy-x-y+1+Bβ2S*P*xz-xzy-y+1+BhT*μ-y-μy+1+Ck1H*y-z-yz+1+CgT*μ-z-μz+1+Dk4P*z-μ-zμ+1=-1-x2xd+μ1S*+B-1β1S*H*xy+B-1β2S*P*xz+β1S*H*-Bβ1S*H*-Bβ2S*P*hhhh-BhT*+Ck1H*y+β2S*P*-Ck1H*z+(-CgT*+Dk4P*)z+BhT*+CgT*-Dk4P*μ+(1+B)β1S*H*+(1+B)β2S*P*+BhT*+Ck1H*+CgT*+Dk4P*-β1S*H*x-Bβ1S*H*x-β2S*P*x-xzyBβ2S*P*-BhT*μy-yzCk1H*-μzCgT*-zμDk4P*.
We choose the coefficients of xy, xz, y, z, and μ which are equal to zero; that is,
(36)B-1=0,BhT*+CgT*-Dk4P*=0,β2S*P*-Ck1H*-CgT*+Dk4P*=0,β1S*H*-Bβ1S*H*-Bβ2S*P*-BhT*+Ck1H*=0.
From 36, we have
(37)B=1,C=β2S*P*+hT*k1H*,D=k1H*hT*+(β2S*P*+hT*)gT*k1k4H*P*.
From 37, the derivative of V is given by
(38)V˙=-1-x2xd+μ1S*+2β1S*H*+2β2S*P*+hT*+Ck1H*+CgT*+Dk4P*-β1S*H*x-β1S*H*x-β2S*P*x-xzyβ2S*P*-hT*μy-yzCk1H*-μzCgT*-zμDk4P*=2-x-1xβ1S*H*+3-1x-xzy-yzβ2S*P*+3-yz-μy-zμhT*+2-μz-zμCgT*≤0.
According to [24], since the arithmetical mean is greater than or equal to the geometrical mean, then 2-x-1/x≤0 for x>0 and 2-x-1/x=0 if and only if x=1; 3-1/x-xz/y-y/z=0 for x>0, y>0, and z>0 if and only if x=y=z=1, and 3-y/z-μ/y-z/μ=0 if and only if μ=y=z. Futhermore, V˙=0 if and only if x=y=z=μ=1. The maximum invariant set of system 3 on the set {(x,y,z,μ):V˙=0} is the singleton (1,1,1,1). Thus, for system 3, the endemic equilibrium is globally asymptotically stable if R0>1 by LaSalle invariance principle.Remark 5.
Biologically speaking,k3=0 implies that the HIV-positive individuals in the asymptomatic stage of HIV infection do not receive treatment. Theorem 4 tells us that when HIV-positive individuals in the asymptomatic stage of HIV infection do not receive treatment, the endemic equilibrium E* is globally asymptotically stable. This means that AIDS is persistent. When k3≠0, the same result is validated by numerical simulations (see Figure 4).


## 3.1. Disease-Free Equilibrium and the Reproduction Number

The model has a disease-free equilibrium given by(15)E0=Λμ1+d,0,0,0,0,Λμ1dμ1+d.
In the following, the basic reproduction number of system 3 will be obtained by the next generation matrix method formulated in [17, 18].Lettingx=(H,P,A,T,S,R)T, then system 3 can be written as
(16)dxdt=F(x)-V(x),
where
(17)F(x)=β1SH+β2SP00000,V(x)=aH-hTbP-gT-k1HnA-k2P-α5TmT-k3H-k4Pβ1SH+β2SP+dS+μ1S-ΛdR-μ1S.
The Jacobian matrices of F(x) and V(x) at the free equilibrium E0 are, respectively,
(18)DF(E0)=F4×400000000,DV(E0)=V4×404×104×1β1Λd+μ1β2Λd+μ100d+μ1-μ10000d,
where
(19)F=β1Λd+μ1β2Λd+μ100000000000000,V=a00-h-k1b0-g0-k2n-α5-k3-k40m.
The model reproduction number, denoted by R0, is thus given by
(20)R0=ρ(FV-1)=β1Λ(bm-k4g)+β2Λ(k1m+k3g)(d+μ1)(abm-k1k4h-k3hb-k4ga).


## 3.2. Global Stability ofE0

Theorem 3.
The disease-free equilibriumE0 is globally asymptotically stable for R0<1 and unstable otherwise.Proof.
According to [19, 20], we introduce the Lyapunov candidate function
(21)V(t)=H+ΘP+cT,Θ>0,c>0.
Differentiating V(t) with respect to time yields
(22)V˙=H˙+ΘP˙+cT˙.
Substituting the model equations into 22, we get
(23)V˙=β1SH+β2SP-aH+hT+Θ(k1H-bP+gT)+ck3H+k4P-mT.
With S0=Λ/(d+μ1) and R¯0=Λμ1/d(d+μ1), we have
(24)V˙≤β1ΛHd+μ1+β2ΛPd+μ1-aH+hT+Θk1H-bP+gT+ck3H+k4P-mT=β1Λd+μ1-a+k1Θ+ck3H+β2Λd+μ1-bΘ+ck4P+h+Θg-cmT.
In order to get the value of c and Θ, we can use
(25)β2Λd+μ1-bΘ+ck4=0,h+Θg-cm=0.
From 25, we can calculate
(26)c=β2Λg+hb(d+μ1)(bm-k4g)(d+μ1),Θ=β2Λm+k4h(d+μ1)(bm-k4g)(d+μ1).
Taking c and Θ into 24, we can get
(27)V˙≤β1Λd+μ1-a+k1β2Λm+k1k4h(d+μ1)(bm-k4g)(d+μ1)hhh+β2Λgk3+hbk3d+μ1bm-k4gd+μ1H=β1Λd+μ1(bm-k4g)+β2Λd+μ1(mk1+gk3)hhhhh-abm+ak4g+k1k4h+k3hbβ2Λd+μ1hhh·bm-k4g-1β1Λd+μ1H.
Since
(28)R0=ρ(FV-1)=β1Λ(bm-k4g)+β2Λk1m+k3gd+μ1abm-k1k4h-k3hb-k4ga,
we can obtain
(29)R0-1=β1Λd+μ1(bm-k4g)+β2Λd+μ1(mk1+gk3)hhhh-abm+ak4g+k1k4h+k3hbβ1Λd+μ1·abm-k1k4h-k3hb-k4ga-1,
when R0≤1,
(30)V˙≤abm-k1k4h-k3hb-k4gaR0-1bm-k4gH≤0.
WhileV˙=0, only H=T=P=0. Substituting H=T=P=0 into system 3, we can have S→Λ/(d+μ1), A→0, and R→Λμ1/d(d+μ1) as t→∞. According to LaSalle’s invariance principle [21], the disease-free equilibrium is globally asymptotically stable.


## 3.3. Endemic Equilibrium

### 3.3.1. Existence of the Endemic Equilibrium

IfR0>1, system 3 has a unique endemic equilibrium E*(S*,H*,P*,A*,T*,R*), where
(31)S*=ΛR0(μ1+d),H*=(bm-gk4)P*gk3+mk1,P*=(R0-1)(μ1+d)(gk3+mk1)β1(bm-gk4)+β2(gk3+mk1),A*=k2P*+α5T*n,T*=(bk3+k1k4)P*gk3+mk1,R*=μ1ΛdR0(μ1+d).


### 3.3.2. Global Stability of the Endemic Equilibrium

Theorem 4.
IfR0>1, k3=0, the endemic equilibrium E* is globally asymptotically stable.Proof.
To study the globally asymptotic stability of the endemic equilibrium, motivated by [22, 23], we use a Lyapunov function V as follows:
(32)V=(S-S*ln⁡S)+B(H-H*ln⁡H)+C(P-P*ln⁡P)+DT-T*ln⁡T.
Applying the replacements
(33)SS*=x,HH*=y,PP*=z,TT*=μ,
the derivative of V is given by
(34)V˙=dSdt-S*dSSdt+BdHdt-H*dHHdt+CdPdt-P*dPPdt+DdTdt-T*dTTdt=dSdt1-S*S+B1-H*HdHdt+C1-P*PdPdt+D1-T*TdTdt=1-1xΛ-β1SH-β2SP-dS-μ1S+B1-1yβ1SH+β2SP-aH+hT+C1-1zk1H-bP+gT+D1-1μk4P-mT.
Since E(t) satisfies the equation of the model with system 3,
(35)Λ-β1S*H*-β2S*P*-dS*-μ1S*=0,β1S*H*+β2S*P*-aH*+hT*=0,k1H*-bP*+gT*=0,k2P*+α5T*-nA*=0,k3H*+k4P*-mT*=0,μ1S*-dR*=0,V˙=1-1xβ1S*H*+β2S*P*+dS*+μ1S*hhhhhhhhhh-β1S*H*xy-β2S*P*xz-dS*x-μ1S*x+B1-1yβ1S*H*xy+β2S*P*xzhhhhhhhhhhhhh-aH*y+hT*μ+C1-1zk1H*y-bP*z+gT*μ+D1-1μk4P*z-mT*μ=-1-x2x(d+μ1)S*+β1S*H*1-xy+β2S*P*1-xz·1-1x+B1-1yβ1S*H*xy+β2S*P*xzβ1S*H*+β2S*P*+hT*H*hhhhhhhhhhhhh-β1S*H*+β2S*P*+hT*H*H*yhhhhhhhhhhhhhβ1S*H*+β2S*P*+hT*H*+hT*μ+C1-1zk1H*y-k1H*+gT*P*P*z+gT*μ+D1-1μk4P*z-k4P*T*T*μ=-1-x2xd+μ1S*+β1S*H*1-xy-1x+y+β2S*P*1-xz-1x+z+B1-1y·β1S*H*xy-y+β2S*P*xz-yhhhh+hT*μ-y+C1-1zk1H*y-z+gT*μ-z+Dk4P*1-1μz-μ=-1-x2xd+μ1S*+β1S*H*1-xy-1x+y+β2S*P*1-xz-1x+z+Bβ1S*H*xy-x-y+1+Bβ2S*P*xz-xzy-y+1+BhT*μ-y-μy+1+Ck1H*y-z-yz+1+CgT*μ-z-μz+1+Dk4P*z-μ-zμ+1=-1-x2xd+μ1S*+B-1β1S*H*xy+B-1β2S*P*xz+β1S*H*-Bβ1S*H*-Bβ2S*P*hhhh-BhT*+Ck1H*y+β2S*P*-Ck1H*z+(-CgT*+Dk4P*)z+BhT*+CgT*-Dk4P*μ+(1+B)β1S*H*+(1+B)β2S*P*+BhT*+Ck1H*+CgT*+Dk4P*-β1S*H*x-Bβ1S*H*x-β2S*P*x-xzyBβ2S*P*-BhT*μy-yzCk1H*-μzCgT*-zμDk4P*.
We choose the coefficients of xy, xz, y, z, and μ which are equal to zero; that is,
(36)B-1=0,BhT*+CgT*-Dk4P*=0,β2S*P*-Ck1H*-CgT*+Dk4P*=0,β1S*H*-Bβ1S*H*-Bβ2S*P*-BhT*+Ck1H*=0.
From 36, we have
(37)B=1,C=β2S*P*+hT*k1H*,D=k1H*hT*+(β2S*P*+hT*)gT*k1k4H*P*.
From 37, the derivative of V is given by
(38)V˙=-1-x2xd+μ1S*+2β1S*H*+2β2S*P*+hT*+Ck1H*+CgT*+Dk4P*-β1S*H*x-β1S*H*x-β2S*P*x-xzyβ2S*P*-hT*μy-yzCk1H*-μzCgT*-zμDk4P*=2-x-1xβ1S*H*+3-1x-xzy-yzβ2S*P*+3-yz-μy-zμhT*+2-μz-zμCgT*≤0.
According to [24], since the arithmetical mean is greater than or equal to the geometrical mean, then 2-x-1/x≤0 for x>0 and 2-x-1/x=0 if and only if x=1; 3-1/x-xz/y-y/z=0 for x>0, y>0, and z>0 if and only if x=y=z=1, and 3-y/z-μ/y-z/μ=0 if and only if μ=y=z. Futhermore, V˙=0 if and only if x=y=z=μ=1. The maximum invariant set of system 3 on the set {(x,y,z,μ):V˙=0} is the singleton (1,1,1,1). Thus, for system 3, the endemic equilibrium is globally asymptotically stable if R0>1 by LaSalle invariance principle.Remark 5.
Biologically speaking,k3=0 implies that the HIV-positive individuals in the asymptomatic stage of HIV infection do not receive treatment. Theorem 4 tells us that when HIV-positive individuals in the asymptomatic stage of HIV infection do not receive treatment, the endemic equilibrium E* is globally asymptotically stable. This means that AIDS is persistent. When k3≠0, the same result is validated by numerical simulations (see Figure 4).


## 3.3.1. Existence of the Endemic Equilibrium

IfR0>1, system 3 has a unique endemic equilibrium E*(S*,H*,P*,A*,T*,R*), where
(31)S*=ΛR0(μ1+d),H*=(bm-gk4)P*gk3+mk1,P*=(R0-1)(μ1+d)(gk3+mk1)β1(bm-gk4)+β2(gk3+mk1),A*=k2P*+α5T*n,T*=(bk3+k1k4)P*gk3+mk1,R*=μ1ΛdR0(μ1+d).


## 3.3.2. Global Stability of the Endemic Equilibrium

Theorem 4.
IfR0>1, k3=0, the endemic equilibrium E* is globally asymptotically stable.Proof.
To study the globally asymptotic stability of the endemic equilibrium, motivated by [22, 23], we use a Lyapunov function V as follows:
(32)V=(S-S*ln⁡S)+B(H-H*ln⁡H)+C(P-P*ln⁡P)+DT-T*ln⁡T.
Applying the replacements
(33)SS*=x,HH*=y,PP*=z,TT*=μ,
the derivative of V is given by
(34)V˙=dSdt-S*dSSdt+BdHdt-H*dHHdt+CdPdt-P*dPPdt+DdTdt-T*dTTdt=dSdt1-S*S+B1-H*HdHdt+C1-P*PdPdt+D1-T*TdTdt=1-1xΛ-β1SH-β2SP-dS-μ1S+B1-1yβ1SH+β2SP-aH+hT+C1-1zk1H-bP+gT+D1-1μk4P-mT.
Since E(t) satisfies the equation of the model with system 3,
(35)Λ-β1S*H*-β2S*P*-dS*-μ1S*=0,β1S*H*+β2S*P*-aH*+hT*=0,k1H*-bP*+gT*=0,k2P*+α5T*-nA*=0,k3H*+k4P*-mT*=0,μ1S*-dR*=0,V˙=1-1xβ1S*H*+β2S*P*+dS*+μ1S*hhhhhhhhhh-β1S*H*xy-β2S*P*xz-dS*x-μ1S*x+B1-1yβ1S*H*xy+β2S*P*xzhhhhhhhhhhhhh-aH*y+hT*μ+C1-1zk1H*y-bP*z+gT*μ+D1-1μk4P*z-mT*μ=-1-x2x(d+μ1)S*+β1S*H*1-xy+β2S*P*1-xz·1-1x+B1-1yβ1S*H*xy+β2S*P*xzβ1S*H*+β2S*P*+hT*H*hhhhhhhhhhhhh-β1S*H*+β2S*P*+hT*H*H*yhhhhhhhhhhhhhβ1S*H*+β2S*P*+hT*H*+hT*μ+C1-1zk1H*y-k1H*+gT*P*P*z+gT*μ+D1-1μk4P*z-k4P*T*T*μ=-1-x2xd+μ1S*+β1S*H*1-xy-1x+y+β2S*P*1-xz-1x+z+B1-1y·β1S*H*xy-y+β2S*P*xz-yhhhh+hT*μ-y+C1-1zk1H*y-z+gT*μ-z+Dk4P*1-1μz-μ=-1-x2xd+μ1S*+β1S*H*1-xy-1x+y+β2S*P*1-xz-1x+z+Bβ1S*H*xy-x-y+1+Bβ2S*P*xz-xzy-y+1+BhT*μ-y-μy+1+Ck1H*y-z-yz+1+CgT*μ-z-μz+1+Dk4P*z-μ-zμ+1=-1-x2xd+μ1S*+B-1β1S*H*xy+B-1β2S*P*xz+β1S*H*-Bβ1S*H*-Bβ2S*P*hhhh-BhT*+Ck1H*y+β2S*P*-Ck1H*z+(-CgT*+Dk4P*)z+BhT*+CgT*-Dk4P*μ+(1+B)β1S*H*+(1+B)β2S*P*+BhT*+Ck1H*+CgT*+Dk4P*-β1S*H*x-Bβ1S*H*x-β2S*P*x-xzyBβ2S*P*-BhT*μy-yzCk1H*-μzCgT*-zμDk4P*.
We choose the coefficients of xy, xz, y, z, and μ which are equal to zero; that is,
(36)B-1=0,BhT*+CgT*-Dk4P*=0,β2S*P*-Ck1H*-CgT*+Dk4P*=0,β1S*H*-Bβ1S*H*-Bβ2S*P*-BhT*+Ck1H*=0.
From 36, we have
(37)B=1,C=β2S*P*+hT*k1H*,D=k1H*hT*+(β2S*P*+hT*)gT*k1k4H*P*.
From 37, the derivative of V is given by
(38)V˙=-1-x2xd+μ1S*+2β1S*H*+2β2S*P*+hT*+Ck1H*+CgT*+Dk4P*-β1S*H*x-β1S*H*x-β2S*P*x-xzyβ2S*P*-hT*μy-yzCk1H*-μzCgT*-zμDk4P*=2-x-1xβ1S*H*+3-1x-xzy-yzβ2S*P*+3-yz-μy-zμhT*+2-μz-zμCgT*≤0.
According to [24], since the arithmetical mean is greater than or equal to the geometrical mean, then 2-x-1/x≤0 for x>0 and 2-x-1/x=0 if and only if x=1; 3-1/x-xz/y-y/z=0 for x>0, y>0, and z>0 if and only if x=y=z=1, and 3-y/z-μ/y-z/μ=0 if and only if μ=y=z. Futhermore, V˙=0 if and only if x=y=z=μ=1. The maximum invariant set of system 3 on the set {(x,y,z,μ):V˙=0} is the singleton (1,1,1,1). Thus, for system 3, the endemic equilibrium is globally asymptotically stable if R0>1 by LaSalle invariance principle.Remark 5.
Biologically speaking,k3=0 implies that the HIV-positive individuals in the asymptomatic stage of HIV infection do not receive treatment. Theorem 4 tells us that when HIV-positive individuals in the asymptomatic stage of HIV infection do not receive treatment, the endemic equilibrium E* is globally asymptotically stable. This means that AIDS is persistent. When k3≠0, the same result is validated by numerical simulations (see Figure 4).


## 4. Numerical Simulation

In this section, we present some numerical simulations of system3 to illustrate our results. To illustrate the analytic results obtained above, we give some simulations using the parameter values in Table 1.First, whenk3≠0, we choose k3=0.6, k4=0.4, numerical simulation gives R0=0.95047<1, and the disease-free equilibrium E0 is globally asymptotically stable (see Figure 2). Figure 2 indicates the effect of treatment on stability of the equilibria of the HIV/AIDS infection model. At first, The number of patients who do not receive treatment is increasing constantly. After treatment, The number of patients is decreasing and tends to be zero at the end, which indicates that the spread of infection can be slowed down.Figure 2
Whenk3=0.6, k4=0.4, and R0=0.95047<1, the disease-free equilibrium is globally asymptotically stable.Second, we choosek3=0.51086, k4=0.4, numerical simulation gives R0=1, and the disease-free equilibrium E0 is globally asymptotically stable (Figure 3). If we give a proper treatment, the number of individuals with AIDS/HIV will be controlled to some extent. Hence, the disease will be eradicated in the future.Figure 3
Whenk3=0.51086, k4=0.4, and R0=1, the disease-free equilibrium is globally asymptotically stable.Figure 4
Whenk3=0.05, k4=0.4, and R0=1.412>1, the endemic equilibrium is globally asymptotically stable.Third, we choosek3=0.05, k4=0.4, numerical simulation gives R0=1.412, and the endemic equilibrium is globally asymptotically stable (Figure 4). Figure 4 indicates that when the basic reproduction is greater than one, as time goes on, the number of people finally tends to be constant, which means that the endemic equilibrium is globally asymptotically stable.Finally, we show the relation betweenk3, k4, and R0. When the value of k3 and k4 is increasing at the time, the basic reproduction is reducing sharply. What is more, Figure 5 confirms that, with the increase in the rate of treatment, the number of individuals with HIV-positive and full-blown AIDS but not receiving treatment and who are receiving ARV treatment tends to zero, which indicates that treatment is important and can prolong their lives.Figure 5
The relationship amongR0, k3, and k4.


## 5. Discussion

A mathematical model is proposed and analyzed to study the spread of HIV/AIDS with treatment. In this paper we have analyzed a stage structured model for HIV and the effect of treatment has also been studied. In the presented model, we get two equilibria: the disease-free equilibrium and the endemic equilibrium. We further consider global asymptotic stability of the disease-free equilibrium by using the well-known Lyapunov-LaSalle invariance principal. It is found that the disease-free equilibrium is globally asymptotically stable when the basic reproduction number is less than one. When the basic reproduction number is greater than one andk3=0, the endemic equilibrium is globally asymptotically stable. Our results show that early treatment for individuals in asymptomatic stage of HIV infection or the pre-AIDS stage is very important.



---

*Source: 630503-2015-05-14.xml*