# Adomian Decomposition Method Combined with Padé Approximation and Laplace Transform for Solving a Model of HIV Infection of CD4+T Cells

Authors: Fang Chen; Qing-Quan Liu
Journal: Discrete Dynamics in Nature and Society
 (2015)

---

## Abstract

The classical Adomian decomposition method (ADM) is implemented to solve a model of
HIV infection of CD4+T cells. The results indicate that the approximate solution by using
the ADM is the same as that by using the Laplace ADM, but it can be obtained in a more
efficient way. We also use Padé approximation and Laplace transform as a posttreatment
technique to obtain the result of the ADM. The advantage of the posttreatment is illustrated
by numerical experiments.


---

## Body

## 1. Introduce

Because of human immunodeficiency virus (HIV), more and more people are infected with more and more casualties especially in Africa. CD4
+T cells are the most abundant white blood cells of the immune system in body. Though HIV infects also other cells, it seriously damages the CD4
+T cells in blood and this decreases the immune ability. Mathematical modeling is an important tool established for understanding HIV infection concerning CD4
+T cells. In 1989, Perelson established a simple and important model in the field of HIV infection [1]. In 1993, Perelson et al. proposed a fitted model based on an old model which appeared in [2]. This HIV infection model of CD4
+T cells is given by the system of nonlinear differential equations: (1)
d
T
d
t
=
p
-
α
T
+
r
T
1
-
T
+
I
T
m
a
x
-
k
V
T
,
d
I
d
t
=
k
V
T
-
β
I
,
d
V
d
t
=
N
β
I
-
γ
V
,with the initial conditions (2)
T
0
=
r
1
,
I
0
=
r
2
,
V
0
=
r
3
.This model neglects the proliferation of infected CD4
+T cells. Here T
(
t
), I
(
t
), and V
(
t
) describe the concentration of susceptible CD4
+T cells, CD4
+T cells infected by the HIV, and free HIV particles in the blood, respectively; α, β, and γ represent natural turnover rates of uninfected T cells, infected T cells, and virus particles, respectively; k
(
>
0
) is the infection rate; p denotes a rate at which the body produces CD4
+T cells from precursors in the bone narrow and thymus; r represents a rate at which T cells multiply through mitosis when the T cells are stimulated by antigen or mitogen; N means the virus particles that each infected CD4
+T cell produces during its life, including all its daughter cells; T
m
a
x indicates the maximum CD4
+T cells concentration in the body. Also, (
1
-
(
T
+
I
)
/
T
m
a
x
) describes the logistic growth of healthy CD4
+T cells. In this paper, all parameters are the same as those in [3, 4]; that is, p
=
0.1, α
=
0.02, β
=
0.3, r
=
3, γ
=
2.4, k
=
0.0027, T
m
a
x
=
1500, N
=
10, r
1
=
0.1, r
2
=
0, and r
3
=
0.1.In the recent twenty years, all kinds of models of HIV infection of CD4
+T cells have been established and studied. These models are usually described by the nonlinear differential equations. When the analytic solution can not be obtained in close form, the numerical and/or approximate solutions need to be computed. Adomian decomposition method (ADM) [5] and Laplace ADM (LADM) [6] are some typical approaches, which give the approximate solution of linear or nonlinear problems [7]. Padé approximation and Laplace transform have been regarded as posttreatment tools that are used to obtain some results of nonlinear differential equations [8–14]. Gauss elimination, LU factorization, and chasing methods are classical direct methods [15], and Jacobi, Gauss-Seidel, SOR, (preconditioned) Krylov subspace, and (preconditioned) HSS methods are iteration methods for solving systems of linear equations [15–19]. Besides, Newton, Euler, Runge-Kutta, Newton-HSS, and Newton-Krylov subspace methods are numerical methods for solving systems of nonlinear equations [20–23]. In this paper, we deeply study ADM and the posttreatment technique for solving the model of HIV infection (1).The type of HIV infection model (1) has been deeply discussed in [3, 4]. Arafa et al. proposed generalized Euler method (GEM) to solve the fractional order model in [24]. Ongun used the LADM to solve the system of nonlinear differential equations (1) in [4]. Doğan proposed the multistep LADM to solve the system of nonlinear differential equations (1) in [3]. We will use ADM to solve HIV infection model (1), and we will utilize Padé approximation and Laplace transform for obtaining some new results.The remainder of this paper is organized as follows. in Section2, we apply ADM for solving HIV infection model (1). In Section 3, we use Padé approximation and Laplace transform as posttreatment technique to improve the results of ADM. In Section 4, we give numerical results for the system of nonlinear differential equations (1) and compare these results with those of the classical fourth-order Runge-Kutta (RK4) method.


