Computational Analysis of the Model Describing HIV Infection of CD4+T Cells

An analysis of the model underpinning the description of the spread of HIV infection of CD4+T cells is examined in detail in this work. Investigations of the disease free and endemic equilibrium are done using the method of Jacobian matrix. An iteration technique, namely, the homotopy decomposition method (HDM), is implemented to give an approximate solution of nonlinear ordinary differential equation systems. The technique is described and illustrated with numerical examples. The approximated solution obtained via HDM is compared with those obtained via other methods to prove the trustworthiness of HDM. Moreover, the lessening and simplicity in calculations furnish HDM with a broader applicability.


Introduction
The cause of the progressive depletion of CD4 + T cells in HIV-infected people is one of the most fundamental and controversial issues in AIDS research. HIV infects and kills CD4 + T cells. The infection results in high T-cell activation and turnover. An immediately intuitive assumption is that HIV-mediated destruction of CD4 + cells directly reduces the number of these cells and that the high turnover rates of T cells and the slow progression to AIDS reflect a long but eventually lost struggle of the immune system to replace killed cells in its effort to maintain T-cell homeostasis [1][2][3][4]. However, HIV mainly infects activated CD4 + cells, and activated cells normally follow different dynamics than cells that belong to resting populations whose numbers are controlled by homeostatic mechanisms.
In this study, we consider that the HIV infection model of CD4 + T cells is examined in [5]. This model is given by the components of the basic three-component model which are the concentration of susceptible CD4 + T cells, CD4 + T cells infected by the HIV viruses, and free HIV virus particles in the blood. CD4 + T cells are also called leukocytes or T helper cells. These with order cells in human immunity systems fight against diseases. HIV use cells in order to propagate. In a healthy person, the number of CD4 + T cells is 800/1200 mm 3 . This model is characterized by a system of the nonlinear differential equations: subject to the initial conditions Here, is any positive constant, ( ), ( ), and ( ) show the concentration of susceptible CD4 + T cells, CD4 + T cells infected by the HIV viruses, and free HIV virus particles in the blood, respectively, , , and stand for natural turnover rates of uninfected Tells, infected T cells, and virus particles, respectively, (1−( + )/ max ) describes the logistic growth of the healthy CD4 + T cells, and proliferation of infected CD4 + T cells is neglected. For > 0 is the infection rate, the term describes the incidence of HIV infection of healthy CD4 + T cells. Each infected CD4 + T cell is assumed to produce l virus particles during its lifetime, including any of its daughter cells. The body is believed to produce CD4 + T cells from precursors in the bone marrow and thymus at a constant rate . T cells multiply through mitosis with a rate when T cells are stimulated by antigen or mitogen. max denotes the maximum CD4 + T cell concentration in the body [6][7][8][9]. As Andrianov and Manevitch wrote in the foreword of a popular science book Asymtotology: Ideas, Methods, and Applications, asymptotic methods belong to the, perhaps, most romantic area of modern mathematics [10]. Though computer science is growing very fast and numerical simulation is applied everywhere, nonnumerical issues will still play a large role [11][12][13][14]. There exist some alternative analytical asymptotic approaches, such as the nonperturbative method, modified Lindstedt-Poincaré method [13], variational iteration method [14], Adomian decomposition method [15], homotopy perturbation method [16,17], and book-keeping artificial parameter perturbation method [18].
The purpose of this paper is to derive analytical solutions of model for HIV infection of CD4 + T cells (1) via the relatively new analytical solution homotopy decomposition method (HDM). The HDM was recently used to solve oneand two-dimensional fractional heat-like problem, Burgers equation, the Klein-Gordon equation, a coupled Burgers equation, 2D and 3D Poisson equations and biharmonic equations, the groundwater flow equation [19], the Hamilton-Jacobi-Bellman equation, high even-order differential equation, parabolic equations and coupled Van der Pol's nonlinear partial differential equation, and coupled attractor onedimensional Keller-Segel equations. This method displays some advantages over the existing methods.

