1. Introduction
AIDS is the most deadly disease caused by a human immunodeficiency virus (HIV). The virus destroys all the immune system and leaves individuals susceptible to any other infections. It multiplies inside lymphocytes and finally destroys them. When the lymphocytes are reduced to a certain numbers, the immune system stops functioning correctly. Therefore, the individual can catch any kind of disease that might kill him easily because of the failure of the immune system. However, there exist drugs that can slow down the evolution of the virus. HIV is usually transmitted in three different ways: sexual contacts, blood transfusion, and exchange between mother and child during pregnancy, childbirth and breastfeeding.
Many mathematical models are used to study the impact of preventive control strategies on the spread of HIV-AIDS in given populations (cf. [1] -[11] , etc.). Some of these models showed that a change in risky behaviour was necessary to prevent the spread of HIV even in the presence of a treatment (see for example [12] -[16] ). Thus, it is instructive to study models that focus on non-pharmaceutical interventions, such as the use of public health sensitization campaign.
The models developed in [15] [17] -[19] study the impact of public health sensitization campaign on the spread of HIV-AIDS. In this paper we propose and study a mathematical model to estimate the impact of public health sensitization campaign on the spread of HIV-AIDS in Mali. We divide for it the population into two classes: “class with high-risk behavioral or class without public health sensitization” and “class with low risk behavioral or class with public health sensitization”. Every class consists of susceptible individuals and infected individuals. The class of the individuals at high-risk behavioral is split into susceptibles individuals, indi-
viduals who are in stage 1 of the infection, individuals who are in stage 2 of the infection and
individuals who are in stage AIDS while the class of the individuals at low-risk behavioral is split into susceptibles individuals, individuals who are in stage 1 of the infection, individuals who are in
stage 2 of the infection and individuals who are in stage AIDS. The total population (Figure 1) at
time t is denoted by and can be expressed as the following sum:
(1)
Our model is given by the following system of ODEs with constant coefficients:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
Figure 1. Behavioral representation of the HIV-AIDS model.
(9)
with
(10)
By adding of (2) to (9), we obtain: where the parameters of the model are defined in Table 1.
Our mathematical model is an extension of the models developed in [15] [17] -[19] . In our model, we suppose:
H1: that he mode of transmission of the virus is the horizontal transmission;
H2: that every individual is susceptible at high risk before his recruitment in the compartment and that the rate of mortality induced by the HIV is neglected;
H3: that, and are in, the parameters, , , , , , , , , , , are in and that.
2. Analysis of the Complete Model
2.1. Existence of Solutions
To show that the model is mathematically and biologically possible, we begin by rewriting it in terms of proportions. So, we introduce the following scalings:
(11)
Consequently:
(12)
By using what precedes, the rates of infection (10) become:
(13)
If we introduce the following parameters:
(14)
In the new variables, (2)-(9) reduces to:
(15)
(16)
(17)
(18)
(19)
(20)
(21)
(22)
We suppose that the initial conditions belong in where
(23)
Now we can enounce the following result:
Theorem 1. For any initial condition in, the system has a unique solution globally defined and which stays in for any time.
Before giving the proof of this theorem, we give at first a technical result which we shall use after.
Lemma 1. Let a(t) and y(t) be n X n matrices of bounded measurable functions on, if
(24)
with for and then for all.
Proof. Indeed, this follows from the integrated form of the differential Equation (24),
□
We rewrite the system (15)-(22) in the form
(25)
Now we can give the proof of the theorem 1.
Proof. Step 1: Local existence of the solutions.
The local existence of the solutions ensues directly from the regularity of the function which is of class in.
Step 2: We show that is positively invariant.
A., , , , , , , ,. Let’s suppose that it exists such, and and let’s rewrite the Equations (15), (17)-(22) in the form:
(26)
(27)
(28)
(29)
(30)
(31)
By Lemma 1, , , , , et for.
We next show that and remains positive for all.
Proceding by contradiction:
We suppose that and for and.
Then, and. By considering the Equation (26) in time, we have:
(32)
(33)
(34)
(35)
(36)
(37)
By Lemma 1, , , , , and for.
Now we consider the Equations (16) and (20) in time
(38)
(39)
which is a contradiction, Consequently and remains positive for all.
B. The following inequalities hold:
Adding all the equations of (15), we obtain:
(40)
By integrating (40) between 0 and t, we have:
So if the initial condition verifies
then the relation
will be verified for all.
This second stage shows that the solutions are limited for everything. We can conclude that the solutions of the model exist globally in. □
2.2. Desease Free Equilibruim
There will be absence of desease in the population if the proportions, , , , et are nil. Let be (resp.) a desease free equilibruim of the model (15)-(22) (resp. (2)-(9)). The following theorem gives us the existence and the uniqueness of this desease free equilibruim. Given that the models (2)-(9) and (15)-(22) are equivalent, then and are also equivalent.
Theorem 2. The model of HIV-SIDA (2)-(9) or (15)-(22) possesses a unique desease free equilibruim in where
(41)
(42)
and
(43)
(44)
Proof. Let be a desease free equilibruim of the model (15)-(22). There will be absence of desease in the population if. If we substitute these useless values in, , we find that the unique free equilibruim for in from (15) is, the unique free equilibruim for in from (19) is. Consequently the unique desease free equilibruim in is.
Let be desease free equilibruim for the model (2)-(9). By substituting in (11), we obtain:
where
□
2.3. Local Stability of Disease Free Equilibrium
By using the method of Van den Drissche and Watmough, we denote by F the rate of appearance of new infections in compartments of the infectious, and by Vs the rate of transfer of individuals in and out the compartments of the infectious by all other means. Then
The next-generation matrix is defined by:
where, , , , , and are defined by the equations of (45) to (51).
(45)
(46)
(47)
(48)
(49)
(50)
(51)
Proposition 3. The basic reproduction ratio for HIV-SIDA model (15)-(22) is explicitly given by the formula (52) where, and are explicitly defined by equations (45), (48) and (51):
(52)
Theorem 4. The disease free equilibrium of the model (15) is locally-asymptotically stable if and unstable if.
2.4. Global Stability of the Disease Free Equilibrium
We have the following theorem.
Theorem 5. For the system (15)-(22), if then the disease free equilibrium is globally asymptotically stable.
Proof. We begin by rewriting the system (15)-(22) in the form:
(53)
(54)
(55)
(56)
(57)
(58)
(59)
(60)
where and, the disease free equilibrium for the system (53)-(60) corresponds to the point.
Now, let us consider the following function:
If, , , , , and are positives. Consequently the function V is positive, and it nulle at the disease free equilibrium. The derivative of this Lyapunov function V along the trajectories of the ordinary differentiel system is:
We can also write
Algebraic manipulations give
or
If, and are positives, consequently is negative definite along the trajectories. This ends the proof of the theorem. □
2.5. Numerical Simulations
Before closing this section, we verify numerically the theoretical results obtained in subsections 2, 2 and 2 for an initial condition, , , , , , ,. For numerical simulations, the system (15) (22) is discretized with a Runge- Kutta’s method (ODE45). We collect a set of values of biological parameters for the model corresponding to the data on the spread of the HIV-SIDA in two cases:
First case: the disease goes extinct in the population (see Figure 2).
Second case: the disease persists in the population (see Figure 3).
These parameters are obtained in the literature and are summarized in the Table 2.
3. Model without Public Health Sensitization
In this section all sensitization-related parameters and variables are fixed to zero in order to understand the dynamic behavior of the population without public health sensitization campaign.
So, we pose. the model (15)-(22) reduces to:
(61)
(62)
(63)
(64)
where
with. For this sub-model by using the same reasoning in the theorem 1, we demonstrate that for any initial condition in, the system has a unique solution globally defined and which stays in for any time where
(65)
3.1. Local Stability of Disease Free Equilibrium
The desease free equilibruim of the sub-model (61)-(64) is:
By using the method of Van den Drissche and Watmough, we denote by F the rate of appearance of new infections in compartments of the infectious, and by Vs the rate of transfer of individuals in and out the compartments of the infectious by all other means. Then:
the next-generation matrix is defined by:
Proposition 6. The basic reproduction ratio for the sub-model (61)-(64) is given by the formula (66):
(66)
Theorem 7. The disease free equilibrium of the sub-model (61)-(64) is locally-asymptotically stable if and unstable if.
3.2. Global Stability of the Disease Free Equilibrium
Theorem 8. For the system (61)-(64), if then the disease free equilibrium is globally asymptotically stable.
Proof. We begin by rewriting the system (61)-(64) in the form:
(67)
(68)
(69)
(70)
where.
The disease free equilibrium for the system (67)-(70) corresponds to the point.
Now, let us consider the following function:
The function V is positive, and it nulle at the disease free equilibrium. The derivative of this Lyapunov function V along the trajectories of the ordinary differentiel system is:
We can also write
Algebraic manipulations give
or
If, then is negative along the trajectories. This ends the proof of the theorem. □
3.3. Existence and Uniqueness of an Endemic Equilibrium
It is found that an unique endemic equilibrium of (61)-(64) for. Thus, we solve the system:
(71)
(72)
(73)
(74)
From (72), we have:
(75)
From (73), we have:
(76)
From (74) et (76), we have:
(77)
From (75), (76) and (77), we have:
(78)
(76), (77) and (78) in (71) give:
(79)
(79) dans (76) give:
(80)
(79) in (77) give
(81)
If, the system (61)-(64) admits a unique endemic equilibrium.
3.4. Numerical Simulations
Before closing this section, we verify numerically the theoretical results obtained in this section for an initial condition, , ,. For numerical simulations, the system (61)-(64) is discretized with a Runge-Kutta’s method (ODE45). We collect a set of values of biological parameters for the sub-model (61)-(64) corresponding to the data on the spread of the HIV-SIDA in two cases:
First case: the disease goes extinct in the population (see Figure 4).
Second case: the disease persists in the population (see Figure 5).
These parameters are obtained in the literature and are summarized in the Table 3.
4. Evaluation of Impact of Public Health Sensitization
Before using the model (15)-(22) to evaluate the impact of public health sensitization in combatting HIV-AIDS spread in a population, it is instructive to evaluate the behaviour of the model under the worst case scenario (i.e., the case where no public health sensitization is provided in the population). By setting all sensitization related parameters to zero (i.e.,) and using the data in Table 4 and Table 5, simulations of the model (15)-(22) show that in Mali the proportion of infected individuals would reach approximately 0.0686 (let 499550 cas) in 9 years from 2001 (Figure 6(b)). These projections of the model are compatible with the EDSM III projections over the year 2010 which predicted that by the year 2010 in Mali, if measures are not taken to control the epidemic of the HIV-AIDS, about 50000 people could be infected by the virus (see Figure 6).
We resume in Table 4 and Table 5, data of Mali concerning the spread of the HIV-AIDS.
(a) (b)
Figure 4. Dynamics of the system (61)-(64) in case where the disease goes extinct in the population. With the parameters of the Table 3 (first case), we have Rh0 = 0.3077. Figure 4(a) shows the evolution of susceptibles individuals, whereas Figures 4(b)-(d) show the evolution of infected individuals. The system converges towards the desease free equilibruim (1, 0, 0, 0). The simulation was realized with the MATLAB logiciel.
(a) (b)
Figure 5. Dynamics of the system (61)-(64) in case where the disease persists in the populationn. With the parameters of the Table 3 (second case), we have Rh0 = 1.6. Figure 5(a) shows the evolution of susceptibles individuals, whereas Figures 5(b)-(d) show the evolution of infected individuals. The system converges towards the endemic equilibrium (0.6298, 0.038, 0.0195, 0.0111). The simulation was realized with the MATLAB logiciel.
Table 4. Mali epidemiological data for the model (15).
Table 5. Mali demographic data of 2001 used as initial conditions.
(a) (b)
Figure 6. Dynamics of the system (15)-(22) in the nose of the cases, Rh0 = 1.6. Figure 6(a) shows the evolution of susceptibles individuals, whereas Figure 6(b) shows the evolution of infected individuals. We use the parameters of the Table 4. The simulation was realized with the MATLAB logiciel.
We evaluate now the behavior of the model (15)-(22) by considering the impact of Public Health sensitization compaign on the spread of HIV-AIDS in Mali. Using the data in Table 4 and Table 5, simulations of the model (15)-(22) show that in Mali the proportion of infected individuals would reach approximately 0.0137 (soit 100020 cas) in 9 years from 2001 (Figure 6(b)). These projections of the model are compatible with the data found in the literature (Source CIA factbook). According to this source, in Mali the individuals infected by the HIV-AIDS in 2010 were 100000 (see Figure 7).