Optimal Control and Bifurcation Analysis of HIV Model

In this study, a very crucial stage of HIV extinction and invisibility stages are considered and a modified mathematical model is developed to describe the dynamics of infection. Moreover, the basic reproduction number R0 is computed using the next-generation matrix method whereas the stability of disease-free equilibrium is investigated using the eigenvalue matrix stability theory. Furthermore, if R0 ≤ 1, the disease-free equilibrium is stable both locally and globally whereas if R0 > 1, based on the forward bifurcation behavior, the endemic equilibrium is locally and globally asymptotically stable. Particularly, at the critical point R0 = 1, the model exhibits forward bifurcation behavior. On the other hand, the optimal control problem is constructed and Pontryagin's maximum principle is applied to form an optimality system. Further, forward fourth-order Runge–Kutta's method is applied to obtain the solution of state variables whereas Runge–Kutta's fourth-order backward sweep method is applied to obtain solution of adjoint variables. Finally, three control strategies are considered and a cost-effective analysis is performed to identify the better strategies for HIV transmission and progression. In advance, prevention control measure is identified to be the better strategy over treatment control if applied earlier and effectively. Additionally, MATLAB simulations were performed to describe the population's dynamic behavior.


Introduction
Human immunodeficiency virus (HIV) is a virus causing HIV infection and results in economic and life devastating crisis if immediate action is not taken to halt further prevalence of the infection [1][2][3][4][5][6][7][8][9][10][11]. Moreover, this infection has no curing medication, but antiretroviral therapy (ART) or its combination is used for halting further progression of infection by inhibiting the virus in human blood, but if left untreated leads to a sever stage called acquired immunodeficiency syndrome (AIDS) [12][13][14][15][16][17][18][19][20][21][22][23]. HIV infection progress through stages as follows: (i) primary stage (asymptomatic stage): this stage faces human individual where the virus is in the blood cannot be diagnosed with medical instruments; (ii) asymptomatic stage: this stage is symptomless stage of HIV infection but diagnosable with medical test; (iii) symptomatic stage: in this stage, the symptom of HIV infection like tiredness, loss of weight, and extreme loss of water starts to manifest in the life of HIV infected individuals; and (iv) AIDS stage: this is advanced stage of HIV infection where it is difficult for treatment and leads to death soon if special care is not taken [24][25][26][27][28][29][30][31][32]. Modes of HIV transmission are through unsafe sexual practices with HIV-infected person, through contacts of normal blood with HIV infected blood, mother to child through breast or birth time, and any contacts of HIV contains fluids of human's with HIV-negative human fluids [33][34][35]. However, in safe practice, the risk of HIV transmission can be reduced by using principle of abstinence-be faithful-use condom (ABC) [36][37][38]. Optimal control intervention through public health education, using of condom and treatment benefits a lot for continuing human life in safe [1-3, 7, 11, 15, 19, 25-27, 33, 39-44].
A mathematical model is a crucial scientific representation of both biological and physical problems in the form of mathematical equations [45][46][47][48][49][50][51][52]. Moreover, different mathematical models have been developed to describe the dynamics of HIV with optimal controls. Particularly, a SICA mathematical model of HIV with ART effect on HIV patients is described in [9]. But, a mathematical model of HIV infection transmission with undetectable behavior of individuals has not been studied. Hence, in this study, we are motivated to highlight the significance of adherence of ART with optimal control intervention and cost-effectiveness analysis in a modified HIV model. Further, the effect of early starting of ART that leads to an undetectable level of viral load in the blood of humans is considered and the SWIUA mathematical model of HIV with optimal control is developed. Organization of the paper is follows: In Section 2, we formulate the SWIUA model of HIV that describes the stated assumptions, variables, and parameters and in Section 3, a mathematical analysis of the model without control is carried out. Moreover, in Section 4, optimal control problem is discussed. Finally, in Section 5, numerical simulations, results and discussion, and the conclusion are presented.