Stability Analysis
The aim of this section is to present a stability analysis of the system equation (1) that will be achieved via the use the eigenvalues obtained via the Jacobian matrix. We will start by providing the equilibrium point and also the desease free equilibrium. To obtain the equilibrium points, we assume that the system does not depend on the parameter ; this will further imply that After the first manipulations, we obtain that From the above equation, we can obtain However, replacing (5) into the first equation of system (3) and rearranging it, we obtain the following: Therefore, the equilibrium points are given as ) .
The above equilibrium points are valid if the following conditions are satisfied: The disease free equilibrium is obtained by solving the following equation: which has the following solution: Therefore, the diseases free equilibrium is given as With the above information on hand, we will now find the eigenvalues associated with this problem, which will allow us to give stability of the system. The Jacobian matrix associated with this problem is given by Now, using the equilibrium free point, we obtain the following: In order to find the eigenvalues associate, we will solve the following equation: Using the standard method, we obtain the following solution: With the Eigen values in hand, we conclude that there is stability if and only if the following conditions are observed:

Some Useful Information regarding the Methodology of (HDM)
To illustrate the basic idea of this method, we consider a general nonlinear nonhomogeneous differential equation with initial conditions of the following form (see also in [16,19]): Subjected to the initial condition, is the order of the derivative. is a known function, is the general nonlinear differential operator, represents a linear differential operator, and is the order of the derivative. The method's first step here is to apply the inverse operator / on both sides of (17) to obtain The multi-integral in (17) can be transformed to so that (19) can be reformulated as Using the homotopy scheme, the solution of the above integral equation is given in series form as follows: and the nonlinear term can be decomposed as where ∈ (0, 1] is an embedding parameter. H ( ) is the He's polynomials that can be generated by 4

BioMed Research International
The homotopy decomposition method is obtained by the graceful coupling of decomposition method with He's polynomials and is given by with Comparing the terms of the same powers of , we obtain solutions of various orders. The initial guess of the approximation is ( ); this is actually the Taylor series of the exact solution of order . Note that this initial guess insures the uniqueness of the series decompositions.

Applications
In this section, we applied the HDM to solve the system of the nonlinear differential equations (1). To be consistent in comparison with the existing methods, we chose the following initial condition and parameters as in [20,21]: Following carefully the steps of the HDM, we obtain the following equations: Comparing the terms of the same power of , we obtain the following integral equations that are much easier to solve. Note that, with the homotopy perturbation method (HPM), one will obtain a set of ordinary differential equations, after comparing the terms of the same power of , which something is very heavy to compute in the case of high order ODE. Consider We will give the general algorithm in order to accommodate scholars using computers program; the algorithm will then be BioMed Research International 5 used to derive special solution of the system of equation for a given set of theoretical parameters.
We will now make use of the above algorithm to derive the special solution. We therefore obtain the following series solutions: Using the package Mathematica, in the same manner, one can obtain the rest of the components. But, in this case, 9 terms were computed and the asymptotic solution is given by (32)

Numerical Applications.
To test the effectiveness and the accuracy of the HDM for solving this type of problem, we compare the approximated solutions obtained via other methods and HDM, and the results are shown in Tables 1, 2, and 3. Numerical comparison shows that the approximated solutions obtained via LADM-Padé [20] are in good agreement with the results obtained via HDM. More precisely, the approximated solutions are exactly the same. However, in the technique used in [20], one needs first to apply the Laplace transform on the system, following by the ADM, and finally take the inverse Laplace transform to obtain the approximated solutions, which is much time consuming and sometimes can lead to a very difficult situation, for example, if the inverse Laplace transform cannot be obtained.

Conclusions
In this paper, homotopy decomposition method has been developed for finding approximate solutions of HIV infection model of CD4 + T which is a class of nonlinear ordinary differential equation systems. We have demonstrated the accuracy and efficiency of the present technique with an example and comparison of the approximate solution obtained via   Comparing the methodology HDM to homotopy perturbation method (HPM), Adomian decomposition method (ADM), variational iteration method (VIM), and homotopy analysis method (HAM) has advantages. Disparate the ADM, the HDM is free from the need of Adomian polynomials.
In this method, we do not need the Lagrange multiplier, correction functional, stationary conditions, or calculation of heavy integrals; the solution obtained is noise free, which eliminate the complications that exist in the VIM [23]. In contract to HPM, we do not need to continuously deform a difficult problem to another that is easier to solve. In contrast to LADM, the technique does not need any change of space, which sometime can be very difficult situation to handle, for example, in case the inverse Laplace transform cannot be found. We can easily conclude that the homotopy decomposition method is an efficient tool to solve approximate solution of nonlinear system partial differential equations.