BMRI BioMed Research International 2314-6141 2314-6133 Hindawi Publishing Corporation 10.1155/2014/618404 618404 Research Article Computational Analysis of the Model Describing HIV Infection of CD4+T Cells Atangana Abdon 1 http://orcid.org/0000-0001-6520-1039 Doungmo Goufo Emile Franc 2 Oukouomi Noutchie Suares Clovis 1 Institute for Groundwater Studies Faculty of Natural and Agricultural Sciences University of the Free State Bloemfontein 9300 South Africa ufs.ac.za 2 Department of Mathematical Sciences University of South Africa Florida Sciences Campus Florida 0003 South Africa unisa.ac.za 2014 1672014 2014 05 06 2014 17 06 2014 17 7 2014 2014 Copyright © 2014 Abdon Atangana and Emile Franc Doungmo Goufo. 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.

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.

1. 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 . 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 . 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 mm3. This model is characterized by a system of the nonlinear differential equations: (1) d T d t = p - α T + r T ( 1 - T + I T max ) - k V T , d I d t = k V T - β I , d V d t = N β I - γ V , subject to the initial conditions (2) T ( 0 ) = T 0 , I ( 0 ) = I 0 , V ( 0 ) = V 0 . Here, r is any positive constant, T ( t ) , I ( t ) , and V ( t ) 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 - ( T + I ) / T max ) describes the logistic growth of the healthy CD4+T cells, and proliferation of infected CD4+T cells is neglected. For k > 0 is the infection rate, the term K V T 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 p . T cells multiply through mitosis with a rate r when T cells are stimulated by antigen or mitogen. T max denotes the maximum CD4+T cell concentration in the body .

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 . Though computer science is growing very fast and numerical simulation is applied everywhere, nonnumerical issues will still play a large role . There exist some alternative analytical asymptotic approaches, such as the nonperturbative method, modified Lindstedt-Poincaré method , variational iteration method , Adomian decomposition method , homotopy perturbation method [16, 17], and book-keeping artificial parameter perturbation method .

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 one- and 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 , 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 one-dimensional Keller-Segel equations. This method displays some advantages over the existing methods.

2. 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 t ; this will further imply that (3) 0 = p - α T - + r T - ( 1 - T - + I T max ) - k V - T - , 0 = k V - T - - β I    - , 0 = N β I - - γ V - . After the first manipulations, we obtain that (4) I - = k V - T - β = γ V - N β . From the above equation, we can obtain (5) T - = γ N β . However, replacing (5) into the first equation of system (3) and rearranging it, we obtain the following: (6) I - = p + ( ( r γ - γ 2 ) / N β ) - ( r γ 2 / ( N β 2 ) T max ) ( r γ / N β T max ) + β , V - = N β γ p + ( ( r γ - γ 2 ) / N β ) - ( r γ 2 / ( N β 2 ) T max ) ( r γ / N β T max ) + β . Therefore, the equilibrium points are given as (7) ( I - , V - , T - ) = ( p + ( ( r γ - γ 2 ) / N β ) - ( r γ 2 / ( N β 2 ) T max ) ( r γ / N β T max ) + β ,    N β γ p + ( ( r γ - γ 2 ) / N β ) - ( r γ 2 / ( N β 2 ) T max ) ( r γ / N β T max ) + β , γ N β p + ( ( r γ - γ 2 ) / N β ) - ( r γ 2 / ( N β 2 ) T max ) ( r γ / N β T max ) + β ) . The above equilibrium points are valid if the following conditions are satisfied: (8) p + r γ - γ 2 N β - r γ 2 ( N β ) 2 T max 0 . The disease free equilibrium is obtained by solving the following equation: (9) - r T max T - 2 + ( r - α ) T - + p = 0 , which has the following solution: (10) T - = r - α ( r - α ) 2 + ( 4 r / T max ) 2 ( r / T max ) . Therefore, the diseases free equilibrium is given as (11) ( I - , V - , T - ) = ( 0,0 , r - α + ( r - α ) 2 + ( 4 r / T max ) 2 ( r / T max ) ) .

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 (12) J = ( r - α - k V - 2 r T T max - T T max - k T k V - β k T 0 N β - γ ) . Now, using the equilibrium free point, we obtain the following: (13) J = ( r - α - 2 r T free T max - T free T max - k T 0 - β k T free 0 N β - γ ) . In order to find the eigenvalues associate, we will solve the following equation: (14) det ( J - λ I ) = ( r - α - 2 r T free T max - λ - T free T max - k T k V - β - λ k T free 0 N β - γ - λ ) = 0 . Using the standard method, we obtain the following solution: (15) λ 1 = r - α - 2 r T T max , λ 2 = - β - γ - ( β - γ ) 2 + 4 N β T free 2 , λ 3 = - β - γ + ( β - γ ) 2 + 4 N β T free 2 . With the Eigen values in hand, we conclude that there is stability if and only if the following conditions are observed: (16) r - α - 2 r T T max < 0 , - β - γ + ( β - γ ) 2 + 4 N β T free < 0 .

