Journal of Applied Mathematics

Volume 2013 (2013), Article ID 795015, 7 pages

http://dx.doi.org/10.1155/2013/795015

Research Article

A New Modification of Adomian Decomposition Method for Volterra Integral Equations of the Second Kind

Lie-jun Xie

Department of Mathematics, Faculty of Science, Ningbo University, Ningbo, Zhejiang 315211, China

Received 5 September 2013; Accepted 4 November 2013

Academic Editor: Roberto Barrio

Copyright © 2013 Lie-jun Xie. 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.

Abstract.

We propose a new modification of the Adomian decomposition method for Volterra integral equations of the second kind. By the Taylor expansion of the components apart from the zeroth term of the Adomian series solution, this new technology overcomes the problems arising from the previous decomposition method. The validity and applicability of the new technique are illustrated through several linear and nonlinear equations by comparing with the standard decomposition method and the modified decomposition method. The results obtained indicate that the new modification is effective and promising.

1. Introduction

The Adomian decomposition method (ADM) [1, 2] has been efficiently used to solve linear and nonlinear problems such as differential equations and integral equations. The method provides the solution as an infinite series in which each term can be easily determined. The rapid convergence of the series obtained by this method is thoroughly discussed by Cherruault et al. [3].

Recently, Wazwaz [4] proposed a reliable modified technique of ADM that accelerates the rapid convergence of decomposition series solution. The modified decomposition needs only a slight variation from the standard decomposition method. Although the modified decomposition method may provide the exact solution by using two iterations only and sometimes without using the so-called Adomian polynomials, its effectiveness is based on the assumption that the function can be divided into two parts, and thus the success of the modified method depends on the proper choice of and .

More recently, another new modification of decomposition method was proposed by Wazwaz and El-Sayed [5]. In the new modification, the process of dividing into two components was replaced by a Taylor series of infinite components. The ideas have been shown to be computationally efficient in applying the proposed technique in several differential and integral equations. However, as will be seen from the examples below, the modified decomposition method will sometimes encounter obstacles to obtain the exact solution. After that, several authors further proposed expressing the function by the orthogonal polynomial series, such as Chebyshev series [6], Legendre series [7, 8], and Jacobi and Gegenbauer series [9].

In this paper, we introduce a new modification of decomposition method and make further progress beyond the achievements made so far in this regard. Several examples concerning integral equations are tested, and the results suggest that this new idea proposes a promising tool for the computation of integral equations both linear and nonlinear.

2. Description of the New Modification

Since the beginning of the 1980s, the Adomian decomposition method has been applied to a wide class of integral equations. To illustrate the procedure, consider the following Volterra integral equations of the second kind given by

where the kernel and the function are given real-valued functions, is a parameter, and are linear and nonlinear operators of , respectively, and the unknown function is the solution to be determined.

In the decomposition method, we usually express the solution of (1) in a series form defined by

Moreover, the decomposition method identifies the nonlinear term by the decomposition series

where is the so-called Adomian polynomials, which can be evaluated by the following formula [10, 11]:

Substituting (2) and (3) into both sides of (1) gives

In the following we outline the basic feature of the standard decomposition method [1, 2] and the modified decomposition method [4, 5].

The Standard Adomian Decomposition Method (SADM). By the standard decomposition method, the components , , of the solution of (1), are completely determined in the following recurrence manner:

Having determined the components , , , the solution in a series form defined by (2) follows immediately.

The Modified Decomposition Method 1 (MADM1). The standard decomposition method by Adomian [1, 2] was modified by Wazwaz [4]. The modification is based on the assumption that the function can be divided into two parts, namely, and . Under this assumption we set

Accordingly, a slight variation was proposed only for the components and . The suggestion was that only the part will be assigned to the zeroth component , whereas the remaining part will be combined with the other terms given into (6b) to define . Consequently, the modified recursive relation

was developed.

It is important to note here that there were some conclusions made in [4]. First, the slight variation in reducing the number of terms of will result in a reduction of the computational work and will accelerate the convergence. Second, this slight variation in the definition of the components and may provide the solution by using two iterations only. Third, there is no need sometimes to evaluate the so-called Adomian polynomials required for the nonlinear equations.

The Modified Decomposition Method 2 (MADM2). In [5], it is pointed that sometimes it may be useful to express in Taylor series as follows:

Meanwhile, a new recursive relationship in the following form can be obtained:

According to (10), the terms of the solution follow immediately, and the solution can be obtained using (2).

