Computational and Mathematical Methods in Medicine

Volume 2014 (2014), Article ID 831506, 19 pages

http://dx.doi.org/10.1155/2014/831506

Research Article

Optimal Control of HIV/AIDS in the Workplace in the Presence of Careless Individuals

Baba Seidu¹ and Oluwole D. Makinde²

¹Applied Mathematics Department, Faculty of Mathematical Sciences, University for Development, Navrongo, Ghana

²Faculty of Military Science, Stellenbosch University, Private Bag X2, Saldanha 7395, South Africa

Received 17 April 2014; Revised 29 May 2014; Accepted 31 May 2014; Published 26 June 2014

Academic Editor: Chung-Min Liao

Copyright © 2014 Baba Seidu and Oluwole D. Makinde. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A nonlinear dynamical system is proposed and qualitatively analyzed to study the dynamics of HIV/AIDS in the workplace. The disease-free equilibrium point of the model is shown to be locally asymptotically stable if the basic reproductive number, , is less than unity and the model is shown to exhibit a unique endemic equilibrium when the basic reproductive number is greater than unity. It is shown that, in the absence of recruitment of infectives, the disease is eradicated when , whiles the disease is shown to persist in the presence of recruitment of infected persons. The basic model is extended to include control efforts aimed at reducing infection, irresponsibility, and nonproductivity at the workplace. This leads to an optimal control problem which is qualitatively analyzed using Pontryagin’s Maximum Principle (PMP). Numerical simulation of the resulting optimal control problem is carried out to gain quantitative insights into the implications of the model. The simulation reveals that a multifaceted approach to the fight against the disease is more effective than single control strategies.

1. Introduction

HIV/AIDS is one of the diseases that have claimed and continue to claim the lives of millions of people worldwide. Over the past three decades alone, HIV/AIDS has claimed the lives of more than 25 million people, most of whom were from Sub-Saharan Africa where 1 in every 20 adults is living with HIV. In 2011 alone, for example, about 34 million people, globally, were living with HIV/AIDS, about 23.5 million of them were from Sub-Saharan Africa and about 1.7 million people died from the disease globally [1]. The disease places so many burdens not only on families as some bread winners are lost but also on governments who have to spend millions of dollars in the purchase of antiretroviral drugs and on other intervention schemes. In 2011, for example, there was a total global expenditure of about US$16.8 billion in the fight against HIV/AIDS [1]. In his address of state of the nation this year, the president of Ghana spoke of the government’s commitment to providing about 5 million Dollars to local pharmaceutical companies to help in the production of antiretroviral drugs in the country. It is these effects of the disease that call for continuous research into the prevention and control of the disease.

Mathematical models have played a major role in increasing our understanding of the dynamics of infectious diseases. Several models have been proposed to study the effects of some factors on the transmission dynamics of these infectious diseases including HIV/AIDS and to provide guidelines as to how the spread can be controlled. Among these models include those of Anderson et al. [2] who presented a preliminary study of the transmission dynamics of HIV by proposing a model to study the effects of various factors on the transmission of the disease, Stilianakis et al. [3]. who proposed and gave a detailed analysis of a dynamical model that describes the pathogenesis of HIV, and Tripathi et al. [4] who proposed a model to study the effects of screening of unaware infective on the transmission dynamics of HIV/AIDS. Several other models proposed to study dynamics of HIV/AIDS can be found in ([5–13], and the references therein).