## 2. Application of ADM

In this section, the classical ADM is used to solve the system of nonlinear differential equations (1), whose results are the same as that using the LADM [4]. For this class of problems, ADM is a simple and convenient method.Now, let us use ADM to solve the system of nonlinear differential equations (1). Rewriting HIV infection model (1) in an operator form, (3)
L
t
T
=
p
-
α
-
r
T
-
r
T
m
a
x
T
2
-
r
T
m
a
x
T
I
-
k
V
T
,
L
t
I
=
k
V
T
-
β
I
,
L
t
V
=
N
β
I
-
γ
V
,where the differential operator L
t denotes (4)
L
t
=
∂
∂
t
.Letting the inverse operator L
t
-
1 be an integral operator of the form (5)
L
t
-
1
=
∫
0
t
d
tand applying L
t
-
1 on both sides of the equations in (3), we obtain (6)
T
=
T
0
+
p
t
-
α
-
r
L
t
-
1
T
-
r
T
m
a
x
L
t
-
1
T
2
-
r
T
m
a
x
L
t
-
1
T
I
-
k
L
t
-
1
V
T
,
I
=
I
0
+
k
L
t
-
1
V
T
-
β
L
t
-
1
I
,
V
=
V
0
+
N
β
L
t
-
1
I
-
γ
L
t
-
1
V
.Using the ADM, we express the unknown items as the infinite series (7)
T
=
Σ
n
=
0
∞
T
n
,
I
=
Σ
n
=
0
∞
I
n
,
V
=
Σ
n
=
0
∞
V
nand rewrite the nonlinear items as the Adomian polynomials (8)
T
2
=
Σ
n
=
0
∞
A
n
,
T
I
=
Σ
n
=
0
∞
B
n
,
V
T
=
Σ
n
=
0
∞
C
n
.Based on the above expressions and the initial conditions, we can obtain the following recursive relationship: (9)
T
0
=
r
1
+
p
t
,
T
n
+
1
=
-
α
-
r
L
t
-
1
T
n
-
r
T
m
a
x
L
t
-
1
A
n
-
r
T
m
a
x
L
t
-
1
B
n
-
k
L
t
-
1
C
n
,
I
0
=
r
2
,
I
n
+
1
=
k
L
t
-
1
B
n
-
β
L
t
-
1
I
n
,
V
0
=
r
3
,
V
n
+
1
=
N
β
L
t
-
1
I
n
-
γ
L
t
-
1
V
n
,
n
=
0,1
,
2
,
…
.By the aid of the software “Mathematica,” we can write the first seven terms from the zeroth to the sixth components ofT
(
t
), I
(
t
), and V
(
t
) and then give their sums as follows: (10)
T
q
6
≔
∑
i
=
0
6
T
i
t
=
0.1
+
0.397953
t
+
0.592849
t
2
+
0.588719
t
3
+
0.438295
t
4
+
0.260863
t
5
+
0.129195
t
6
+
0.0136354
t
7
-
0.278755
×
1
0
-
4
t
8
+
1.00825
×
1
0
-
8
t
9
-
1.34165
×
1
0
-
12
t
10
+
8.06031
×
1
0
-
17
t
11
-
2.22614
×
1
0
-
21
t
12
+
2.29896
×
1
0
-
26
t
13
,
I
q
6
≔
∑
i
=
0
6
I
i
t
=
2.7
×
1
0
-
5
t
+
1.72737
×
1
0
-
5
t
2
-
8.40515
×
1
0
-
6
t
3
+
6.14728
×
1
0
-
6
t
4
-
2.83586
×
1
0
-
6
t
5
+
1.1533
×
1
0
-
6
t
6
-
3.92561
×
1
0
-
7
t
7
+
4.57012
×
1
0
-
11
t
8
-
2.67794
×
1
0
-
14
t
9
+
7.55134
×
1
0
-
19
t
10
+
1.54535
×
1
0
-
22
t
11
-
6.38153
×
1
0
-
27
t
12
,
V
q
6
≔
∑
i
=
0
6
V
i
t
=
0.1
-
0.24
t
+
0.288041
t
2
-
0.230415
t
3
+
0.138243
t
4
-
0.0663528
t
5
+
0.0265397
t
6
+
2.48124
×
1
0
-
6
t
7
+
1.28486
×
1
0
-
10
t
8
-
3.58846
×
1
0
-
15
t
9
-
1.44871
×
1
0
-
19
t
10
+
2.57662
×
1
0
-
23
t
11
.We now analyze the expressionsT
q
6, I
q
6, and V
q
6. In fact, the first seven terms in the expressions of T
q
6, I
q
6, and V
q
6 are exact because the other terms will be changed by the remaining components. That is to say, the coefficients from t
0 to t
6 are invariable. So the truncated Taylor series of the approximate solutions T
q
6, I
q
6, and V
q
6 should be of the forms (11)
T
t
6
t
=
0.1
+
0.397953
t
+
0.592849
t
2
+
0.588719
t
3
+
0.438295
t
4
+
0.260863
t
5
+
0.129195
t
6
,
I
t
6
t
=
1
0
-
5
×
2.7
t
+
1.72737
t
2
-
0.840515
t
3
+
0.614728
t
4
-
0.283586
t
5
+
0.11533
t
6
,
V
t
6
t
=
0.1
-
0.24
t
+
0.288041
t
2
-
0.230415
t
3
+
0.138243
t
4
-
0.0663528
t
5
+
0.0265397
t
6
.Note that these results are the same as those of LADM in [4]. We believe that the results by using ADM are the same as those by using LADM, because LADM is essentially the same as that of ADM. LADM firstly uses the Laplace transform to the original equations. Then, based on the idea of ADM, a recursive relationship of the LADM is obtained [4]. And lastly, an approximate solution of the unknown function is derived by using the inverse Laplace transform. From the implementation procedure, we find that LADM is more complicated than ADM, because LADM uses not only Laplace transform but also inverse Laplace transform. However, ADM only uses inverse operator. So we can directly use ADM to solve the system of nonlinear differential equations (1).


