A Computational Approach to a Model for HIV and the Immune System Interaction

: This study deals with the numerical solution of the human immunodeﬁciency virus (HIV) infection model, which is a signiﬁcant problem for global public health. Acquired immunodeﬁciency syndrome (AIDS) is a communicable disease, and HIV is the causative agent for AIDS, which damages the ability of the body to ﬁght against disease and easily usual innocuous infections attack the body. On entering the body, HIV infects a large amount of CD4+ T-cells and disturbs the supply rate of these cells from the thymus. Herein, we consider the model with variable source terms in which the production of these cells is a monotonically decreasing function of viral load. Based on the reproduction number, we describe the stability of free equilibrium. The continuous Galerkin–Petrov method, in particular the cGP(2)-method, is implemented to determine the numerical solutions of the model. The inﬂuence of different parameters on the population dynamics of healthy/infected CD4+ T-cells and free HIV particles are examined, and the results are presented graphically. On the other hand, the model is solved using the fourth-order Runge–Kutta method, and brieﬂy, the RK4-method, and the results of the proposed schemes are compared with those obtained from other classical schemes such as the Bessel collocation method (BCM), Laplace Adomian decomposition method (LADM), perturbation iteration algorithm (PIA), modiﬁed variational iteration method (MVIM), differential transform method (DTM), and exponential Galerkin method (EGM), numerically. Furthermore, absolute errors relative to the RK4 method are computed to describe the accuracy of the proposed scheme. It is presented that the cGP(2)-method gains accurate results at larger time step sizes in comparison with the results of the aforementioned methods. The numerical and graphical comparison reveals that the proposed scheme yields more accurate results relative to other traditional schemes from the literature.