“HIV/AIDS is a major threat to the world of work: it is affecting the most productive segment of the labour force and reducing earnings, and it is imposing huge costs on enterprises in all sectors through declining productivity, increasing labour costs and loss of skills and experience. In addition, HIV/AIDS is affecting fundamental rights at work, particularly with respect to discrimination and stigmatization aimed at workers and people living with and affected by HIV/AIDS. The epidemic and its impact strike hardest at vulnerable groups including women and children, thereby increasing existing gender inequalities and exacerbating the problem of child labour” [14]. Due to the effects of HIV/AIDS on firms, the International Labour Organization sees the field of work as a major stakeholder in the fight against the disease. The ILO envisages a world that sees HIV as a workplace issue like any other disease/sickness. It envisages a world of work that makes efforts to prevent discrimination in any form against people with HIV and also makes efforts to provide healthy work environments through social dialogue, prevention, and care and support for people with HIV. Dixon et al. [15] studied the impact of HIV/AIDS on Africa’s economic development while [16] studied the impact of AIDs on developing economies. Not much research has been done in the study of epidemic models that consider the effect of HIV/AIDS on productivity and how the workplace can contribute to the fight against the disease. Okosun et al. [17] presented a dynamical model that studied the impact of susceptibles and infectives with different levels of productivity on the spread of HIV/AIDS at the workplace. They sought to determine the optimal levels of education, antiretroviral therapy that is required to optimally reduce the spread of the disease and increase productivity. In this paper, we present an extension of the model of Okosun et al. [17] to include susceptibles and infectives with different behaviors towards sex and with varying levels of productivity. Thus, we consider a dynamical system that incorporates the effects Careful-Productive Susceptibles, Careful-Non-Productive Susceptibles, Careless-Productive Susceptibles, Careless-Non-Productive Susceptibles, and similar groups of infectives on the transmission dynamics of HIV/AIDS at the workplace. We study the optimal levels of various intervention strategies needed to optimally reduce the spread of the disease and increase productivity. To do this, we modify our basic model to include various intervention strategies to obtain an optimal control problem which is analyzed qualitatively using the Pontryagin’s Maximum principle. The resulting optimal control problem is also solved numerically to gain more insights into the implications of the interventions. The remainder of the paper is organized as follows. In Section 2, we present the mathematical model describing the dynamics of the disease and some basic properties of the model are also presented. The equilibrium states of the model and some implications are discussed in Section 3. In Section 4, we present a modification of the basic model into an optimal control problem and, finally, we present the results of the numerical simulations of the resulting optimal control problem in Section 5.

2. Formulation of the Model

In this section, we develop a deterministic model that describes the dynamics of HIV/AIDS in a homogeneously mixed workplace of population size . The population is subdivided into nine (9) mutually-exclusive compartments, namely, Careful Productive Susceptibles, ; Careful Non-Productive Susceptibles, ; Careless-Productive Susceptibles, ; Careless-Non-Productive Susceptibles, ; Careful-Productive Infectives, ; Careful-Non-Productive Infectives, ; Careless-Productive Infectives, ; Careless-Non-Productive Infectives, ; and AIDS patients, , so that we have . The schematic diagram of the model is shown in Figure 1.

Figure 1: Flowchart of model (1).

Our model assumes that there is a constant recruitment rate, , into the population with , , , , , , and being the fractions of the respective subpopulations recruited into the population. Susceptible individuals acquire HIV through contact with infected ones with force of infection given by , where is the probability of infection per contact and is the average number of sexual partners per unit time and is a modification parameter due to irresponsibility, which we assume is a factor that increases the chance of an infective transmitting the disease as they may tend to have a negative attitude towards protected sex. Newly infected susceptibles are assumed to be irresponsible and nonproductive as they are often unaware of their HIV status in the early stages and their productivity will reduce due to the infection. This is because, in the asymptomatic phase of the infection, infectives will often experience occasional fevers and general feeling of tiredness and non-feeling-well among others which can negatively impact the productivity. Due to the administration of highly active antiretroviral therapy (HAART), the responsible and irresponsible nonproductive infectives become responsible and irresponsible productives at the rates and , respectively. Responsible and irresponsible nonproductive susceptibles become productive at the rates and , respectively. Careful-Productive Infectives, Careful-Non-Productive Infectives, Careless-Productive Infectives, and Careful-Non-Productive Infectives develop AIDS at the rates , , , and , respectively. There is a positive change in behavior leading to Careless individuals (Productive Susceptibles, Nonproductive Susceptibles, Productive Infectives, and Nonproductive Infectives) becoming careful individuals (Respectively, Productive Susceptibles, Nonproductive Susceptibles, Productive Infectives, and Nonproductive Infectives) at rates , ,, and , respectively. There is a natural death rate of for all individuals in all subgroups and is the disease-induced death rate.