## 3. Posttreatment of ADM

In this section, we first introduce Padé approximation and then combine Laplace transform with Padé approximation to conduct the results of ADM. By using this posttreatment technique, we can obtain better approximate solutions.Assume that a functionf
(
x
) has a Taylor expansion f
(
x
)
=
∑
i
=
0
∞
c
i
x
i at the zero point. We want to compute a Padé approximation to f
(
x
). To this end, let the [
L
/
M
] Padé approximation have the form (12)
L
M
=
P
L
x
Q
M
x
=
a
0
+
a
1
x
+
⋯
+
a
L
x
L
1
+
b
1
x
+
⋯
+
b
M
x
M
.By matching the coefficients of Taylor expansion of f
(
x
) with those of [
L
/
M
] Padé approximation, we have (13)
a
0
+
a
1
x
+
⋯
+
a
L
x
L
=
1
+
b
1
x
+
⋯
+
b
M
x
M
c
0
+
c
1
x
+
⋯
.The coefficients a
0
,
a
1
,
…
,
a
L and b
1
,
b
2
,
…
,
b
M can be computed via the above equality. This means that the power series of f
(
x
) equals the [
L
/
M
] Padé approximation from x
0 to x
L
+
M, or, in other words, it holds that (14)
f
x
-
P
L
x
Q
M
x
=
O
x
L
+
M
+
1
.It is known that [
L
/
M
] Padé approximation of a function is uniquely determined [25]. We want to apply Padé approximation to the approximate solutions of ADM. Because the coefficients of the first seven terms are invariable, to obtain better approximate solutions we should adopt [
3
/
3
] Padé approximation. In [4] the author was aware of the invariance of the coefficients of the first seven terms but used the [
4
/
4
] Padé approximation. According to the characteristic of the Padé approximation, the [
3
/
3
] Padé approximation should be better than the [
4
/
4
] Padé approximation [25].By “Mathematica” software, the[
3
/
3
] Padé approximations of T
q
6
(
t
), I
q
6
(
t
), and V
q
6
(
t
) are, respectively, given by (15)
T
P
3
=
0.1
+
0.249038
t
+
0.0889628
t
2
+
0.0369466
t
3
1
-
1.48915
t
+
0.887248
t
2
-
0.220151
t
3
,
I
P
3
=
1
0
-
5
×
2.7
t
+
3.57914
t
2
+
0.50175
t
3
1
+
0.685843
t
+
0.0583565
t
2
-
0.051507
t
3
,
V
P
3
=
0.1
-
0.119945
t
+
0.0575645
t
2
-
0.011445
t
3
1
+
1.20055
t
+
0.576566
t
2
+
0.115382
t
3
.Usually, Laplace transform can make the coefficients of a power function become large, the constants disappear, and the index of unknown number increase by one. For example, applying Laplace transform toT
t
6
(
t
), I
t
6
(
t
), and V
t
6
(
t
), we obtain (16)
Ψ
T
t
6
t
=
0.1
s
+
0.397953
s
2
+
1.1857
s
3
+
3.53231
s
4
+
10.5191
s
5
+
31.3036
s
6
+
93.0202
s
7
,
Ψ
I
t
6
t
=
1
0
-
5
×
2.7
s
2
+
3.45473
s
3
-
5.04309
s
4
+
14.7535
s
5
-
34.0303
s
6
+
83.0376
s
7
,
Ψ
V
t
6
t
=
0.1
s
-
0.24
s
2
+
0.576081
s
3
-
1.38249
s
4
+
3.31783
s
5
-
7.96234
s
6
+
19.1086
s
7
.For simplicity, we let s
=
1
/
t. Then the above expressions become (17)
Ψ
T
t
6
t
=
0.1
t
+
0.397953
t
2
+
1.1857
t
3
+
3.53231
t
4
+
10.5191
t
5
+
31.3036
t
6
+
93.0202
t
7
,
Ψ
I
t
6
t
=
1
0
-
5
×
2.7
t
2
+
3.45473
t
3
-
5.04309
t
4
+
14.7535
t
5
-
34.0303
t
6
+
83.0376
t
7
,
Ψ
V
t
6
t
=
0.1
t
-
0.24
t
2
+
0.576081
t
3
-
1.38249
t
4
+
3.31783
t
5
-
7.96234
t
6
+
19.1086
t
7
.The highest indices of the Laplace transformations Ψ
(
T
t
6
(
t
)
), Ψ
(
I
t
6
(
t
)
), and Ψ
(
V
t
6
(
t
)
) are 7. Hence, Laplace transform increases the index by one for each component.To combine Laplace transform and Padé approximation, we adopt[
3
/
4
] Padé approximation to approximate Ψ
(
T
t
6
(
t
)
), Ψ
(
I
t
6
(
t
)
), and Ψ
(
V
t
6
(
t
)
). Because the highest indices of these power functions are 7, by Mathematica software again we have (18)
Ψ
T
t
6
t
P
=
0.1
t
-
0.521777
t
2
-
0.622784
t
3
1.0
-
9.1973
t
+
18.5161
t
2
+
0.0436746
t
3
-
0.0324062
t
4
,
Ψ
I
t
6
t
P
=
1
0
-
5
×
2.7
t
2
+
9.19088
t
3
1.0
+
2.1245
t
-
0.850547
t
2
-
0.407785
t
3
-
0.071839
t
4
,
Ψ
V
t
6
t
P
=
0.1
t
-
0.0321052
t
2
-
0.0160753
t
3
1.0
+
2.07895
t
-
0.932089
t
2
-
0.388528
t
3
+
0.000113259
t
4
.Recalling t
=
1
/
s, we obtain (19)
Ψ
T
t
6
t
P
=
0.1
s
3
-
0.521777
s
2
-
0.622784
s
1.0
s
4
-
9.1973
s
3
+
18.5161
s
2
+
0.0436746
s
-
0.0324062
,
Ψ
I
t
6
t
P
=
1
0
-
5
×
2.7
s
2
+
9.19088
s
1.0
s
4
+
2.1245
s
3
-
0.850547
s
2
-
0.407785
s
-
0.071839
,
Ψ
V
t
6
t
P
=
0.1
s
3
-
0.0321052
s
2
-
0.0160753
s
1.0
s
4
+
2.07895
s
3
-
0.932089
s
2
-
0.388528
s
+
0.000113259
.Using the inverse Laplace transform, we have(20)
T
L
P
L
t
≔
Ψ
-
1
Ψ
T
t
6
t
P
=
-
0.0161448
e
-
0.0425685
t
-
0.0174215
e
0.0410798
t
+
0.133575
e
2.97986
t
-
8.63837
×
1
0
-
6
e
6.21892
t
,
I
L
P
L
t
≔
Ψ
-
1
Ψ
I
t
6
t
P
=
1
0
-
5
×
0.420561
e
-
2.41214
t
+
3.48737
e
0.61872
t
-
2
e
-
0.165538
t
×
1.95397
cos
⁡
0.143987
t
-
3.15941
sin
⁡
0.143987
t
,
V
L
P
L
t
≔
Ψ
-
1
Ψ
V
t
6
t
P
=
0.100023
e
-
2.39987
t
-
0.0000742602
e
-
0.272939
t
+
0.0000120428
e
0.000291303
t
+
0.0000387879
e
0.593571
t
,where Ψ
-
1
(
·
) denotes the inverse Laplace transform. Till now, we have stated completely the posttreatment of ADM. In this way, we obtain the approximate solutions T
P
3, I
P
3, and V
P
3 by using the posttreatment Padé approximation and T
L
P
L
(
t
), I
L
P
L
(
t
), and V
L
P
L
(
t
) by using the posttreatment Padé approximation and Laplace transform, respectively.