Introduction
Dynamical systems are well known for their various applications, such as in a population growth model, biomedical, and engineering. For example, in the field of population dynamics, one can study the changes in population with respect to time and make a longterm prediction about population growth for the future. Ding and Ye [1] used a fractional order model of HIV infection. They showed a non-negative solution of that model along with stability analysis. They also expressed numerical simulation and displayed the result of the HIV mathematical model. Burton and Mascola [2] analyzed that AIDS is a result of HIV that minimizes the capability of body scraps beside the disorder and makes it more suitable for natural inoffensive infection. HIV infects a huge quantity of T-cells and divides rapidly while inflowing into the host. Throughout the earliest phases of the disease, blood consists of the soaring power of HIV harmful particles, which extend throughout the body as the blood circulates. These microorganisms swell into liquid substances of the body, e.g., blood, tears, urine, etc. Most HIV-infected individuals will develop AIDS 10 to 15 years after infection, while several people stay fit for more than this period without ailment. Samanta [3] identified that HIV has a stretched evolution and transmissible phase. In the nonexistence of antiretroviral treatment, the usual period of evolution from HIV to AIDS is 9 to 10 years. The average endurance duration following AIDS is 9 months and 20 days.
Furthermore, the ratio of medical infection evolution depends broadly on personalities, from 14 days until 20 years. Several issues upset the ratio of development. These causes embrace the host's capacity to protect in opposition to HIV. Humble entrance to physical care and the presence of associated infections may also inspire quicker illness development. Krischner et al. [4] discussed treating HIV infection and certain chemotherapies have been tested. The most commonly used drugs are AZT, DDI, DDC, and D4T. All these drugs are inhibiters, which work like reverse transcriptase inhibitors. The reverse transcriptase inhibitors reduce the growth of the virus and finally block the infection.
Many researchers have worked in this field to conduct studies and analyze research to control this disease. Chum et al. [5] addressed the rule of sexually transferred infection in HIV-1. Sun and Min [6] studied a modified HIV infection model simulation. Khalil et al. [7] analyzed HIV infection with the drug therapy effect. They also discussed the generalized Euler method (GEM) and the numerical solution of the HIV fractional-order model. Wang et al. [8] examined the HIV pathogenesis model along with its cure rate and logistic growth of healthy and infected CD4+ T-cells. Osman and Abdurrahman [9] studied the stability of the delayed HIV model and investigated its treatment and transmission. They also computed the basic reproduction number R 0 and discussed its stability based on the reproduction number. Wan et al. [10] discussed the behavior of the viral dynamical model and its solution. Wang et al. [11] addressed the global dynamics of HIV infection of CD4+ T-cells and analyzed the mathematical model. Srivastava et al. [12] presented drug therapy for HIV infection. Culshaw et al. [13] demonstrated the differential equation of HIV infection and discussed the modeling and dynamics of CD4+ T-cells.
HIV/AIDS is a universal health challenge; over 70 million people worldwide have been infected with this disease, from which 35 million people have passed away, and 36.7 million patients are still living with this infection [14]. Nelson and Perelson [15] worked on the model of HIV infection and its mathematical analysis based on kinetic parameters. They analyzed intracellular delays, antiretroviral therapy, and dynamics of infected and uninfected CD4+ T-cells. They also illustrated that when the drug effect is less than the loss rate of productively infected cells, the approximated value is ideal and depicts the results of the model. Perelson et al. [16] examined the interaction of HIV and CD4+ T-cells. They discussed the infected, active infected, uninfected, and free virus T-cells. They explained the characteristics of HIV infection in a straightforward manner. They considered two types of HIV models. The first is with a constant source for the growth of CD4+ T-cells, and the second is a decrease in the relative viral load of CD4+ T-cells. Ronga et al. [17] discussed a mathematical model to examine the method of drug resistance during therapy and its appearance. They indicated that drug-resistant viruses contributed significantly to their stability. When antiretroviral treatments are available, both strain and wild-type reproduction are close to each other. Duffin and Tullis [18] worked on a mathematical model for HIV infection and AIDS to understand these diseases. They suggested a mathematical model of targeted cells; for this purpose, they evaluated two models and compared the results with actual observations. Song and Chang [19] described the time between the infection of T-cells and the secretion of viral elements. In addition, they demonstrated the condition under which the model is stable and the model's existence in that solution. Mechee and Haitha [20] investigated the application of Lie symmetry for an HIV-infected model that deals with the initial value of a nonlinear differential equation. Based on Lie symmetry approach, Mechee and Haitha [20] attempted to identify uninfected T-cells in the host body. Zhou et al. [21] and Leenheer et al. [22] studied a model for HIV infection with its cure rate. They proved that when the basic reproduction number, R 0 < 1, there are no infections in the host cells, were as when R 0 > 1, the host cell is infected by HIV infection. They also found conditions for which the system is stable or unstable. Sarivastav and Chandra [23] proposed the dynamics of HIV and CD4+ T-cells in primary infection, also modeling for HIV infection. They analyzed stability during the infection state of the host cells of the body. They further discussed the local stability of the model and showed the results numerically. They determined the basic reproduction number R0 and discussed that when R 0 , the host cells are free of disease; when R 0 > 1, the host cells are infected with the disease. Liu and Li [24] formulated an HIV infection model that depends on the age of infection, elapsed age of disease, and antiretroviral treatment. The model is concerned with two classes, infective (concentration of infection at birth) and AIDS (concentration of disease without birth), whose results give an approximate interval of optimal pulse and proportion of impulsivity. Ho et al. [25] discussed the rapid turnover of plasma virion and HIV-1 infection. They also demonstrated the treatment of individuals infected with HIV-1 infection, which causes AIDS. Perelson et al. [26] discussed the decay characteristics of HIV-1-infected compartments during combination therapy. Kuang [27,28] proposed the delay differential equation with applications in population dynamics. Ongun [29] implemented the Laplace Adomian decomposition method for solving a model for HIV infection of CD4+ T-cells. Yuzbasi [30] employed the numerical approach to solve the model for HIV infection of CD4+ T-cells. Khalid et al. [31] discussed the numerical solution of a model for HIV infection of CD4+ T-Cells. Merdan et al. [32] described the numerical solution of the model for HIV infection of CD4+ T-Cells. Attaullah et al. [33] studied the transmission and dynamical behavior of an HIV/AIDS epidemic model with a cure rate based on mathematical model. Ogunlaran and Noutchie [34] considered the HIV-infected model with two variables, uninfected CD4+ T-cells and incident term of the free virion. Ogunlaran and Noutchie [34] aimed to increase the concentration of uninfected CD4+ T-cells using minimal drug therapies and to stop the reproduction of infected cells. Boukari et al. [35] analyzed a discrete HIV infection model with a time delay and expressed the global stability of the model numerically. They used the backward Euler method and proved that E 0 is globally asymptotically stable when R 0 < 1. Li and Xiao [36] investigated the global dynamics of a virus' immune system in order to establish the HIV load and structured treatment outages. They also addressed the global dynamics of the HIV elimination and infection cell growth model. Espindola et al. [37] investigated macrophages and their role in HIV infection. They also discussed how highly active antiretroviral therapy (HAART) affected HIV infection. Kinner et al. [38] examined the incidence of HIV, hepatitis B, and hepatitis C in older and younger adults and deduced that the incidence is lower in younger adults than in older adults. Angulo et al. [39] demonstrated that the main path of HIV-1 infection is transmitted from a mother to her child. They also found polymorphisms in human leukocyte antigen class-B (HLA-B) concerned with HIV-1 infection. Theys et al. [40] studied HIV-1 impact on host cells and their transmission. They also found a link between the evolution of the host cell and the fitness of the host cell. Hallberge et al. [41] established a developed stage of knowledge on the significance of HIV revelation between partners. They discussed that most HIV infection is transferred from one infected person (female) to another person (male) and the status of both are the main factor of this disease. Ransome et al. [42] analyzed the spread, cure, and prevention of HIV infection in social relationships. They realized that social capital is an important factor in HIV transmission from one person to another. Naidoo et al. [43] studied the care of tuberculosis (TB) and their class in those people that are already infected with HIV. Omondi et al. [44] considered a mathematical model of HIV infection and investigated the transmission between two kinds of different ages. They also showed that males are less infected than their female partners by this infection. Duro et al. [45] illustrated the CD4+ T-cell monitoring in HIV-infected people with the help of CD4+ T-cell counts. They also found the possibility of CD4+ T-cells being maintained during viral suppression by using the Kaplan-Meier technique. Mbogo et al. [46] considered the model of HIV infection with the help of the Stochastic approach and the probability of HIV when clear, which depends on drug cure rate and intracellular delay length, both of which play a vital role in HIV progression. Ghoreishi et al. [47] utilized the homotopy analysis method (HAM) to determine the solution of the HIV-infected model in the form of an infinite series. They used auxiliary parameters to adjust and control the convergence portion of the infinite series. Elaiw [48] presented the HIV model along with two categories of targeted cells, T-cells and macrophages, with an infection rate. They used Lyapunov and Lasalle principles to recognize the global stability of the infected and uninfected state. They computed the basic reproduction number "R 0 " and analyzed that when "R 0 < 1" the uninfected state is globally asymptotically stable, and when "R 0 > 1" the infected state is globally asymptotically stable. Ali et al. [49] computed the solution of the HIVinfected model by using the Adomian decomposition method (ADM) that illustrates the solution of ODE's in terms of infinite series components. Yuzbasi and Karacayir [50] considered a model of HIV infection and determined the solution of the model by using the exponential Galerkin method (EGM). They used a technique of residual correction. The purpose of this technique is to reduce the error of the solution. They also showed his result numerically and compared them with numerous existing methods. Kirschner [51] studied HIV immunological dynamics employing mathematics. In the Chemotherapy of AIDS, Webb et al. [52] described the mathematical model for HIV treatment approach. Attaullah and Sohaib [53] implemented two numerical schemes, namely, continuous Galerkin-Petrov (cGP(2)) and Legendre wavelet collocation method (LWCM), for the approximate solution of the mathematical model which describes the behavior of CD4+ T-cells, infected CD4+ T-cells, and free HIV virus particles after HIV infection. They presented and analyzed the effect of constant and different variable source terms (depending on the viral load) used for the supply of new CD4+ T-cells from the thymus on the dynamics of CD4+ T-cells, infected CD4+ T-cells, and free HIV virus. Furthermore, they also solved the model using the fourth-order Runge-Kutta (RK4) method. They highlight the accuracy and efficiency of the proposed schemes with the other traditional schemes.