Putting the above formulations and assumption leads to the following set of ordinary differential equations representing the model describing the dynamics of HIV/AIDS at the workplace:

By using , , , , , , , , and and keeping , , for convenience, we have

where

In the next section, some basic facts about model (2) are presented.

2.1. Basic Properties of the Model

We show in this section that model (2) is reasonable both mathematically and biologically. This is achieved via the following theorems.

The model is epidemiologically feasible if the following theorem is true.

Theorem 1. Let . If the initial values of the model are nonnegative (i.e., ), then solutions of model (2) remain positive for all time .

In particular,.

2.1.1. Positive Invariant Region of the Model

Let, . We shall analyze model (2) in the domain

The region, , can be shown to be positively invariant (i.e., solutions in will always remain in ).

Theorem 2. The region is positively invariant for model (2) with initial conditions in .

Proof. Let

Then, adding all equations of model (2) we have

Thus, .

A standard comparison theorem [18] can be used to prove that .

In particular, if , then as required.

This shows that the region is positively invariant and that the dynamics of the model can be sufficiently studied in inside which the model is considered to be epidemiologically and mathematically well posed [19]. This means that all solutions of the model starting in will remain in for all time, .

3. Equilibrium Points of the Model

The model exhibits two equilibrium points, namely, the disease-free equilibrium point, , and the endemic equilibrium point, .

3.1. The Disease-Free Equilibrium

The disease-free equilibrium point exists in the absence of the disease and is given by

where

3.1.1. Basic Reproduction Number

We use the next generation matrix method of [20] to calculate the basic reproduction number,. The transmission and transition matrices are, respectively, given by

with and

van den Driessche and Watmough [20] defined the basic reproduction number, , as the largest eigenvalue of the matrix .

Consider

Using theorem (2) of [20], the following theorem is established.

Theorem 3. The disease-free equilibrium point, , of model (2) is locally asymptotically stable if and unstable if .

The basic reproduction ratio is a threshold quantity that measures the average number of secondary infections caused by a single infected individual introduced into a completely susceptible population over its duration of infectivity [19, 21]. Epidemiologically, Theorem 3 implies that a small influx of infectives will not lead to an epidemic if . The theorem also implies that HIV/AIDs can be eradicated when provided that the initial population sizes are within the region of attraction of the disease-free equilibrium point.

3.1.2. Sensitivity Analysis of Model Parameters

In this section, the relative effects of the parameters that determine are presented. We use the normalized forward sensitivity index defined as follows.

Definition 4. Let. Then the normalized forward sensitivity index of relative to is given by .

Thus, . Due to the complex nature of the resulting expressions, the numerical sensitivity indexes of the remaining parameters are presented in Table 2. These indexes are evaluated using the parameter values in Table 1. Quite a number of these parameter values are used mainly for the simulation purposes to illustrate the kind of response expected for the given parameter values and may not be correct epidemiologically.

Table 1: Model parameter descriptions and values used for simulations.

Parameter
Parameter description
Value
Reference

Rate of recruitment
100 People (Year)⁻¹

Fraction of subpopulations recruited
0.04

Rate at which Careless-Productive Susceptibles become Careful
0.4 (Year)⁻¹

Rate at which Careless-Non-Productive Susceptibles become Careful
0.3 (Year)⁻¹

Rate at which Careless-Productive Infectives become Careful
0.6 (Year)⁻¹

Rate at which Careless-Non-Productive Infectives become Careful
(Year)⁻¹

Rate at which Careful-Non-Productive Susceptibles become Productive.
0.6 (Year)⁻¹

Rate at which Careless-Non-Productive Susceptibles become Productive.
0.4 (Year)⁻¹

Rate at which Careful-Productive Infectives lose their Productivity.
0.4 (Year)⁻¹

Rate at which Careless-Productive Infectives lose their Productivity
0.6 (Year)⁻¹

Contact rate between susceptibles and infectives
0.344 (People)⁻¹
[4]

Modification parameter due to careless behavior towards sex
1.2

Rate of progression of Productive Infectives into AIDs
0.100 (Year)⁻¹
[4]