3. 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]): (17) m U ( t ) t m = L ( U ( t ) ) + N ( U ( t ) ) + f ( t ) , m = 1,2 , 3 , . Subjected to the initial condition, (18) i U ( 0 ) t i = y i , i = 0,1 , 2 , , m - 1 . m is the order of the derivative.

f is a known function, N is the general nonlinear differential operator, L represents a linear differential operator, and m is the order of the derivative. The method’s first step here is to apply the inverse operator m / t m on both sides of (17) to obtain (19) U ( t ) = k = 0 m - 1 t k k ! y k + 0 t 0 t 1 0 t m - 1 L ( U ( τ ) ) + N ( U ( τ ) ) + f ( τ ) d τ d t . The multi-integral in (17) can be transformed to (20) 0 t 0 t 1 0 t m - 1 L ( U ( τ ) ) + N ( U ( τ ) ) + f ( τ ) d τ d t = 1 ( m - 1 ) ! 0 t ( t - τ ) m - 1 L ( U ( τ ) ) + N ( U ( τ ) ) + f ( τ ) d τ so that (19) can be reformulated as (21) U ( t ) = k = 0 m - 1 t k k ! y k + 1 ( m - 1 ) ! × 0 t ( t - τ ) m - 1 L ( U ( τ ) ) + N ( U ( τ ) ) + f ( τ ) d τ . Using the homotopy scheme, the solution of the above integral equation is given in series form as follows: (22) U ( t , p ) = n = 0 p n U n ( t ) , U ( t ) = lim p 1 U ( t , p ) , and the nonlinear term can be decomposed as (23) N U ( t ) = n = 1 p n H n ( U ) , where p ( 0 , 1 ] is an embedding parameter. H n ( U ) is the He’s polynomials that can be generated by (24) H n ( U 0 , , U n ) = 1 n ! n p n [ N ( j = 0 n p j U j ( t ) ) ] , 1111111111111111111111111111 n = 0,1 , 2 , . The homotopy decomposition method is obtained by the graceful coupling of decomposition method with He’s polynomials and is given by (25) n = 0 p n U n ( t ) = T ( t ) + p 1 ( m - 1 ) ! × 0 t ( t - τ ) m - 1 [ f ( τ ) + L ( n = 0 p n U n ( τ ) ) + n = 0 p n H n ( U ) ] d τ , with (26) T ( t ) = k = 0 m - 1 t k k ! y k . Comparing the terms of the same powers of p , we obtain solutions of various orders. The initial guess of the approximation is T ( t ) ; this is actually the Taylor series of the exact solution of order m . Note that this initial guess insures the uniqueness of the series decompositions.