Model Formulation
In this study, a deterministic mathematical model is formulated to analyze the transmission dynamics of the human immunodeficiency virus through applying prevention and treatment control strategies. The current model is a modification of the SIA model of HIV discussed in [13]. A base model has three compartments of human population and did not explicitly show the effectiveness of ART that reduces the viral load in the blood to the undetectable level for ART users. The present work is done to fill the gap observed in the previous study. The compartments of the current model are described as follows: (i) Susceptible compartment: this compartment is denoted by S and embraces all humans who are free of HIV but have a chance of being infected in the infective environment. All humans in this compartment transfer to the HIV compartment at the transmission rate β, provided that effective contact is done with humans in HIV compartment (ii) HIV untested compartment: this compartment consists of all individuals who are at pre-AIDS stage and not test for HIV. In the formulated mathematical model of HIV, the following assumptions are stated: (i) The total size of population is assumed to be nonconstant (ii) It is assumed that HIV untested individuals get tested for HIV infection before they develop AIDS disease (iii) The total population size at time t is denoted by N ðtÞ is given by Moreover, the notations and description of model variables are given in Table 1, whereas model parameters notations and descriptions are given in Table 2.
Based on flow diagram given in Figure 1 and assumptions considered, without control measures, the dynamics of the populations are represented by the subsequent dynamical system: Computational and Mathematical Methods in Medicine with nonnegative initial conditions are Sð0Þ > 0, Wð0Þ ≥ 0, Ið0Þ ≥ 0, Uð0Þ ≥ 0, Að0Þ ≥ 0.

Invariant Region
+ is a region in which all solutions of the model (2) are bounded provided that initial conditions are bounded. That is, Proof. Consider a total population size, NðtÞ, at time t given by Now, differentiating both sides of equation (4) with respect to time t, we have Reduced to, Integrating both sides of (6) and using comparison theorem [53], we have, Here, from (7) it follows that as t ⟶ ∞,NðtÞ ⟶ λ/μ. That is, for all possibility, the expressions on the right hand side of the inequality (7) either increase to least upper bound or decrease to greatest lower bound as time increase. Hence, NðtÞ ≤ λ/μ. Therefore, the feasible solution set of model (2) is the invariant region Ω, defined as   Ignoring the positive term λ from the preceding equation we get the following inequality.
Solving the forgoing inequality over time interval ½0, t, we get the following inequality.
In the preceding inequality, we know that Hence, the solutions of model (2) are nonnegative. Therefore, from the discussions in subsection it can be concluded that the formulated model is mathematically well-posed and solutions of the model are biologically meaningful.   (2) and (4), we obtain I in terms of S. Finally, using equation (4), setting N = S + W + I + U + A, and expressing all variables in terms of S, we obtain the desired solution. Thus, the computed endemic equilibrium E 1 = ðS * , W * , I * , U * , A * Þ of model (2) is given by

Existence and Uniqueness of Solutions
where 3.5. Basic Reproduction Number. The basic reproduction number is the average number of cases generated by one infected person in wholly susceptible population [45][46][47][48][49][50]. The basic reproductive number R 0 can be determined using the next-generation matrix. In this method, R 0 is defined as the largest eigenvalue of the next-generation matrix [51]. The formulation of this matrix involves classification of all compartments of the model in to two classes: infected and noninfected. Let f be a matrix of newly infected cases and v be a matrix of transition cases in model (2). Consider model (2).
The eigenvalues of next-generation matrix are computed and given by Here, λ 1 is the largest eigenvalue of next-generation matrix. Therefore, the reproduction number, R 0 , is given by 3.6. Bifurcation Analysis. The Hopf bifurcation is studied in [5], and it describes the behavior of the system as system changes its steady state. In this subsection, we verify the possibility of backward and forward bifurcation using the center manifold theory stated in reference [28]. Now, model (2) can be written in the vector form, by renaming the variables as where, Here, from preceding system of nonlinear equation, choosing β 1 as a bifurcation parameter and setting R 0 = 1, we have So that the disease-free equilibrium, E 0 , is locally stable when β 1 < β * 1 and is unstable when β 1 > β * 1 . The linearized matrix of the system around the disease-free equilibrium E 0 and evaluated at β * 1 is given by The eigenvalues of matrix The sign of remaining eigenvalues are determined from characteristic polynomial: where a 0 = 1, We observe that coefficients, a 0 , a 1 , and a 2 , are positive, and all eigenvalues of Jacobian matrix are negative while one eigenvalue is zero. Moreover, let w =