Main Objectives
The main contribution of this paper is to present the implementation and analysis of the cGP(2)-method [53] on the HIV-infection model proposed by Yuzbasi and Karacayir [50] and compare the results with those obtained from other conventional methods existing in the literature. Comparisons verify that the new findings align well with the existing solutions. Since the method mentioned earlier is based on the finite element method, the present scheme can be utilized as an adaptive scheme. Moreover, this variational type of discretization has several advantages over the standard schemes. In the existing literature, the majority of the researchers presumed that the supply rate of new CD4+ T cells from the thymus is constant when assessing the dynamics of HIV. However, HIV is capable of infecting these cells in the thymus, and variable phenomena have been observed as contrasted to constant versions. Therefore, we extended the model considered in [50] by introducing a variable source term, a monotonically decreasing function depending on the concentration of viral load. To understand the stability, we determined the basic reproduction number of the model. The model is essential for mathematically simulating HIV infection of CD4+ T cells. This will be employed to study CD4+ T-cell population dynamics in the involvement and complete lack of HIV, which will assist in the perception of diagnostic and therapeutic AIDS symptoms as well as diminishing the disease. This will be a helpful addition to the current literature on biomathematics. In the remainder of the present study, healthy cells, infected cells, and the virus will be used to mean healthy, infected CD4+ T-cells, and free HIV particles, respectively.