4. 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]: (27) T 0 = 0.1 , I 0 = 0 , V 0 = 0.1 , p = 0.1 , α = 0.02 , β = 0.3 , γ = 0.3 , k = 0.0027 ,    T max = 1500 , N = 10 . Following carefully the steps of the HDM, we obtain the following equations: (28) n = 0 p n T n ( t ) = T ( 0 ) + p 0 t ( p - α n = 0 p n T n ( τ ) + r n = 0 p n T n ( τ ) × ( 1 - n = 0 p n T n ( τ ) + n = 0 p n I n ( τ ) T max ) - k n = 0 p n T n ( τ ) n = 0 p n V n ( τ ) ) d τ , n = 0 p n I n ( t ) =    I ( 0 ) + p 0 t ( k n = 0 p n T n ( τ ) n = 0 p n V n ( τ ) - β n = 0 p n I n ( τ ) ) d τ , n = 0 p n V n ( t ) = V ( 0 ) + p 0 t ( β N n = 0 p n T n ( τ ) - γ n = 0 p n V n ( τ ) ) d τ . Comparing the terms of the same power of p , 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 p , which something is very heavy to compute in the case of high order ODE. Consider (29) p 0 = T 0 ( t ) = T 0 , T 0 ( 0 ) = T 0 , p 0 = I 0 ( t ) = I 0 ,    I 0 ( 0 ) = I 0 , p 0 = V 0 ( t ) = V 0 ,    V 0 ( 0 ) = V 0 , p 1 : T 1 ( t ) = 0 t ( p - α T 0 + r T 0 ( 1 - T 0 + I 0 T max ) - k V 0 T 0 ) d τ ,    T 1 ( 0 ) = 0 p 1 : I 1 ( t ) = 0 t ( k V 0 T 0 - β I 0 ) d τ , I 1 ( 0 ) = 0 p 1 : V 1 ( t ) = 0 t ( N I 0 - γ V 0 ) d τ , V 1 ( 0 ) = 0 p n : T n ( t ) = 0 t ( ( j = 0 n - 1 T j T n - j - 1 + j = 0 n - 1 T j I n - j - 1 ) ( r - α ) T n - 1    - r T max ( j = 0 n - 1 T j T n - j - 1 + j = 0 n - 1 T j I n - j - 1 ) - k j = 0 n - 1 T j V n - j - 1 ) d τ , T n ( t ) = 0 , n 2 p n : I n ( t ) = 0 t ( k j = 0 n - 1 V j T n - j - 1 - β I n - 1 ) d τ , I n ( t ) = 0 ,    n 2 p n : V n ( t ) = 0 t ( N β I n - 1 - γ V n - 1 ) d τ , V n ( t ) = 0 ,    n 2 . We will give the general algorithm in order to accommodate scholars using computers program; the algorithm will then be used to derive special solution of the system of equation for a given set of theoretical parameters.

Algorithm 1.

input: J 1 ( t ) , J 2 ( t ) , and J 3 ( t ) as initial guest;

j -number terms in the rough calculation;

output: T approx ( t ) , I approx ( t ) , and V approx ( t ) the approximate solutions.

Step 1. Put T 0 ( t ) = T 0 , I 0 ( 0 ) = I 0 , V 0 ( t ) = V 0 , T approx ( t ) = T 0 ( t ) ,    I approx ( t ) = I 0 ( t )    and    V approx ( t ) = V 0 ( t ) .

Step 2. From j = 0 to n - 1 , do Steps 3 and Step 4.

Step 3. Compute (30) T n ( t ) = 0 t ( ( r - α ) T n - 1 - r T max ( j = 0 n - 1 T j T n - j - 1 + j = 0 n - 1 T j I n - j - 1 ) - k j = 0 n - 1 T j V n - j - 1 ) d τ I n ( t ) = 0 t ( k j = 0 n - 1 V j T n - j - 1 - β I n - 1 ) d τ , V n ( t ) = 0 t ( N β I n - 1 - γ V n - 1 ) d τ .

Step 4. Compute T approx ( t ) = T approx ( t ) + T n + 1 ( t ) , I approx ( t ) = I approx ( t ) + I n + 1 ( t ) , and V approx ( t )    =    V approx ( t ) + V n + 1 ( t ) .

Stop.

We will now make use of the above algorithm to derive the special solution. We therefore obtain the following series solutions: (31) T 0 ( t ) = T 0 = 0.1 , I 0 ( t ) = I 0 = 0 , V 0 ( t ) = V 0 = 0.1 , T 1 ( t ) = 0.397953 t , I 1 ( t ) = 0.000027 t , V 1 ( t ) = - 0.24 t , T 2 ( t ) = 0.592849 t 2 , I 2 ( t ) = 0.0000172737 t 2 , V 2 ( t ) = 0.288041 t 2 , T 3 ( t ) = 0.588719 t 3 , I 3 ( t ) = - 8.40515 × 10 - 6 t 3 , V 3 ( t ) = - 0.230415 t 3 , T 4 ( t ) = 0.43829519 t 4 , I 4 ( t ) = 6.14727882 × 10 - 6 t 4 , V 4 ( t ) = 0.1382427719 t 4 , T 5 ( t ) = 0.2608633362 t 5 , I 5 ( t ) = 2.835861862 × 10 - 6 t 5 , V 5 ( t ) = - 6.635284213 × 10 - 2 t 5 , T 6 ( t ) = 0.1291947662 t 6 , I 6 ( t ) = 1.153299855 × 10 - 5 t 6 , V 4 ( t ) = 2.653971892 t 6 , T 7 ( t ) = 0.01363544556 t 7 , I 7 ( t ) = - 3.925609032 × 10 - 7 t 7 , V 7 ( t ) = 2.481235575 × 10 - 6 t 7 , T 8 ( t ) = - 2.787485189 × 10 - 5 t 8 , I 8 ( t ) = 4.5660570151 × 10 - 11 t 8 , V 8 ( t ) = 1.28523993 × 10 - 10 t 5 . 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) T ( t ) = T 0 ( t ) + T 1 ( t ) + T 2 ( t ) + T 3 ( t ) + T 4 ( t ) + T 5 ( t ) + T 6 ( t ) + T 7 ( t ) + T 8 ( t ) + I ( t ) = I 0 ( t ) + I 1 ( t ) + I 2 ( t ) + I 3 ( t ) + I 4 ( t ) + I 5 ( t ) + I 6 ( t ) + I 7 ( t ) + I 8 ( t ) + V ( t ) = V 0 ( t ) + V 1 ( t ) + V 2 ( t ) + V 3 ( t ) + V 4 ( t ) + V 5 ( t ) + V 6 ( t ) + V 7 ( t ) + V 8 ( t ) + .