Computational and Mathematical Methods in Medicine
Þ T be the right eigenvector, associated with simple zero eigenvalue and can be obtained by solving equation J E 0 , β * ð Þw = 0. Thus, we have Also Þ be the left eigenvector corresponding to simple zero eigenvalue, obtained by setting and solving vJ In order that the required conditions v:w = 1, the product gives, Again substituting the corresponding components in the preceding equation, we have Thus, the preceding equality is satisfied if we choose Next, we compute all second-order partial derivatives of functions on the right hand of ( * ) as given in the appendix.
Next, we compute bifurcation coefficients a and b: Since a < 0 and b > 0, the model exhibits forward bifurcation at R 0 = 1.
Next, the following procedures given in (Biswas et al., 2020), we compute the bifurcation coefficients a and b, to identify the direction of bifurcation at R 0 = 1. Thus, we have Since all parameters in model (2) are nonnegative and additionally v 2 and w 2 are positive, we conclude a < 0 and b > 0. Thus, according to [5], model (2) exhibits a supercritical (forward) bifurcation, when R 0 crosses the threshold R 0 = 1. That is, there exist locally asymptotically stable endemic equilibrium point E 1 = S * , W * , I * , U * , A * ð Þfor R 0 > 1. Based on the results of the above discussion and [7], the following theorem is stated. Theorem 4. The trans-critical bifurcation of model (2) that occurs at R 0 = 1 is a supercritical (forward) bifurcation. That is, there exists locally asymptotically stable endemic equilib- Remark 5. Theorem 4 shows that if R 0 > 1, then a few inflow of infectious individuals in fully susceptible population can result in persistence of the HIV in the population.

Stability Analysis of Equilibrium Points.
In absence of the infectious disease, model (2) has a unique disease-free steady-state E 0 , and in the presence of disease, model (2) has unique endemic equilibrium E 1 .

Local Stability of Disease-Free Equilibrium
Proof. Consider model (2), so that the Jacobian matrix of the system at DFE is given by The eigenvalues of a preceding Jacobian matrix J is computed and given by Clearly, the first three eigenvalues are negative whereas the fourth and fifth eigenvalues are negative if R 0 < 1.
Further, by a stability analysis of a point using Jacobian matrix, we conclude that the disease-free equilibrium point is locally asymptotically stable if R 0 < 1 and unstable otherwise.  [4]. Accordingly, let X ∈ ℝ 1 denote individuals in uninfected compartment ðSÞ and Y ∈ ℝ 4 denotes individuals in infected compartments W, I, U, A ð Þ. Hence, we write model (2) in the form: Also, the disease-free equilibrium is given by Here, X 0 is the disease-free equilibrium of the foregoing system.
To guarantee global asymptotic stability of the diseasefree equilibrium point, the technique we employed must met the following two conditions H1 and H2.
H1: for dX/dt = HðX, 0Þ, X 0 is globally asymptotically stable H2: GðX, YÞ = PY −ĜðX, YÞ,ĜðX, YÞ ≥ 0 for ðX, YÞ ∈ Ω Here, P = D Y GðX, 0Þ is a Metzler matrix and Ω is a region where solutions of the model are biologically feasible. Theorem 8. The disease-free equilibrium point E 0 of model (2) is globally asymptotically stable in a region Ω if R 0 < 1 as the fourth and fifth eigenvalues are negative if R 0 < 1 and unstable whenever R 0 > 1 provided that the above two conditions H1 and H2 are satisfied, where Ω is a feasible solution region of model (2) Proof. From model (2), we have Putting HðX, 0Þ = 0 and solving, we obtain S = λ/μ. Hence, X 0 = ðλ/μ, 0Þ. Clearly, X 0 is globally asymptotically stable equilibrium point of equation: From infected compartments of model (2), we have In the preceding computation, let so that at disease-free equilibrium it reduces to Again, PY gives Therefore, Since S ≤ S + W + I + U + A, we have ðβ 1 W + β 2 IÞð1 − S /ðS + W + I + U + AÞÞ ≥ 0; the preceding matrix is a nonnegative matrix. Therefore, it can be written aŝ which is a Metzler matrix as all off diagonal elements of a matrix are nonnegative.
Therefore, the disease-free equilibrium E 0 is globally asymptotically stable if R 0 < 1 as the fourth and fifth eigenvalues are negative if R 0 < 1.