## 4. Numerical Results

In this section, we demonstrate the effectiveness of the posttreatment technique, that is, the Padé approximation and the Laplace transform. For simplicity, if the Padé approximation is used to conduct the results of ADM, we call it ADM-Padé; if we use Padé approximation and Laplace transform to obtain the results of ADM, we call it ADM-LPL. We use the classical RK4 method, ADM, ADM-Padé, and ADM-LPL to solve HIV infection model (1). Numerical results of these methods are displayed by graphics.Figures1, 2, and 3 show the results of KR4, ADM-[3/3] Padé approximation, and ADM for functions T
(
t
), I
(
t
), and V
(
t
), respectively. From Figure 1, we see that the result of ADM is closer to that of RK4 than that of ADM-[3/3] Padé approximation when the time is growing. From Figures 2 and 3, we see that the results of ADM-[3/3] Padé approximation are closer to those of RK4 than those of ADM for I
(
t
) and V
(
t
) when the time is growing. From Figures 1, 2, and 3, we observe that Padé approximation is a better tool that can be used to obtain the results of ADM. This coincides with the theoretical analysis about Padé approximation. But, for HIV infection model (1), Figures 1, 2, and 3 show that the results of ADM-[3/3] Padé approximation are the same as those of ADM in a small time interval.Figure 1
Curves of KR4, ADM-[3/3] Padé approximation, and ADM forT
(
t
).Figure 2
Curves of KR4, ADM-[3/3] Padé approximation, and ADM forI
(
t
).Figure 3
Curves of KR4, ADM-[3/3] Padé approximation, and ADM forV
(
t
).Figures4, 5, and 6 show the results of KR4, ADM-[4/4] Padé approximation [4], and ADM-[3/3] Padé approximation for the functions T
(
t
), I
(
t
), and V
(
t
), respectively. We see that the results of ADM-[3/3] Padé approximation are closer to those of KR4 than those of ADM-[4/4] Padé approximation when the time is growing, for functions T
(
t
), I
(
t
), and V
(
t
) whatever. This coincides with the theoretical analysis in Section 3. Therefore, for HIV infection model (1), Figures 4, 5, and 6 show that the results of ADM-[3/3] Padé approximation are the same as those of ADM-[4/4] Padé approximation in a small time interval.Figure 4
Curves of KR4 and ADM-[3/3] and ADM-[4/4] Padé approximations forT
(
t
).Figure 5
Curves of KR4 and ADM-[3/3] and ADM-[4/4] Padé approximations forI
(
t
).Figure 6
Curves of KR4 and ADM-[3/3] and ADM-[4/4] Padé approximations forV
(
t
).Figures7, 8, and 9 show the results of KR4, ADM-[3/3] Padé approximation, and ADM-LPL for the functions T
(
t
), I
(
t
), and V
(
t
), respectively. Figure 7 shows that the result of ADM-LPL is almost the same as that of RK4 for function T
(
t
). And Figure 8 shows that ADM-[3/3] Padé approximation outperforms ADM-LPL when the time is growing. Because of the scalar of the longitudinal coordinates in this figure, there is no much error for I
(
t
). Figure 9 shows that the result of ADM-LPL is almost the same as that of RK4 for the function V
(
t
), too. Therefore, for HIV infection model (1), Figures 7, 8, and 9 show that the results of ADM-[3/3] Padé approximation are the same as those of ADM-LPL in a small time interval.Figure 7
Curves of KR4, ADM-[3/3] Padé approximation, and ADM-LPL forT
(
t
).Figure 8
Curves of KR4, ADM-[3/3] Padé approximation, and ADM-LPL forI
(
t
).Figure 9
Curves of KR4, ADM-[3/3] Padé approximation, and ADM-LPL forV
(
t
).In summary, the results of ADM, ADM-Padé, and ADM-LPL are the same when they are used to solve HIV infection model (1) in a small time interval. From the results of these methods, we find that ADM-LPL is the best method for solving HIV infection model (1). Generally speaking, we need to find good posttreatment to conduct the result of ADM, which makes the results better when the time is growing. Padé approximation and Laplace transform are such good tools as needed. For HIV infection model (1), our numerical experiments have shown the effectiveness of posttreatment by using Padé approximation and Laplace transform.



---

*Source: 584787-2015-09-30.xml*