4.1. 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 of T ( t ) for N = 8 .

0 0.1 0.1 0.1 0.1 0.1 0.1
0.2 0.2088072731 0.2088072731 0.2088080833 0.2088080868 0.2088073214 0.2038616561
0.4 0.4061052625 0.4061052625 0.4062405393 0.4062407949 0.4061346587 0.3803309335
0.6 0.7611467713 0.7611467713 0.7644238890 0.7644287245 0.7624530350 0.6954623767
0.8 1.3773198590 1.3773198590 1.4140468310 1.4140941730 1.3978805880 1.2759624442
1 2.3291697610 2.3291697610 2.5915948020 0.2088080868 2.5067466690 2.3832277428

Numerical comparison of I ( t ) for N = 8 .

0 0 0 0 0.1 × 10 - 13 0 0
0.2 6.03270728 · 10−6 6.03270728 · 10−6 6.032702150 · 10−6 6.032701651 · 10−6 6.032634366 · 10−6 6.247872100 · 10−6
0.4 1.31591617 · 10−5 1.31591617 · 10−5 1.315834073 · 10−5 1.315830167 · 10−5 1.314878543 · 10−5 1.293552225 · 10−5
0.6 2.12683688 · 10−5 2.12683688 · 10−5 2.122378506 · 10−5 2.1223310013 · 10−5 2.101417193 · 10−5 2.035267183 · 10−5
0.8 3.00691867 · 10−5 3.00691867 · 10−5 3.017741955 · 10−5 3.0174509323 · 10−5 2.795130456 · 10−5 2.837302120 · 10−5
1 3.98736542 · 10−5 3.98736542 · 10−5 4.003781468 · 10−5 4.0025404050 · 10−5 2.431562317 · 10−5 3.690842367 · 10−5

Numerical comparison of V ( t ) for N = 8 .

0 0.1 0.1 0.1 0.1 0.1 0.1
0.2 0.06187996025 0.06187996025 0.06187984331 0.06187990876 0.06187995314 0.06187991856
0.4 0.03831324883 0.03831324883 0.03829488788 0.03829595768 0.03830820126 0.03829493490
0.6 0.02439174349 0.02439174349 0.02370455014 0.02371029480 0.02392029257 0.02370431860
0.8 0.009967218934 0.009967218934 0.01468036377 0.01470041902 0.01621704553 0.01467956982
1 0.003305076447 0.003305076447 0.009100845043 0.009157238735 0.01608418711 0.02370431861

Numerical comparison shows that the approximated solutions obtained via LADM-Padé  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 , 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.

5. 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 the technique with those obtained with other methods. 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 . 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.

Conflict of Interests

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