Sensitivity
Analysis. Sensitivity analysis is used to determine the sensitivity of the variable with respect to the parameters involved in it [8,13,28]. The normalized forward sensitivity index of a particular variable, R, with respect to a parameter, p, is defined as It is already shown that the explicit expression of R 0 is given by The normalized forward sensitivity indices of R 0 with respect to parameters in it are given by where, based on parametric values given in Table 3, we obtain From Table 3, it can be observed that parameters β 1 and β 2 have positive sensitivity indices and the values of the remaining two parameters ξ, η, and μ get negative sensitivity indices. Further, these parameters affect the value of a 8 Computational and Mathematical Methods in Medicine reproduction number that helps in the analysis of virus extinction or persistence of the disease in the population. On the other hand, an increase of positive parameter value will increase the value of R 0 ; this implies that disease transmission increase with significant amount. Also, an increase in the magnitude of negative parameter value will cause a value of reproduction number to decrease in some amount, which means the disease transmission significantly decreases in some amount.

Optimal Control Problem
In this section, we extend SWIUA model by considering control strategies. The optimal control analysis assists to identify the best control strategies to eradicate or control the disease in the community at a specified period of time [3,15,19]. In the optimal control problem, the following three control measures are used.
(i) Preventive control measure u 1 that protects susceptible population from getting the disease (ii) Treatment control measure u 2 that is used by patients to reduce viral load in the body and slow/ stop the progression of the virus Now, including these control measures u 1 and u 2 in model (2), we get the following optimal control model: with nonnegative initial conditions

Computational and Mathematical Methods in Medicine
To study optimal levels of the controls, we define the Lebesgue measurable control set U as follows: Our goal is to find the optimal controls u * 1 , and u * 2 optimal solutions S * , W * , I * , U * , and A * by fixing the terminal time t f that minimize the objective functional J given by where c 1 , c 2 , w 1 , and w 2 are constants. The expressions 0:5w 1 u 2 1 , 0:5w 2 u 2 2 , and 0:5w 3 u 2 3 are costs associated with controls. The form of cost is quadratic because we assumed it to be nonlinear in nature [3]. Our goal is to minimize the population size of WðtÞ, IðtÞ, and AðtÞ through intervention with control measures u 1 ðtÞ, and u 2 ðtÞ along with the costs associated with them. For two optimal controls u * 1 , and u * 2 , we have where U = fðu 1 , u 2 Þ: 0 ≤ u 1 ≤ 1, 0 ≤ u 2 ≤ 1g, u 1 , u 2 , and u 3 are measurable controls.

Existence of the Optimal Control.
To show the existence of optimal control, we use the approach used in [3]. We have already proved that the HIV model (2) is bounded, so this result can be used to prove the existence of optimal control over finite time interval as applied in [3,52]. To ensure the existence of optimal control, we need to check if the following conditions are satisfied: 10 Computational and Mathematical Methods in Medicine (v) Further, from restriction on control measures, we have Also, considering the preceding inequality, the integrand can be written as where k 1 = min ðw 1 /2, w 2 /2Þ, k 2 = w 2 /2, k = 2. Therefore, there exists optimal control measures u 1 and u 2 that minimize the objective functional J u 1 , u 2 , ð Þ.