Rate of progression of Nonproductive Infective into AIDs
0.100 (Year)⁻¹
[4]

Rate of progression of Productive Infective into AIDs
0.100 (Year)⁻¹
[4]

Rate of progression of Productive Infective into AIDs
0.100 (Year)⁻¹
[4]

Natural Death rate
0.020 (Year)⁻¹
[4]

AIDs related death rate
1.000 (Year)⁻¹
[4]

Table 2: Sensitivity indexes of .

Parameter
Parameter description
Sensitivity index

Rate of recruitment
+1.000

Contact rate between Susceptibles and Infectives
+1.000

Average number of sexual partners of an infective per unit time
+1.000

Rate of progression of Productive Infectives into AIDs
−0.460

Rate of progression of Nonproductive Infective into AIDs
−0.015

Rate of progression of Productive Infective into AIDs
−0.033

Rate of progression of Productive Infective into AIDs
−0.082

Natural Death rate
−1.138

Rate at which Careful-Productive Infectives lose their Productivity
+0.000

Rate at which Careless-Productive Infectives lose their Productivity
−0.350

Modification parameter due to careless behavior towards sex
+0.312

Rate at which Careless-Non-Productive Infectives become Careful
+0.180

Rate at which Careless-Productive Infectives become Careful
−0.102

The sensitivity indexes reflect the percentage change in the dependent variable (in this case ) as a result of a percentage change in the independent variable. Thus, a 10% increase (or decrease) in the transmission probability, , leads to a 10% increase (or decrease) in the basic reproduction number, while a 10% increase (or decrease) in the rate of progression of the Productive Infectives into AIDS leads to 4.6% decrease (or increase) in the basic reproduction number. It is observed from Table 2 that the most sensitive parameter is followed by , , and , which are equally sensitive. Thus, these parameters should be those that can be used to control the spread of the disease.

3.2. The Endemic Equilibrium

In the presence of the infection, the system exhibits the endemic equilibrium point, , given by

where

After some algebraic manipulations, it can be shown that satisfies the fourth order polynomial

where

When there is no recruitment of infectives (i.e., ), we have , and, hence, the polynomial has two roots, namely, which corresponds to the disease-free equilibrium and the other being which is positive if and only if . Thus, in the absence of recruitment of infectives, the endemic equilibrium point exists only when and is given by , where

Thus, the following theorem is established.

Theorem 5. In the absence of recruitment of infectives:

(a)if , model (2) has exactly one equilibrium point which is the disease-free equilibrium;(b)if , model (2) has two equilibria, namely, the disease-free equilibrium (8) and the endemic equilibrium point (16), coexisting.

By Theorem 5, a necessary and sufficient condition for eradication of the disease in the absence of recruitment of infectives is that the basic reproduction number, , be less than unity.

Theorem 6. In the presence of recruitment of infectives, model (2) has a unique positive equilibrium irrespective of the sign of. Thus, in the presence of recruitment of infectives, the model does not exhibit backward bifurcation.

Since all model parameters are nonnegative, then clearly and, hence, the discriminant of the quadratic equation, , is positive. By the Descartes rule of signs, the polynomial has two real roots of opposite signs. Hence, the model has a unique endemic (positive) equilibrium irrespective of the sign of .

4. Extended Model with Controls

In this section, an optimal control problem is formulated by incorporating four intervention strategies into our basic model (2). The following interventions are incorporated into the basic model:

(i) is the control effort aimed at reducing the infection of susceptible individuals;(ii) is the control effort aimed at treating infected individuals;(iii) is the control effort aimed at changing behavior. That is, is the control effort aimed at making Careless Susceptibles (both Productive and Nonproductive) and Infectives (both Productive and Nonproductive) Careful;(iv) is the control effort aimed at reducing nonproductivity at the workplace.

Thus, the basic model becomes

Our main aim in developing this extended model is to seek optimal levels of the intervention strategies needed to minimize the number of nonproductive workers and the cost of implementing the control strategies. We choose a functional given by