Ho D. D. Neumann A. U. Perelson A. S. Chen W. Leonard J. M. Markowitz M. Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection Nature 1995 373 6510 123 126 10.1038/373123a0 2-s2.0-0028874048 Mohri H. Bonhoeffer S. Monard S. Perelson A. S. Ho D. D. Rapid turnover of T lymphocytes in SIV-infected rhesus macaques Science 1998 279 5354 1223 1227 10.1126/science.279.5354.1223 2-s2.0-0032549086 Wei X. Ghosh S. K. Taylor M. E. Johnson V. A. Emini E. A. Deutsch P. Lifson J. D. Bonhoeffer S. Nowak M. A. Hahn B. H. Saag M. S. Shaw G. M. Viral dynamics in human immunodeficiency virus type 1 infection Nature 1995 373 6510 117 122 10.1038/373117a0 2-s2.0-0028811974 Perelson A. S. Neumann A. U. Markowitz M. Leonard J. M. Ho D. D. HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time Science 1996 271 5255 1582 1586 10.1126/science.271.5255.1582 2-s2.0-0029967721 Perelson A. S. Kirschner D. E. De Boer R. Dynamics of HIV infection of CD4+ T cells Mathematical Biosciences 1993 114 1 81 125 10.1016/0025-5564(93)90043-A 2-s2.0-0027457101 Perelson A. S. Nelson P. W. Mathematical analysis of HIV-1 dynamics in vivo SIAM Review 1999 41 1 3 44 10.1137/S0036144598335107 MR1669741 2-s2.0-0001918317 Wang L. Li M. Y. Mathematical analysis of the global dynamics of a model for HIV infection of CD4+ T cells Mathematical Biosciences 2006 200 1 44 57 10.1016/j.mbs.2005.12.026 MR2211927 2-s2.0-33644556174 Asquith B. Bangham C. R. M. The dynamics of T-cell fratricide: application of a robust approach to mathematical modelling in immunology Journal of Theoretical Biology 2003 222 1 53 69 10.1016/S0022-5193(03)00013-4 MR2069578 2-s2.0-0037466387 Nowak M. A. May R. M. Mathematical biology of HIV infections: antigenic variation and diversity threshold Mathematical Biosciences 1991 106 1 1 21 10.1016/0025-5564(91)90037-J 2-s2.0-0025871325 Andrianov I. V. Manevitch L. I. Asymptotology: Ideas, Methods, and Applications 2003 Kluwer Academic Publishers Andrianov I. Awrejcewicz J. Construction of periodic solutions to partial differential equations with non-linear boundary conditions International Journal of Nonlinear Sciences and Numerical Simulation 2000 1 4 327 332 10.1515/IJNSNS.2000.1.4.327 MR1793840 2-s2.0-0034553955 Bender C. M. Milton K. A. Pinsky S. S. Simmons J. A new perturbative approach to nonlinear problems Journal of Mathematical Physics 1989 30 7 1447 1455 10.1063/1.528326 MR1002247 ZBL0684.34008 2-s2.0-0000901190 Delamotte B. Nonperturbative method for solving differential equations and finding limit cycles Physical Review Letters 1993 70 22 3361 3364 10.1103/PhysRevLett.70.3361 MR1218585 2-s2.0-0001241171 Wu G. Baleanu D. Variational iteration method for fractional calculus—a universal approach by Laplace transform Advances in Difference Equations 2013 2013 article 18 10.1186/1687-1847-2013-18 MR3020186 2-s2.0-84874035986 Adomian G. A review of the decomposition method in applied mathematics Journal of Mathematical Analysis and Applications 1988 135 2 501 544 10.1016/0022-247X(88)90170-9 ZBL0671.34053 2-s2.0-0041185368 Atangana A. Secer A. The time-fractional coupled-Korteweg-de-Vries equations Abstract and Applied Analysis 2013 2013 8 947986 MR3035318 10.1155/2013/947986 Wu G. C. Baleanu D. Variational iteration method for the Burgers' flow with fractional derivatives—new Lagrange multipliers Applied Mathematical Modelling 2013 37 9 6183 6190 10.1016/j.apm.2012.12.018 MR3039000 2-s2.0-84875378030 Wu G. C. New trends in the variational iteration method Communications in Fractional Calculus 2011 2 59 75 Atangana A. Botha J. F. Analytical solution of the groundwater flow equation obtained via homotopy decomposition method Journal of Earth Science & Climatic Change 2012 3 115 10.4172/2157-7617.1000115 Ongun M. Y. The Laplace adomian decomposition method for solving a model for HIV infection of CD4+ T cells Mathematical and Computer Modelling 2011 53 5-6 597 603 10.1016/j.mcm.2010.09.009 MR2769431 2-s2.0-78751604821 Atangana A. Alabaraoye E. Solving a system of fractional partial differential equations arising in the model of HIV infection of CD4+ cells and attractor one-dimensional Keller-Segel equations Advance in Difference Equation 2013 2013, article 94 10.1186/1687-1847-2013-94 Ş. Yüzbaşı A numerical approach to solve the model for HIV infection of CD4+T cells Applied Mathematical Modelling 2012 36 12 5876 5890 10.1016/j.apm.2011.12.021 MR2968433 2-s2.0-84864986117 Atangana A. New class of boundary value problems Applied Mathematics & Information Sciences 2012 1 2 67 76