The Hamiltonian and Optimality
System. The Pontryagin maximum principle stated the necessary conditions which are satisfied by optimal pair. Hence, by this principle, we obtained the Hamiltonian function ðHÞ defined as where, λ i , i = 1, 2, 3, 4, 5 are the adjoint variables corresponding to state variables S, W, I, U, and A, respectively, and to be determined using Pontryagin's maximal principle for the existence of optimal pairs. Theorem 9. Let S, W, I, U, and A are optimal state solutions with associated optimal control variables u 1 , u 2 , and u 3 for the optimal control model, there exist co-state variables λ 1 , λ 2 , λ 3 , λ 4 , and λ 5 that satisfy With transversality or final time conditions, λ 1 ðt f Þ = λ 2 ð t f Þ = λ 3 ðt f Þ = λ 4 ðt f Þ = λ 5 ðt f Þ = 0 and where H is Hamiltonian function given in ( * ). Furthermore, the optimal controls u * 1 , u * 2 , and u * 3 are

Computational and Mathematical Methods in Medicine
Over the constraints, Proof. Pontryagin's maximum principle gives the standard form of adjoint equation with transversality conditions [40]. Now, differentiating the Hamiltonian function with respect to state variables S, W, I, U and A , respectively, the adjoint equations can be written as Further, the characterization of optimal controls u * 1 , and u * 2 shows that It follows that the optimal solution subject to constraints Therefore, considering the bounds of the control, the optimal control variables are given by

Computational and Mathematical Methods in Medicine
In compact form, the optimal controls can be written as Next, we write the optimality system using a state variables system of equations with initial conditions, a costate variables system of equations with final time conditions and optimal control solution.