where the s are positive weights which measure relative costs of implementing the respective intervention strategies over the period , whilst the terms measure the cost of the intervention strategies. We chose a quadratic cost functional in line with several other literatures on models of epidemic control [22–25]. Thus, we seek an optimal control quadruple such that

where such that are measurable with ; is the set of admissible controls.

Pontryagin’s Maximum Principle [26] provides the necessary condition for optimality of the controls. Using this principle, (17) and (18) are converted into a problem of minimizing, with respect to the controls s, the Hamiltonian given by

The s, are the adjoint variables or costate variables which determine the adjoint system, which together with the state system (17) describes the optimality system.

Pontryagin’s Maximum principle [26] and the existence result for optimal control from [27] can be used to obtain the following proposition.

Proposition 7. The optimal control 4-tuple minimizes the functional if there exist adjoint variables , that satisfy the adjoint system given by

With transversality conditions , , where

Further more

Proof. We obtain the existence of the optimal controls from [27, Corollary 4.1] due to the convexity of the integrand of the functional with respect to the quadruple , a prior boundedness of the state solutions, and the Lipschitz property of the state system with respect to the state variables. Using Pontryagin’s Maximum Principle, the adjoint or costate equations (21) are obtained by differentiating the Hamiltonian partially with respect to the state variables. Thus, we have

Since the Hamiltonian is minimized at the optimal controls, the optimality conditions at are met. These optimality conditions can be used to obtain expressions for . By standard control arguments involving the bounds on the controls, (24) is obtained, concluding the proof.

5. Numerical Simulations

5.1. Methodology

The solution of the optimal control problem is obtained by solving the optimality system which consists of the state and adjoint systems (17) and (21), respectively. For computational illustration, the values of parameters in Table 1 were employed and the solution is obtained by using the following iterative scheme.

Step 1. Make a guess of the controls.

Step 2. Use the values of the controls together with the initial conditions to solve the state equations, using a forward numerical scheme.

Step 3. Using the current solution of the state system together with the transversality conditions, solve the adjoint equations using a backward numerical scheme. We use a backward scheme for the costate system because the transversality conditions are final time conditions.

Step 4. Update the controls using the characterizations in (24).

Step 5. Repeat Steps 2 to 4 until the values of the unknowns at the current iteration are very close to those of the previous iteration [28].

We note here that human resource departments could only be concerned with reducing nonproductivity or seek to reduce the effect of HIV or combine both efforts. To compare the effects of these options, we consider the following combinations of the controls.

 Strategy 1: Implementing all controls (i.e., ) Strategy 2: Implementing the controls aimed at reducing infection and treating infected individuals (i.e., , ) Strategy 3: Implementing only the control effort aimed at reducing Nonproductivity (i.e., , ) Strategy  4: Implementing only and  Strategy 5: Implementing only .

5.2. Results

In this section, we present the results of the numerical simulation of our optimal control problem by discussing the implications implementing the five intervention schemes above.

We examine the effects of applying the intervention schemes in each of the strategies. Thus, we aim to determine the optimal levels of the controls that will minimize the objective functional . To observe the effects of the intervention strategies, we plot results from simulation of the uncontrolled model (2) and that from the controlled one (17) together in Figures 2 to 6. It is observed in Figures 2(a)–2(d), 3(a)–3(d), 4(a)–4(d), 5(a)–5(d) and 6(a)–6(d) that the number of susceptives remains higher for the controlled problem than for the uncontrolled problem. That means that each of the intervention strategies will lead to saving more people from being infected. It is also observed in Figures 2(e)–2(h), 3(e)–3(h), 4(e)–4(h), 5(e)–5(h) and 6(e)–6(h) that implementing the controls in each strategy will lead to a reduction in the number of people infected with the disease and also reduces the number of Nonproductive individuals.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 2: Simulations of Basic model (2) and the Optimal Control Problem (17) showing the effect of implementing all the intervention strategies on the dynamics of HIV/AIDS transmission.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 3: Simulations of the Optimal Control Problem (17) showing the effect of control Strategy 2 on the dynamics of HIV/AIDS transmission.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 4: Simulations of the Optimal Control Problem (17) showing the effect of control Strategy 3 on the dynamics of HIV/AIDS transmission.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 5: Simulations of the Optimal Control Problem (17) showing the effect of control Strategy 4 on the dynamics of HIV/AIDS transmission.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 6: Simulations of the Optimal Control Problem (17) showing the effect of control Strategy 5 on the dynamics of HIV/AIDS transmission.