It is evident that algorithm (10) reduces the number of terms involved in each component, and hence the size of calculations is minimized compared to the standard Adomian decomposition method. Moreover, this reduction of terms in each component facilitates the construction of Adomian polynomials for nonlinear operators [5].

The New Modified Decomposition Method. As indicated earlier, MADM1 may provide the exact solution by using two iterations only and sometimes without any need for Adomian polynomials, but its effectiveness depends on the proper choice of and , which may need quite a little computational work. Regarding SADM and MADM2, sometimes one will inevitably encounter the difficulties as will be seen from the examples below. In this new modification, we obtain the first two terms by

which are the same as the results by SADM. However, sometimes is too complicated to continue or is analytically impossible. Here, we suggest that may be expressed in Taylor series form as

next, we set

where the Adomian polynomial can be evaluated by (4) with and , which will be obtained by (11a) and (12), respectively. Also, we suggest that the same series expansion can be acted upon , that is,

Similarly, we set

and their Taylor expansions given by

By combining (11a), (12), (14), and (16), the solution can be obtained using the series (2).

In practice, for some problems, we have to determine a few terms in the series such as by truncating the series at certain term, regardless of the standard decomposition method and the modified decomposition method. Because of the uniformly convergence property of the infinite series, a few terms will attain the maximum accuracy.

To give a clear overview of this study, we have chosen several linear and nonlinear Volterra integral equations to testify the validity and applicability of the new modification. For comparison, examples will also be discussed by using SADM, MADM1, and MADM2.

3. Application for Linear Integral Equations

Example 1. Consider the linear Volterra integral equation [12]

Following suggestions (6a) and (6b) of SADM, the second component of the solution cannot be obtained because of the impossible integration, which also appears by using MADM2. If we turn to MADM1, the six different choices of and need quite a little computational work only if we set and directly. In the following part, we solve (17) by the new modification with the help of computer algebraic system Maple 12.

According to suggestions (11a) and (11b), we obtain

We then set the Taylor expression for by

And also, we set the Taylor expansion for , because the evaluation of integrals, obtained by applying (13) and (15) in this example, is analytically impossible. In this way, we have

Thus, the solution in a series form is given by

Combining (21) with the Taylor expansion of

the solution in a closed form is given by

Example 2. Consider the linear Volterra integral equation [5, 13]

In this example, both SADM and MADM2 can fulfil the exact solution effectively. Whereas, MADM1 needs quite a little computational work to obtain the exact solution only if we set and directly. In the following part, we solve (24) by the new modification.

Proceeding as before, we set

We then set the Taylor expression for by

and also

The solution in a series form is given by

and in a closed form by

Example 3. Consider the system of linear Volterra integral equations [12]

In this example, it is both difficult to obtain the exact solution by using SADM and MADM2. And also, like the problem mentioned in preceding examples, MADM1 needs much computational work to obtain the exact solution only if we set and for the first equation and and for the second equation, directly. In the following parts, we solve (30) by the new modification.

Proceeding as before, the system (30) is transformed into a set of recursive relations given by

We then set the Taylor expansion for and by

Similarly,

The solution in a series form is given by

and in a closed form by

4. Application for Nonlinear Integral Equations

Example 4. Consider the nonlinear Volterra integral equation [12, 14]

Following suggestions (6a) and (6b) of , the second component described via the sum of twenty-seven terms is too complicated for no confidence to continue. Meanwhile, there does not exist noise terms [15], which leads to the result quickly. By MADM1, the six different choices of and need quite a little computational work as the examples given before, unless we set and directly. The solution can be obtained by MADM2 exactly. In the following part, we solve (36) by the new modification.

According to suggestions (11a) and (11b), we obtain

We then set the Taylor expression for by

And also, we set the Taylor expansion for , because of the complicated integration. In this way, we have

where the Adomian polynomials for nonlinear term were calculated by using (4). The solution in a series form is given by

Combining (40) with the Taylor expansion of

the solution in a closed form is given by

Example 5. Consider the nonlinear Volterra integrodifferential equation of the first kind [12]

By integrating the second integral by parts and using the initial conditions, (43) can be converted to Volterra integral equation of the second kind given by

Following suggestions (6a) and (6b) of SADM, consists of fourteen terms and twenty-nine terms. It is so complicated that we can give up. The similar obstacle derived from the selections of and will be encountered if we use MADM1. By using MADM2, it is of certain difficulty to obtain the exact solution through the first finite terms of the series solution. In the following part, we solve (43) by the new modification.