Numerical Simulations
Next, we investigate qualitatively the effect of optimal control strategies on the spread of HIV in a population. Hence, we categorized strategies as follows: (i) Strategy 1: only prevention control effort ðu 1 Þ (ii) Strategy 2: only treatment control effort ðu 2 Þ (iii) Strategy 3: both treatment and prevention control efforts ðu 1 and u 2 Þ Solving the optimality system yields the better control. For solving the optimality system, an iterative scheme is used. We begin by using the fourth order Runge-Kutta scheme to solve the state equations with a guess for the controls over simulated time. The adjoint equations are solved using the current iterations solutions of the state equation using the backward fourth-order Runge-Kutta scheme because of the transversality conditions. The controls are then updated by combining the previous controls with the value from the characterizations. In this section, numerical simulations are done to illustrate the analytical results obtained in the above analysis. The initial values for variables of model (2) are Sð0Þ = 6000, Wð0Þ = 200, Ið0Þ = 1000 , Uð0Þ = 300, Að0Þ = 100, c 1 = 20, c 2 = 15, c 3 = 25, w 1 = 1, and w 2 = 1. Moreover, the value and the source of parameters used in the simulations are given in the Table 4.
In Figure 2, the endemic equilibrium changes its stability from unstable to stable and disease-free equilibrium changes its stability from stable to unstable at the bifurcation point R 0 = 1. Moreover, the point R 0 = 1 is the critical point where forward bifurcation behavior of model is exhibited.
13 Computational and Mathematical Methods in Medicine Epidemiologically, interpretation gives the disease persists in the population if R 0 > 1 and extinct if R 0 ≤ 1 In Figure 3, the susceptible population dynamics is simulated considering the prevention and treatment control interventions. The susceptible population size increases due to intervention with prevention control. However, without control, more people get infected and the size of the susceptible population decreases.
In Figure 4, the asymptomatic population dynamics is simulated considering the prevention and treatment control interventions. The asymptomatic population size decrease due to intervention with prevention control. However, without control, more people get infected and the size of the asymptomatic population increases with time.
In Figure 5, the symptomatic HIV population dynamics is simulated considering the prevention and treatment con-trol interventions. The symptomatic HIV population size decreases due to intervention with control measures. However, without control, more people get infected and the size of the HIV population increases with time.
In Figure 6, the undetectable population dynamics is simulated considering the prevention and treatment control interventions. The undetectable population size decreases without intervention with controls due to replication of virus without control. However, with effective adherence to controls, the size of the undetectable population increases with time.
In Figure 7, the AIDS population dynamics is simulated considering the prevention and treatment control interventions. The AIDS population size decreases due to intervention with controls. However, without controls, more people severely attacked with HIV and progress to advanced stage so that the size of the AIDS population increases.
In Figure 8, the susceptible population dynamics is simulated considering only the prevention control. The susceptible population size increase due to intervention with prevention control. However, without prevention control, people get infected and the size of the susceptible population decreases with time.
In Figure 9, the asymptomatic population dynamics is simulated considering the prevention control intervention. The asymptomatic population size decreases due to intervention with prevention control as less number of individuals get infected in the presence of prevention control. However, without control, more people get infected and the size of the asymptomatic population increases with time.
In Figure 10, the symptomatic HIV population dynamics is simulated considering the prevention control intervention. The symptomatic HIV population size increases without intervention with prevention control. However, the prevention control reduces the number of individuals infected whose impact reduces the number of symptomatic HIV population.   14 Computational and Mathematical Methods in Medicine In Figure 11, the undetectable population dynamics is simulated considering the only prevention control intervention. The undetectable population size remains the same without or with presence of prevention control due to absence of treatment control intervention.
In Figure 12, the AIDS population dynamics is simulated considering the prevention control intervention. The AIDS population size increases without treatment control intervention control. However, with control, less people progress to the advanced stage of HIV and the size of the AIDS population decreases with time.
In Figure 13, the susceptible population dynamics is simulated considering only treatment control intervention. The HIV population size increases with treatment control intervention control. However, without treatment control, more people get infected and the size of susceptible population decreases. Figure 14 shows that the absence of the prevention control measure increases the number of HIV-infected individuals in the asymptomatic stage whereas the intervention with only treatment reduces the number HIV-infected individuals in the asymptomatic stage.
In Figure 15, the symptomatic HIV-infected individuals increase the HIV population without control measures, but the early intervention with treatment reduces the number of individuals that get serious with virus.
In Figure 16, the dynamics of undetectable stage HIVinfected individuals are described. The simulation shows that the number of undetectable individuals increases continuously with continuous using of treatment, but in the absence or stopping the medicine causes the number of individuals at the undetectable individuals to become detectable. Figure 17 shows that the number of individuals at AIDS stage increases with time in the absence of the medicine that inhibits the replication of the viruses in the body and prevention control measures. However, early intervention with treatment reduces the number of individuals at the AIDS stage effectively.
In Figure 18, the prevention and the treatment control profiles are simulated. The simulation results indicate that early and continuous application of the prevention and treatment control measures is effective in the controlling the HIV infection.
In Figure 19, the simulation of adjoint variables are performed to show the conditions required in the analysis of optimal control problems. Moreover, at the final time, the adjoint variables yield zero values.
In Figure 20, the simulation of SðtÞ versus WðtÞ is performed; the size of SðtÞ decreases to 2000 as the size of Wð tÞ increases to 1600. However, for size of SðtÞ ≤ 2000 and WðtÞ ≤ 1600, both population sizes increase (decrease) together.
In Figure 21, the simulation of SðtÞ versus IðtÞ is performed. Initially, the size of SðtÞ decrease approximately to 300 as the size of WðtÞ increase to 4700; however, for size    Figure 18: Simulation of prevention and treatment control profiles. 15 Computational and Mathematical Methods in Medicine of SðtÞ ≤ 300 and WðtÞ ≤ 4700, both population sizes increase (decrease) together.
In Figure 22, the simulation of SðtÞ versus AðtÞ is performed. Initially, the size of SðtÞ decrease approximately to 300 as the size of AðtÞ increase to 430; however, for size of SðtÞ ≤ 300 and AðtÞ ≤ 430 both population sizes increase (decrease) together.
In Figure 23, the phase portrait of SðtÞ, WðtÞ, and IðtÞ is performed. The simulation shows the initial occurrence of HIV prevalence increases. However, as time increases, less number of individuals becomes susceptible with the presence of HIV-infected individuals.

Cost-Effectiveness Analysis.
Incremental costeffectiveness ratio (ICER) used to compare the differences between the costs and health outcomes of two alternative intervention strategies that compete for the same resources and is generally described as the additional cost per addi-tional health outcome [54,55].

ICER =
Differences between the costs Difference between health outcomes : ð67Þ In ICER, when comparing two competing intervention strategies incrementally, one intervention should be compared with the next less effective alternative.
The ICER is computed as follows: Comparing the computations, the ICER ð1Þ < ICER ð3Þ.