Mathematical Formulation of HIV Model
This section is concerned with the HIV infection model described by Yuzbasi and Karacayir [50], in which the total population is divided into three compartments, i.e., con-centration of healthy T-cells T(t), concentration of infected T-cells I(t), and concentration of free HIV V(t) at time t. The description of the model is as follows: Initial conditions for this system are T (0) = T 0 , I(0) = I 0 , and V(0) = V 0 , such that T(t) indicates the concentration of healthy T-cells, I(t) indicates the population of infected T-cells, and V(t) indicates the dynamics of free HIV at time t.

Modified Formulation of HIV Model
This section is concerned with the extension of the model (Yuzbasi and Karacayir [50]) by introducing a variable source term depending on the viral load instead of using a constant source for the supply of new cells from the thymus, i.e., "0.5η + 5η 1+V(t) " [51,52] where "η" is the supply rate of healthy cells. The modified model is as follows: The initial conditions and parameters involved in the model are explained in Table 1.

Uninfected Steady State
In this condition, viruses are not present in the human cells and in other living things. In this situation, "E 0 " is stable and "E" is unstable in the uninfected steady state. In the uninfected steady condition, we consider that T = T 0 , I = 0, and V = 0. Then, we discover "T 0 " from the uninfected condition from System (1), which after some simplification, is as follows: then, using the quadratic formula to find the value of "T 0 " we have: , and C = −5.5η. Using the above we get:

Infected Steady State
In an infected state, viruses are in human cells and in other living things. In an infected steady state, we consider that T = T, I = I, and V = V in System (1). To obtain the value of T, I, and V, we get the following after simplification:

Reproduction Number
In the cure rate of infection, the virus can easily be controlled when the value of the reproduction number is smaller because the viral infection speed is less. However, the disease cannot be controlled easily without a cure rate because the viral infection speed is larger, which shows a faster speed of the virus. We determined that the reproduction number R 0 , and observe that the local and global dynamics are absolutely resolute by the significance of R 0 . When R 0 ≤ 1, then the disease-free equilibrium "E 0 " is locally and globally asymptotically stable, while R 0 > 1, then the disease-free equilibrium "E 0 " becomes locally and globally asymptotically unstable. The reproduction number is

Jacobian Matrix
The Jacobian matrix corresponding to System (7) about E(T, I, V) is provided by:
Proof. The Jacobian matrix for the given System (1) about E 0 (M, 0, 0) is as follows: the characteristic equation of the above matrix is det(λI − J(E 0 )) = 0, where "I" is the unit matrix. Then, to expand this matrix with respect to the 1st column we obtained It is clear from the above equation that one root of this equation is λ = −b. Then, to determine the other roots, we consider the equation as: Hence, by Hurwitz Criterion [28] all roots of the given equation have a negative real part.
Proof. E 0 (M, 0, 0), given from the Jacobian matrix of System (7), is presented below: the characteristic equation of the above matrix as det(λI − J(E 0 )) = 0, and "I" is the unit matrix. Then, we expand this matrix with respect to the 1st column and obtained: It is clear from the above equation that one root of this equation is λ = −b. Then, to determine the other roots, we consider the equation: Hence by Hurwitz-Criterion [28] all roots of this scheme have a positive real part.