Following as before, we set

where Adomian polynomials for nonlinear term were calculated by using (4). The solution in a series form is given by

Combining (46) with the Taylor expansion of

the solution in a closed form is given by

5. Discussion and Conclusion

In this work, we have carefully developed a new modification of Adomian decomposition method based on the Taylor expansion of the components apart from the zeroth term for Volterra integral equations. The new technology overcomes such obstacles as the complicated integrals and the impossible integrals analytically, which arise from either SADM or MADM2. If the exact solution exists in the zeroth component, MADM1 needs extensive workload for the suitable choice of and , while the new modification can proceed in a fixed manner. Therefore, it is appropriate for one to implement through the mathematical software Maple without manual tedious calculations.

The idea has been shown to be computationally efficient in applying the proposed technique to several linear and nonlinear Volterra integral equations and integrodifferential equations. Also, through the examples, it is seen that the new modification suits for those integral equations where the exact solution appears as part of the zeroth term. It is also interesting to point out that one can obtain the approximate solution of integral equations as mentioned in [16, 17] simply by slightly varying the new modification of this work.

Acknowledgment

The work was supported by the Scientific Research Fund of Zhejiang Provincial Education Department of China (no. Y201120655), K.C. Wong Magna Fund in Ningbo University.

References

G. Adomian, “A review of the decomposition method and some recent results for nonlinear equations,” Mathematical and Computer Modelling, vol. 13, no. 7, pp. 17–43, 1990.

G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, vol. 60 of Fundamental Theories of Physics, Kluwer Academic, Dordrecht, The Netherlands, 1994.

Y. Cherruault, G. Saccomandi, and B. Some, “New results for convergence of Adomian's method applied to integral equations,” Mathematical and Computer Modelling, vol. 16, no. 2, pp. 85–93, 1992.

A.-M. Wazwaz, “A reliable modification of Adomian decomposition method,” Applied Mathematics and Computation, vol. 102, no. 1, pp. 77–86, 1999.

A.-M. Wazwaz and S. M. El-Sayed, “A new modification of the Adomian decomposition method for linear and nonlinear operators,” Applied Mathematics and Computation, vol. 122, no. 3, pp. 393–405, 2001.

M. M. Hosseini, “Adomian decomposition method with Chebyshev polynomials,” Applied Mathematics and Computation, vol. 175, no. 2, pp. 1685–1693, 2006.

Y. Liu, “Adomian decomposition method with orthogonal polynomials: Legendre polynomials,” Mathematical and Computer Modelling, vol. 49, no. 5-6, pp. 1268–1273, 2009.

W.-C. Tien and C.-K. Chen, “Adomian decomposition method by Legendre polynomials,” Chaos, Solitons & Fractals, vol. 39, no. 5, pp. 2093–2101, 2009.

Y. Çenesiz and A. Kurnaz, “Adomian decomposition method by Gegenbauer and Jacobi polynomials,” International Journal of Computer Mathematics, vol. 88, no. 17, pp. 3666–3676, 2011.

A.-M. Wazwaz, “A new algorithm for calculating Adomian polynomials for nonlinear operators,” Applied Mathematics and Computation, vol. 111, no. 1, pp. 33–51, 2000.

E. Babolian and Sh. Javadi, “New method for calculating Adomian polynomials,” Applied Mathematics and Computation, vol. 153, no. 1, pp. 253–259, 2004.

A.-M. Wazwaz, Linear and Nonlinear Integral Equations: Methods and Applications, Higher Education Press, Beijing, China, 2011.

A.-M. Wazwaz, A First Course in Integral Equations, World Scientific Publishing, River Edge, NJ, USA, 1997.

M. Rahman, Integral Equations and Their Applications, WIT Press, Southampton, UK, 2007.

A.-M. Wazwaz and A. M. Wazwaz, “Necessary conditions for the appearance of noise terms in decomposition solution series,” Applied Mathematics and Computation, vol. 81, no. 2-3, pp. 265–274, 1997.

E. Babolian and A. Davari, “Numerical implementation of Adomian decomposition method,” Applied Mathematics and Computation, vol. 153, no. 1, pp. 301–305, 2004.

E. Babolian and A. Davari, “Numerical implementation of Adomian decomposition method for linear Volterra integral equations of the second kind,” Applied Mathematics and Computation, vol. 165, no. 1, pp. 223–227, 2005.