Computational and Mathematical Methods in Medicine
This implies that strategy (2) is more effective and economical over strategy (6). Therefore, discard strategy 3 and select strategy 1 as a better intervention for optimal control of HIV transmission dynamics. Therefore, based on the results obtained from Tables 5-7, the prevention control strategy is more cost-effective strategy than other control strategies to overcome the burden of HIV on human life.

Result and Discussion.
In this section, we present the result and discuss the results obtained from analytical analysis and numerical simulations of model (2). Particularly, we discuss the impact of ART on HIV transmission as a result of varying the value of transfer rate parameter. In Figure 1, the dynamics of population is described by schematic diagram and the model is developed using it. In Figure 2, we observe that the endemic equilibrium exhibits unstable state behavior if R 0 < 1. But, it exhibits global asymptotic stable state behavior if R 0 > 1. On the other hand, DFE is stable for R 0 < 1 and unstable if R 0 > 1. Moreover, it can be observed that the system shifts from the stable DFE to unstable at the bifurcation point R 0 = 1. That is, there is forward bifurcation at the critical point. Further, in Figures 3-17, the impact of optimal controls is described clearly. Moreover, the presence of optimal control measures increases the number of susceptible individuals. However, the presence of optimal control measures decreases the number of individuals in the W, I, and A classes. In Figure 18, the prevention and the treatment control profiles are simulated. The simulations results indicate early and continuous application of the prevention and treatment control measures are effective in the controlling the HIV infection, whereas in Figure 19 the simulation of adjoint variables are performed to show the conditions required in the analysis of optimal control problems. Moreover, at the final time, the adjoint variables yield zero values. Further, in Figure 20, the simulation of SðtÞ versus WðtÞ is performed; the size of SðtÞ decreases to 2000 as the size of WðtÞ increases to 1600. However, for size of SðtÞ ≤ 2000 and WðtÞ ≤ 1600, both population sizes increase (decrease) together, and in Figure 21, the simulation of SðtÞ versus Iðt Þ is performed. Initially, the size of SðtÞ decreases approximately to 300 as the size of WðtÞ increases to 4700; however, for size of SðtÞ ≤ 300 and WðtÞ ≤ 4700, both population sizes increase (decrease) together; in Figure 22, the simulation of SðtÞ versus AðtÞ is performed. Initially, the size of SðtÞ decreases approximately to 300 as the size of AðtÞ increases to 430; however, for size of SðtÞ ≤ 300 and AðtÞ ≤ 430, both population sizes increase (decrease) together. In Figure 23, the phase portrait of SðtÞ, WðtÞ, and IðtÞ is performed. The simulation shows the initial occurrence of HIV prevalence increases. However, as time increases, less number of individuals becomes susceptible with presence of HIV-infected individuals.
Also, the investigation of cost-effectiveness analysis with possible combination with prevention and treatment control measures for HIV infection shows that early applying of prevention control measures is a better strategy than applying only ART or combined strategy with prevention and treatment control measures.

Conclusion.
In this study, we have modified the HIV model by including an undetectable compartment that emphasizes the impact of using ART properly and continuously in the entire period of HIV infection. Also, the study identified the impact of faulty using ART by varying the parameter values that control right and faulty of using ART. From sensitivity analysis, it is concluded that increasing proper and continuous use of ART help significantly to reduce the number of reproduction number. Moreover, the analytical and numerical simulation results show that the disease-free equilibrium is both locally and globally stable for R 0 ≤ 1: But, the endemic equilibrium is both locally and    17 Computational and Mathematical Methods in Medicine globally stable for R 0 > 1. Furthermore, at the threshold point R 0 = 1, the system experiences a forward branch. Analysis results show that HIV infection can be controlled with proper administration and use of ART. Particularly, with prevention and treatment as the optimal control measures, it reduces the number of individuals living with HIV and the progression of the HIV patient to the AIDS stage. Also, examining cost-benefit analyses of possible combinations of preventive and therapeutic management measures for HIV infection shows that preventive control measure is a better strategy if applied earlier effectively than using ART alone or combined.