Theorem 3.
For any positive solution (T(t), taking the derivative of Equation (5), with the solution of System (1), to obtain: . Then, N 1 > 0 exists, depending on the parameter of the given system. Therefore, L(t) < N 1 , for all t. Then, T(t) and I(t) are both bounded above. Then, we also have V(t) bounded above from the 3rd equation of System (1), and the maximum of this is N. The proof is complete.
Then, the equilibrium E T, I, V is locally asymptotically stable.
Proof. For equilibrium E T, I, V , the given system after simplification reduces to where, To find the second part of this Theorem 4, we have By Routh-Hurwitz criterion [28], i.e., Then, E T, I, V is locally asymptotically stable (LAS) .

Theorem 5. System (1) is a competitive system.
Proof. Consider a matrix "R" and Jacobian matrix of System (1), such that For some partial order defined as K = (T, I, V) ∈ R 3 : T ≤ 0, I ≥ 0, V ≥ 0. After simplification, we obtained: This proof is complete.

The Continuous Galerkin-Petrov Method
Numerical methods are widely used to simulate complex real-world problems. Jiwari [54][55][56] proposed the local radial basis function-finite difference based algorithms for the singularly perturbed Burgers' model. Mittal et al. [57] suggested a cubic B-spline quasiinterpolation algorithm to capture the pattern formation of coupled reaction-diffusion models. Pandit [58] discussed the local radial basis functions and scale-3 Haar wavelet operational matrices based on numerical algorithms for a generalized regularized long wave model. Mittal et al. [59] considered a new scale-3 haar wavelet algorithm for numerical simulation of second-order ordinary differential equations. Nowadays, the cGP-method has been successfully employed to solve many types of non-linear problems in science and engineering, for example [27,[33][34][35][36][37][38][39]. In this paper, we applied this approach to the HIV infection model. The system of ODEs for HIV Model (1) can be considered as: where where T(0), I(0) and V(0) are the initial conditions given in Table 1.
In order to find the approximate solution of System (1), we partitioned the time interval The symbol τ = t m − t m−1 is used to represent the maximum time step size. For the derivation of the cGP-method, the system of equations in System (1) is multiplied with suitable test functions (see [33][34][35]39] for more details) and integrated over I m . The discrete solution u 1 /I m can be represented by the polynomial ansatz where U m j are the members of the function space V and the basis functions ϕ m, j ∈ P k (I m ) are chosen as Lagrange basis functions w. r. t. the k + 1 points t m, j ∈ I m with the following assumption where δ i, j the usual Kronecker delta. We choose the points as t m, 0 = t m−1 and t m, 1 = t m the (k + 1)-quadrature points of the Gauß-Lobatto formula on each time interval. In this way, the initial condition can be written as The basic functions ϕ m, j ∈ P k (s m ) of (7) are defined using the reference transformations (see [33][34][35]39] for more details). Similarly, the test basis functionsψ i ∈ P k−1 (ŝ) are defined with the appropriate choice in order to compute the coefficients (see [33][34][35]39] for details). Finally, the cGP(k)-method reads: where U m 0 = U m−1 k for m > 1 and U 1 0 = u 0 for n = 1 are the initial values α i, j and β i and are defined as: Once the above system is solved, the initial condition for the next time interval I − m+1 is set to U m+1 0 = U m k . For k = 2, the coefficients α i, j and β i, j of the cGP(2)-method are computed as follows: The cGP(2) Method Thus, the system to be solved for U m 1 , U m 2 ∈ V from the known U m 0 = U m−1 2 becomes: where U m 0 represent the initial condition at the current time interval.

The Classical Explicit Runge-Kutta Method
This the fourth order Runge-Kutta method is very well known and was developed by Kutta [54] (see [55] for more details).

Numerical Results and Discussions
In this section, we investigate the numerical solutions of the extended HIV Model (1) by using the cGP(2) method. Initial values of the variables, different parameters, and their detailed descriptions are provided in Table 1. Various parameters present in the model are investigated with the help of graphs to understand their behavior. Figures 1-3 depict the healthy cells, infected cells, and virus particles for different values of the growth rate of healthy cells "β T ", respectively. It is concluded that the concentration of healthy cells decreases, whereas the distribution of infected cells and the virus increases. Furthermore, it could be seen clearly from the previously mentioned figures that all the distributions display decaying oscillatory behavior.                                            Figures 16-19, respectively.   Figures 13-15 indicate that when the death rate of healthy cells " " increases, the concentration of healthy/infected cells and virus decreases. From the graph, we obtained that the number of healthy cells increases initially, but after some days these cells decrease due to the continuous death of healthy cells and eventually reaches the minimum stage, which shows that the virus attacks the healthy cells. From the aforementioned discussion, all the distributions show decaying oscillatory behaviors that demonstrate the fight among the immune system of the body and the HIV virus. The amplitude of the oscillations gradually decreasing shows that the immunity of the body decreases with the passage of time during infection to fight against the infectious virus. The phase diagram of I Figures 16-19, respectively.

