# A Comparison between Adomian’s Polynomials and He’s Polynomials for Nonlinear Functional Equations

Authors: Hossein Jafari; Saber Ghasempoor; Chaudry Masood Khalique
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

We will compare the standard Adomian decomposition method and the homotopy perturbation method applied to obtain the solution of nonlinear functional equations. We prove analytically that the two methods are equivalent for solving nonlinear functional equations. In Ghorbani (2009), Ghorbani presented a new definition which he called as He’s polynomials. In this paper, we also show that He’s polynomials are only the Adomian polynomials.


---

## Body

## 1. Introduction

The Adomian decomposition method (ADM) and the homotopy perturbation method (HPM) are two powerful methods which consider the approximate solution of a nonlinear equation as an infinite series usually converging to the accurate solution. These methods have been used in obtaining analytic and approximate solutions to a wide class of linear and nonlinear, differential, and integral equations.Öziş and Yıldırım compared Adomian’s method and He’s homotopy perturbation method [1] for solving certain nonlinear problems. Li also has shown that the ADM and HPM for solving nonlinear equations are equivalent [2]. In [3], Ghorbani has presented a definition which he called it as He’s polynomials.Consider the following nonlinear functional equation:(1)u=f+N(u),
where N is a nonlinear operator from Hilbert space H to H, u is an unknown function, and f is a known function in H. We are looking for a solution u of (1) belonging to H. We will suppose that (1) admits a unique solution. If (1) does not possess a unique solution, the ADM and HPM will give a solution among many (possible) other solutions. However, relatively few papers deal with the comparison of these methods with other existing techniques. In [4], a useful comparison between the decomposition method and the perturbation method showed the efficiency of the decomposition method compared to the tedious work required by the perturbation techniques. In [5], the advantage of the decomposition method over the Picard’s method has been emphasized. Sadat has shown that the Adomian decomposition method and perturbation method are closely related and lead to the same solution in many heat conduction problems [6]. In [7, 8] the HPM has compared with Liao’s homotopy analysis method and showed the HPM is special case of HAM, and the advantage of the HAM over the HPM has been emphasized.In this paper, we want to prove that He’s polynomials are only Adomian’s polynomials. We will also show that the standard Adomian decomposition method and the standard HPM are equivalent when applied for solving nonlinear functional equations.


## 2. Adomian’s Decomposition Method (ADM)

Let us consider the nonlinear equation (1) which can be written in the following canonical form:
(2)u=f+N(u).
The standard ADM consists of representing the solution of (1) as a series
(3)u(x)=∑i=0∞‍ui(x),
and the nonlinear function as the decomposed form:
(4)N(u(x))=∑i=0∞‍Ai,
where An, n=0,1,2,… are the Adomian polynomials of u0,u1,…,un given by [9, 10]
(5)An=1n!dndpn[N(∑i=0n‍uipi)]p=0.
Substituting (3) and (4) into (1) yields
(6)∑i=0∞‍ui(x)=f+∑i=0∞‍Ai.
The convergence of the series in (6) gives the desired relation
(7)u0=f,un+1=An,n=0,1,2,….
It should be pointed out that A0 depends only on u0, A1 depends only on u0 and u1; A2 depends only on u0, u1, and u2, and so on. The Adomian technique is very simple in its principles. The difficulties consist in proving the convergence of the introduced series.


## 3. Homotopy Perturbation Method (HPM)

This is a basic idea of homotopy method which is to continuously deform a simple problem easy to solve into the difficult problem under study.In this section, we apply the homotopy perturbation method [11–13] to the discussed problem. To illustrate the homotopy perturbation method (HPM), we consider (1) as
(8)L(v)=v(x)-f(x)-N(v)=0,
with solution u(x). The basic idea of the HPM is to construct a homotopy H(v;p):R×[0,1]→R which satisfies
(9)ℋ(v;p)=(1-p)F(v)+pL(v)=0,
where F(v) is a proper function with known solution which can be obtained easily. The embedding parameter p monotonically increases from 0 to 1 as the trivial problem F(v)=0 is continuously transformed to the original problem v-f-N(v)=0. From ℋ(v;p)=0, we have H(v;0)=F(v)=0 and H(v;1)=v-f-N(v)=0.It is better to takeF(v) as a deformation of L(v). For example, in (9), F(v)=v-f(x). By selecting F(v)=v-f(x) we can define another convex homotopy ℋ(v;p) by
(10)ℋ(v;p)=v(x)-f(x)-pN(v)=0.
The embedding parameter p∈(0,1] can be considered as an expanding parameter [14, 15]. The HPM uses the embedding parameter p as a “small parameter,” and writes the solution of (10) as a power series of p, that is,
(11)v=v0+v1p+v2p2+⋯.
Setting p=1 results in the approximate solution of (10):
(12)u=limp→1v=v0+v1+v2+⋯.
Substituting (11) into (10) and equating the terms with identical powers of p, we can obtain a series of equations of the following form:
(13)p0:v0-f(x)=0,p1:v1-H(v0)=0,p2:v2-H(v0,v1)=0,p3:v3-H(v0,v1,v2)=0,⋮,
where H(v0,v1,…,vj) depend upon v0,v1,…,vj. In view of (10) to determine H(v0,v1,…,vj), we use [16]
(14)H(v0,v1,…,vj)=1j!∂j∂pjN(∑i=0j‍vipi)|p=0.
It is obvious that the system of nonlinear equations in (13) is easy to solve, and the components vi, i≥0 of the homotopy perturbation method can be completely determined, and the series solutions are thus entirely determined. For the convergence of the previous method we refer the reader to the work of He [12, 17, 18].


## 4. Equivalence between ADM and HPM