To compare the various strategies, we also plot the results of all the strategies on same graphs as in Figure 7. It is observed from the graphs in Figure 7 that the strategy that involves implementing all the controls leads to higher susceptible populations and lower infectives populations. This implies that the fight against HIV/AIDS should be multifaceted in order to achieve maximum benefits.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 7: Simulations of the Optimal Control Problem (17) showing the effect of control Strategy 1 on the dynamics of HIV/AIDS transmission.

6. Conclusion

In this paper, a nonlinear dynamical model has been proposed to study the dynamics of HIV/AIDS in the workplace. The model assumes that there is no discrimination against people living with HIV and that, thus, allows for recruitment of both susceptible and infected individuals by the human resource department. Disease-free and endemic equilibrium states are shown to exist for certain parameter values of the model.

It is shown that the model cannot have a disease-free equilibrium point when infectives are recruited, which is in agreement with [29, 30]. A sensitivity analysis of the basic reproduction number indicates that rate of recruitment, death rate, transmission probability, and number of sexual partners of infected persons are the most sensitive parameters that can be used to control the spread of the disease. Thus, these parameters are those that should be targeted most by policymakers in the fight against the disease. Due to the International Labour Organization’s campaign for no discrimination on the basis of ones’ HIV status at the workplace, using the recruitment rate might be compromised, but using the other parameters can still be of immense help. The model is extended to an optimal control problem by incorporating time-varying controls into the basic model and the conditions for optimality are derived using the Pontryagin’s Maximum Principle [26]. Finally, numerical simulations of the resulting control problem are carried out to determine the effectiveness of various combinations of the controls. It is revealed from the simulation of the control problem that the strategy that employs all the control efforts is most effective in the fight against the disease. Thus, there is the need for a multifaceted approach in the fight against the spread of HIV/AIDS.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

Joint United Nations Programme on HIV/AIDS (UNAIDS), Report on the Global HIV/AIDS Epidemic, Joint United Nations Programme on HIV/AIDS, 2012.

R. M. Anderson, G. F. Medley, R. M. May, and A. M. Johnson, “A preliminary study of the transmission dynamics of the human immunodeficiency virus (HIV), the causative agent of AIDS,” IMA Journal of Mathematics Applied in Medicine and Biology, vol. 3, no. 4, pp. 229–263, 1986.

N. I. Stilianakis, K. Dietz, and D. Schenzle, “Analysis of a model for the pathogenesis of AIDS,” Mathematical Biosciences, vol. 145, no. 1, pp. 27–46, 1997.

A. Tripathi, R. Naresh, and D. Sharma, “Modelling the effect of screening of unaware infectives on the spread of HIV infection,” Applied Mathematics and Computation, vol. 184, no. 2, pp. 1053–1068, 2007.

B. D. Aggarwala, “On two ode models for HIV/AIDS development in Canada and a logistic seir model,” Far East Journal of Applied Mathematics, vol. 6, no. 1, pp. 25–70, 2002.

J. T. Bertrand, K. O'Reilly, J. Denison, R. Anhang, and M. Sweat, “Systematic review of the effectiveness of mass communication programs to change HIV/AIDS-related behaviors in developing countries,” Health Education Research, vol. 21, no. 4, pp. 567–597, 2006.

J. Chin, “Current and future dimensions of the HIV/AIDS pandemic in women and children,” The Lancet, vol. 336, no. 8709, pp. 221–224, 1990.

M. I. Daabo, O. D. Makinde, and B. Seidu, “Modelling the spread of HIV/AIDS epidemic in the presence of irresponsible infectives,” African Journal of Biotechnology, vol. 11, no. 51, pp. 11287–11295, 2012.