Comparison between the Results of Proposed Method and other Classical Methods
In this section, we describe the comparison of the solutions of the proposed scheme with those solutions already obtained in the literature. In Tables 2-4, the approximate solutions obtained from the proposed scheme are compared with those obtained from the Laplace Adomian decomposition method with Pade approximation [29], the Bessel

Comparison between the Results of Proposed Method and other Classical Methods
In this section, we describe the comparison of the solutions of the proposed scheme with those solutions already obtained in the literature. In Tables 2-4, the approximate solutions obtained from the proposed scheme are compared with those obtained from the Laplace Adomian decomposition method with Pade approximation [29], the Bessel Collocation method [31], modified variational iteration method [33], an exponential Galerkin method (EGM) [50], and classical Runge-Kutta method of the fourth order. In order to illustrate the solutions, the graphs of the RK4-method solutions and the proposed method solutions are depicted in Figures 20-22, for T(t), I(t), and V(t), respectively. Lastly, for the purpose of having an idea about the accuracy among the approximate solutions obtained from the aforesaid methods relative to the RK4-method, estimations of the absolute errors are demonstrated in Tables 5-7. A close inspection of the results is displayed in Tables 2-7 and Figures 20-22. The results suggest that our approximate solutions are closer to solutions of the RK4-method compared with the solutions of the aforementioned schemes. It is observed that the proposed scheme is reliable for finding approximate solutions to real-world problems.  Figure 19. The phase diagram of ( ), ( ), and ( ).

Comparison between the Results of Proposed Method and other Classical Methods
In this section, we describe the comparison of the solutions of the proposed scheme with those solutions already obtained in the literature. In Tables 2-4, the approximate solutions obtained from the proposed scheme are compared with those obtained from the Laplace Adomian decomposition method with Pade approximation [29], the Bessel Collocation method [31], modified variational iteration method [33], an exponential Galerkin method (EGM) [50], and classical Runge-Kutta method of the fourth order. In order to illustrate the solutions, the graphs of the RK4-method solutions and the proposed method solutions are depicted in Figures 20-22, for T(t), I(t), and V(t), respectively. Lastly, for the purpose of having an idea about the accuracy among the approximate solutions obtained from the aforesaid methods relative to the RK4-method, estimations of the absolute errors are demonstrated in Tables 5-7. A close inspection of the results is displayed in Tables 2-7 and Figures 20-22. The results suggest that our approximate solutions are closer to solutions of the RK4-method compared with the solutions of the aforementioned schemes. It is observed that the proposed scheme is reliable for finding approximate solutions to real-world problems.

Conclusions
The aim of the present study was to modify the HIV model and introduce a variable source term for the generation of healthy cells dependent on the virus and numerically approximate the density of healthy cells, infected cells, and free HIV virus after infection. We implemented a new numerical technique, namely, the continuous Galerkin-Petrov scheme for efficient and accurate solutions with the aforementioned model In addition, the impact of various physical parameters was analyzed, and all the observations were presented graphically. The main findings of the current study are summarized as follows: i.
Increasing growth rate of healthy cells, (β T ), shows a decreasing effect in the population dynamics of healthy cells, while showing an increasing effect in the population dynamics of infected cells and HIV particles. All the profiles showed a decaying oscillatory behavior. ii.
The healthy cells and infected cells show an increasing effect, while free virus distribution shows a decreasing behavior with an increase in the values of the virus death rate (α V ). iii.
It is noticed that the virus particles released by infected cells (φ) show significant variations in the population distributions of healthy cells, infected cells, and the virus. By increasing the value of "φ", the healthy cells, infected cells, and the virus increases. iv.
The graphical trends illustrate increased decay in distributions of all dependent variables with an increase in the death rate of infected cells (r I ). v.
The decrease in the density of healthy cells, infected cells, and free HIV particles is observed by increasing "α T ".
Moreover, we performed an analysis for the reproduction of number "R 0 " and concluded that when "R 0 ≤ 1" the disease-free equilibrium state "E 0 " is locally as well as globally stable, while when "R 0 > 1", then "E 0 " is locally and globally unstable. Furthermore, the well-known classical method for the initial value problem, namely, the RK-method, is utilized for the model and all the outcomes were compared with those from the proposed scheme and other schemes available in the open literature. It was observed that the suggested scheme is efficient and accurate in comparison with the other schemes. The proposed scheme could be applicable for solving complex real-world problems. In the future, the presented work is extendable in different directions. For instance, the authors are interested to determine the most proper and effective method of vaccination and treatment using the model.