Pure and Applied Mathematics Journal
Volume 5, Issue 4, August 2016, Pages: 97-102

Reduction of Mortality Rate Due to AIDS When Treatment Is Considered

Udoy S. Basak

Department of Mathematics, Pabna University of Science & Technology, Pabna, Bangladesh

Email address:

To cite this article:

Udoy S. Basak. Reduction of Mortality Rate Due to AIDS When Treatment Is Considered. Pure and Applied Mathematics Journal. Vol. 5, No. 4, 2016, pp. 97-102. doi: 10.11648/j.pamj.20160504.12

Received: May 16, 2016; Accepted: May 28, 2016; Published: June 18, 2016

Abstract: AIDS is one of the most threatening diseases for human being that is caused by a virus named HIV. Here the reduction of the death rate after infected by AIDS has been discussed. A mathematical model of HIV has been formulated. Then its positivity and boundedness has been investigated. It has been shown that it is possible to minimize the mortality rate by providing the treatment to the HIV infected people. Moreover, the control of the transfer rate from the infected class to the AIDS class reduces the disease rate. The increasing of the transfer rate from the infected class to the treated class also reduces the mortality rate.

Keywords: HIV, AIDS, Reproduction Number, Endemic Equilibrium Point

1. Introduction

In the last decades or so mathematical models in epidemiology are playing a very important role in analyzing the spread and control of the infectious disease that are causing a lot of human death every year. In the past, human being have been threatened by many infectious disease, most of them are caused by viruses [1]. Among them HIV is responsible for a lot of human death. Till now people are dying by this virus.