D. Mohammed Ibrahim and B. Seidu, “Modelling the effect of irresponsible infective immigrants on the transmission dynamics of HIV/AIDS,” Advances in Applied Mathematical Biosciences, vol. 3, pp. 31–40, 2012.

P. Essunger and A. S. Perelson, “Modeling HIV infection of CD4⁺ t-cell subpopulations,” Journal of Theoretical Biology, vol. 170, no. 4, pp. 367–391, 1994.

B. Flugentius, J. Y. T. Mugisha, and L. S. Luboobi, “An HIV/AIDS model with variable force of infection and its application to the epidemic in uganda,” The American Journal of Applied Sciences, vol. 2, pp. 1274–1278, 2005.

A. B. Gumel, P. N. Shivakumar, and B. M. Sahai, “A mathematical model for the dynamics of HIV-1 during the typical course of infection,” Nonlinear Analysis: Theory, Methods and Applications, vol. 47, no. 3, pp. 1773–1783, 2001.

M. Coffee, M. N. Lurie, and G. P. Garnett, “Modelling the impact of migration on the HIV epidemic in South Africa,” AIDS, vol. 21, no. 3, pp. 343–350, 2007.

ILO, “ILO code of practice on HIV/AIDS and the world of work,” Tech. Rep., International Labour Organization, Geneva, Switzerland, 2001.

S. Dixon, S. McDonald, and J. Roberts, “The impact of HIV and AIDS on Africa's economic development,” British Medical Journal, vol. 324, no. 7331, pp. 232–234, 2002.

F. Fleck, “Developing economies shrink as AIDS reduces workforce,” BMJ, vol. 329, article 129, 2004.

K. O. Okosun, O. D. Makinde, and I. Takaidza, “Analysis of recruitment and industrial human resources management for optimal productivity in the presence of the HIV/AIDS epidemic,” Journal of Biological Physics, vol. 39, no. 1, pp. 99–121, 2013.

V. Lakshmikanthan, S. Leela, and A. A. A. Martyniuk, Stability Analysis of Nonlinear Systems, CRC Press, 1989.

H. W. Hethcote, “The mathematics of infectious diseases,” SIAM Review, vol. 42, no. 4, pp. 599–653, 2000.

P. van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002.

R. Anderson, Infectious Diseases of Humans: Dynamics and Control, Oxford University Press, Oxford, UK, 1991.

B. M. Adams, H. T. Banks, H. Kwon, and H. T. Tran, “Dynamic multidrug therapies for HIV: optimal and STI control approaches,” Mathematical Biosciences and Engineering, vol. 1, no. 2, pp. 223–241, 2004.

E. Jung, S. Iwami, Y. Takeuchi, and T.-C. Jo, “Optimal control strategy for prevention of avian influenza pandemic,” Journal of Theoretical Biology, vol. 260, no. 2, pp. 220–229, 2009.

K. O. Okosun, O. D. Makinde, and G. J. Abiodun, “Transmission Dynamics of HIV/AIDS with optimal control in the presence of carefree suscptibles and treatment,” in BIOMAT Book Series, chapter A, pp. 131–152, World Scientific Publishing, 2011.

K. O. Okosun, O. D. Makinde, and I. Takaidza, “Impact of opti mal control on the treatment of HIV/AIDS and screening of unaware infectives,” Applied Mathematical Modelling, vol. 37, no. 6, pp. 3802–3820, 2013.

L. S. Pontryagin, The Mathematical Theory of Optimal Processes, vol. 4, CRC Press, 1962.

W. H. Fleming and R. W. Rishel, Deterministic and Stochastic Optimal Control, vol. 1, Springer, New York, NY, USA, 1975.

S. M. Lenhart and J. T. Workman, Optimal Control Applied to Biological Models, CRC Press, 2007.

F. Brauer and P. van den Driessche, “Models for transmission of disease with immigration of infectives,” Mathematical Biosciences, vol. 171, no. 2, pp. 143–154, 2001.

O. D. Makinde and K. O. Okosun, “Impact of chemo-therapy on optimal control of malaria disease with infected immigrants,” BioSystems, vol. 104, no. 1, pp. 32–41, 2011.