In this section, we prove that the HPM and the ADM give same solution for solving nonlinear functional equations. We also show that the He polynomials are like the Adomian polynomials. In [3], Ghorbani has presented the following definition.Definition 1 (see [3]).
The He polynomials are defined as follows:(15)Hn(v0,…,vn)=1n!∂n∂pnN(∑i=0n‍vipi)|p=0,n=0,1,2,….Note 1.
Comparison between (5) and (15) has shown that the He polynomials are only Adomian’s polynomials, and it is calculated like Adomian’s polynomials.Theorem 2.
Suppose that nonlinear functionN(u) and the parameterized representation of v are v(p)=∑i=0∞‍vipi, where p is a parameter, then we have
(16)∂nN(v(p))∂pn|p=0=∂nN(∑i=0∞vipi)∂pn|p=0=∂nN(∑i=0nvipi)∂pn|p=0.Proof (see [3, 19]).
In Theorem3, we prove that the He polynomials are the Adomian polynomials.Theorem 3.
The He polynomials which are given by (15) are the Adomian polynomials.Proof.
From Taylor’s expansion ofN(v), we have
(17)N(v)=N(v0)+N′(v0)(v-v0)+12!N′′(v0)(v-v0)2+⋯
substituting (11) in (17) and expanding it in terms of p leads to
(18)N(∑i=0∞‍vipi)=N(v0)+N′(v0)(v1p+v2p2+⋯)+12!N′′(v0)(v1p+v2p2+⋯)2+⋯=N(v0)+N′(v0)v1p+(N′(v0)v2+12!N′′(v0)v12)p2+⋯=H0+H1p+H2p2+⋯,
where Hi, i=0,1,2,…, depends only on v0,v1,…,vi.
In order to obtainHn, we give n-order derivative of both sides of (18) with respect to p and let p=0, that is,
(19)∂nN(v(p))∂pn|p=0=∂n∑i=0∞Hipi∂pn|p=0.
According to Theorem 2(20)∂nN(∑i=0∞vipi)∂pn|p=0=∂nN(∑i=0nvipi)∂pn|p=0,∂n∑i=0∞Hipi∂pn|p=0=∂n∑i=0nHipi∂pn|p=0=n!Hn.
We know that Hi just depends on v0,v1,…,vi so (∂n∑i=0n‍Hipi)/∂pn|p=0=n!Hn. Substituting (20) in (19) leads us to find Hi in the following form:
(21)Hn=1n!∂nN(∑i=0nvipi)∂pn|p=0which is called for the first time by Ghorbani as the He polynomials [3]!Theorem 4.
The homotopy perturbation method for solving nonlinear functional equations is the Adomian decomposition method with the homotopyℋ(v;p) given by
(22)ℋ(v;p)=v-f(x)-pN(v).Proof.
Substituting (11) and (18) into (10) and equating the terms with the identical powers of p, we have
(23)ℋ(v;p)=∑i=0∞‍vipi-f(x)-p∑i=0∞‍Hipi=0,ℋ(v;p)=v0-f(x)+∑i=0∞‍(vi+1-Hi)pi+1=0,(24)p0:v0-f(x)=0,pn+1:vn+1-Hn=0,n=0,1,2,…
From (24) we have
(25)v0=f(x),vn+1=Hn,n=0,1,2,…
According to Theorem 3 we have Hn=An. Substituting (25) in (11) leads us to
(26)v=v0+v1p+v2p2+⋯=f(x)+A0p+A1p2+⋯,
so
(27)limp→1v=f(x)+A1+A2+⋯=f(x)+∑i=0∞‍Ai=∑i=0∞‍ui=u.
Therefore, by letting
(28)ℋ(v;p)=v-f(x)-pN(v),
we observe that the power series v0+v1p+v2p2+⋯ corresponds to the solution of the equation ℋ(v;p)=v-f(x)-pN(v)=0 and becomes the approximate solution of (1) if p→1. This shows that the homotopy perturbation method is the Adomian decomposition method with the homotopy ℋ(v;p) given by (28). The proof of Theorem 4 is completed.These two approaches give the same equations for high-order approximations. This is mainly because the Taylor series of a given function is unique, which is a basic theory in calculus. Thus, nothing is new in Ghorbani’s definition, except the new name “He’s polynomial.” He just employed the early ideas of ADM.Example 5 (see [20]).
Consider the following nonlinear Volterra integral equation:(29)y(x)=x+∫0x‍y2(t)dt,
with the exact solution y(x)=tanx.We apply standard ADM and HPM. For applying standard HPM, we construct following homotopy:(30)ℋ(u;p)=u(x)-x-p∫0x‍[u(t)]2dt=0.
In view of (13), we have
(31)p0:v0(x)-x=0,pn:vn+1(x)-∫0x‍H(v0,v1,…,vn)dt=0,n≥0.
Now if we apply ADM for solving (29), substituting (3) and (4) in (29) leads to
(32)∑i=0∞‍ui(x)=x+∫0x‍∑i=0∞‍Aidt.
In view of (7), we have following recursive formula:
(33)u0(x)=x,un+1(x)=∫0x‍Andtn≥0.
According to Theorem 3, we have An=H(v0,v1,…,vn). By solving (31) and (33), we have
(34)u(x)=∑i=0∞‍ui(x)=limp→1v0+v1p+v2p2+⋯=x+x33+2x515+17x7315+62x92835+⋯=tanx.


## 5. Conclusion

It has been shown that the standard HPM provides exactly the same solutions as the standard Adomian decomposition method for solving functional equations. It has been proved that He’s polynomials are only Adomian’s polynomials with different name.



---

*Source: 943232-2013-06-24.xml*