Human Immunodeficiency Virus (HIV) is a member of Lent Virus sub member family of retrovirus that produces chronic infections in the host and gradually damages the host`s immune system [2]. There are two types of HIV viruses: HIV-1 and HIV-2. Among them, HIV-1 is more dangerous. It easily transmits among the people and causes a lot of HIV infections globally every year [3]. Immediately after HIV infection, people cannot determine his/her infection. In this time, HIV infection may not show up on an HIV test, but people who have it are highly infectious and can spread the infection to others. People with chronic HIV infection may not have any HIV-related symptoms or only mild ones [4].

After being infected, people who aren’t taking medicine to treat HIV, this period can last a decade or longer, but some may progress through this phase faster. People who are taking medicine to treat HIV the right way, every day may be in this stage for several decades because treatment helps keep the virus in check.

HIV is treated using a combination of medicines to fight HIV infection. This is called antiretroviral therapy (ART). ART isn’t a cure, but it can control the virus so that one can live a longer, healthier life and reduce the risk of transmitting HIV to others. ART is recommended for all people with HIV, regardless of how long they’ve had the virus or how healthy they are. If left untreated, HIV will attack the immune system and eventually progress to AIDS.

2. Formulation of the Model

The total sexually-active population at time t, denoted by  is subdivided into mutually-exclusive compartments, namely susceptible , exposed to HIV but show no clinical symptoms of the disease  HIV infected individuals with clinical symptoms of HIV , infected individuals who are taking treatment , individuals having AIDS , so that

Let us suppose that the susceptible population is increasing at a constant rate. Let the rate be

Susceptible individuals acquire HIV infection, following effective contact with people infected with HIV only (i.e. those in the E, I and T classes) at a rate λ, where

Here β is the transmission rate for HIV and  are known as the modification parameter. These modification parameters indicate that the infected individuals are more infectious than the treated class. Among these class, AIDS class is assumed to spread more infection than other two classes. In the HIV exposed class , individuals entering into this class from susceptible class  at a rate λ. In the HIV exposed class  individuals progress to infected class  at a rate . Individual also progresses to the treatment class ) at a rate And Individual also progresses to the AIDS class) at a rate

In the treatment class, individuals progress to AIDS class at a rate. Also individuals progress to the AIDS and HIV treatment class  at a rate. In this class, let  denotes the death rate due to the disease and  be the natural mortality rate.

Combining all the aforementioned assumption and definitions, the model becomes:


Where, .

Schematically this can be shown as follows:

Figure 1. Diagram of the model.

3. Analysis of the Model

3.1. Boundedness Region

Here we have used some parameters (explained in the table 2) to described the model (1). These variables are non-negative. Also the state variables (described in table 1) are non-negative for all time

The total population is denoted by  Total population can be written as:

The equation written above, is changing at a rate


In absence of disease, i.e. for  we have


By the separation of variables we can write

After integration we have


Here  is the integrating constant which is to be determined. Let at  From the Eq. (4) we have

Putting the value of  in Eq. (4) we have



Therefore, the region  contains the feasible solution of the model (1).

3.2. Positivity of Solution

Now we need to check the positivity of the solution of the model (1). Let the initial data be Now we have to show that the solution set  of the system (1) is positive.

Consider the 1st equation of the model (1)

After integration we can write


Similarly from the 2nd, 3rd, 4th and 5th equations respectively we have





Thus the solution set  is positive for

4. Disease-Free Equilibrium Point

Disease free equilibrium point is denoted by. To obtain the disease-free equilibrium points on , we need to put  (since in the disease-free equilibrium there is no disease). The positive DFE for this model is . So we can conclude that the model has a disease-free equilibrium point and the disease-free equilibrium point of this model is .

5. Local Stability of the Disease-Free Equilibrium

The basic reproduction number which is denoted by  plays a vital role in the stability of the model. The number of new infections generated by an infected individual is known as the basic reproduction number. The associated nonnegative matrix for the new infection terms and the non-singular-matrix,  for the remaining transfer terms are given respectively by

And  has the following expression

The basic reproduction number  is the spectral radius of the matrix . The eigenvalues of the matrix  are

Lemma [5]. Local stability of the disease-free equilibrium depends on the basic reproduction number. If the basic reproduction number is bigger than one, then disease-free equilibrium point locally asymptotically unstable and if the basic reproduction number is less than one, then it is locally asymptotically unstable.

6. Global Stability of the Disease-Free Equilibrium

Here we have used a technique established by Castillo-Chavez [6] to investigate the global stability of the disease-free equilibrium point. We rewrite the model as:


Where  and  Here the components of  denotes the infected population. The disease free equilibrium is now denoted as:

The condition must be met to guarantee a local asymptotic stability:

Here,  is globally asymptotically stable (GAS).


Where is an -matrix (the off-diagonal elements of  are non-negative) and  is the region where the model makes biological sense. If the system satisfies the conditions of then the theorem below holds:

Theorem [5]. The fixed point  is a globally asymptotically stable equilibrium of the system (6) provided that and the assumptions in (4) are satisfied.

 Form the model system (1) and (6), we have


Here  &

So the conditions are fulfilled. So is globally asymptotically stable when .

7. Endemic Equilibrium of the Model

We have used a well-known technique known as the Central Manifold Theorem [6] to check the stability of the endemic equilibrium point of the model. To apply this theorem we make the following change of variables.


The model (1) can be rewritten in the form:

Where  and  as



The Jacobean of the system (8) is

To analyze the dynamics of (8), we compute the eigenvalues of the Jacobian of (8) at the disease free equilibrium (DFE). It can be shown that this Jacobian has a right eigenvector given by:


And the left eigenvectors are given by

It can be shown after some manipulations involving the evaluation of the associated non-vanishing partial derivatives of  that


Since  which is less than zero and  which is positive.


Thus we have

Again  are all positive and all the parameters are positive. So

Thus we have established the following theorem:

Theorem. Thus we have the following results: the endemic equilibrium point is locally stable when the basic reproduction number is bigger than one and unstable when basic reproduction number is less than one.

8. Numerical Simulations and Discussions

The description of the variable is shown in the Table 1. In the model we have used some parameters. The values of the parameters that are used in the model are shown in the table 2.

Table 1. Description of variables of the model.

Table 2. The values of the parameters of the model.

Figure 2 shows that the total infected population is increasing when the basic reproduction number  and figure 3 shows that the infected population is decreasing when  Here the value of the basic reproduction number  depends on the transfer rate  from the infected class  to treated class . The values of the reproduction number  decreases with the transfer rate . Thus we can minimize the mortality rate by increasing the number of the treated population. Figure 4 shows that, the prevalence is decreasing when the basic reproduction number  and figure 5 shows that it is increasing when

Figure 2. Total infection when

Figure 3. Total Infection when

Figure 4. Prevalence when

Figure 5. Prevalence when

9. Conclusions

AIDS is one of the devastating diseases for human being. Mathematical models are used to estimate the dynamics of the diseases. Here we have formulated a model describing the HIV with treatment. The findings of the study are listed below:

We can minimize the mortality rate by providing the treatment to the HIV infected people.

The basic reproduction number highly depends on , which is the transfer rate from  to  class. Thus the reduction in the transfer rate from the infected rate  to the AIDS class  will reduce the mortality rate.

The increasing of the transfer rate  from the infected class  to the treated class  also reduces the disease burden.


  1. F. Brauer and C. Castillo-Chavez, Mathematical Models in Population Biology and Epidemiology, Springer-Verlag. New York, 2001.
  2. Hu et al (1996), "The emerging genetic diversity of HIV, the importance of global surveillance for diagnostic research and prevention". JAMA, 275 (3), 210-216.
  3. Reeves JD, Doms RW (2002). "Human immunodeficiency virus type 2". The Journal of general virology 83 (Pt 6): 1253–65. doi:10.1099/vir.0.18253-0 (inactive 2015-01-13). PMID 12029140.
  4. Santiago et al, "Simian Immunodeficiency Virus Infection in Free-Ranging Sooty Mangabeys (Cercocebus atys atys) from the Tai Forest, Cote d'Ivoire: Implications for the Origin of Epidemic Human Immunodeficiency Virus Type 2", Journal of Virology 79 (19): 12515–27.
  5. Chavez et al (2004), "Dynamical model of tuberculosis and their applications". Math. Bioscience., 1: 361-404.
  6. Lawi et al (2011), "Mathematical model for malaria and meningitis co-infection among children". Applied Mathematics Sciences, Vol. 5, 2011, no. 47, 2337-2359.
  7. Weiss RA (May 1993). "How does HIV cause AIDS?". Science 260(5112): 1273-1279.
  8. Douek DC, Roederer M, Koup RA (2009). "Emerging Concepts in the Immunopathogenesis of AIDS". Annu. Rev. Med. 60: 471–84.
  9. Reeves JD, Doms RW (2002). "Human immunodeficiency virus type 2". The Journal of general virology 83 (Pt 6): 1253–65. doi:10.1099/vir.0.18253-0 (inactive 2015-01-13). PMID 12029140.
  10. Garg H, Mohl J, Joshi A (Nov 9, 2012). "HIV-1 induced bystander apoptosis". Viruses 4 (11): 3020–43. doi: 10.3390/v4113020. PMC 3509682. PMID 23202514.
  11. Udoy S. Basak et al. "Mathematical Analysis of an HIV/AIDS Epidemic Model". American Journal of Mathematics and Statistics 2015, 5(5): 253-258 DOI: 10.5923/j.ajms.20150505.05.

Article Tools
Follow on us
Science Publishing Group
NEW YORK, NY 10018
Tel: (001)347-688-8931