Research Article  Open Access
Couple of the Variational Iteration Method and Legendre Wavelets for Nonlinear Partial Differential Equations
Abstract
This paper develops a modified variational iteration method coupled with the Legendre wavelets, which can be used for the efficient numerical solution of nonlinear partial differential equations (PDEs). The approximate solutions of PDEs are calculated in the form of a series whose components are computed by applying a recursive relation. Block pulse functions are used to calculate the Legendre wavelets coefficient matrices of the nonlinear terms. The main advantage of the new method is that it can avoid solving the nonlinear algebraic system and symbolic computation. Furthermore, the developed vectormatrix form makes it computationally efficient. The results show that the proposed method is very effective and easy to implement.
1. Introduction
Nonlinear phenomena are of fundamental importance in applied mathematics and physics and thus have attracted much attention. It is well known that most engineering problems are nonlinear, and it is very difficult to achieve the solution analytically or numerically. The analytical methods commonly used to solve them are very restricted, while the numerical techniques involving discretization of the variables on the other hand give rise to rounding off errors. Considerable attention has been paid to developing an efficient and fast convergent method. Recently, several approximate methods are introduced to find the numerical solutions of nonlinear PDEs, such as Adomian’s decomposition method (ADM) [1–6], homotopy perturbation method (HPM) [7–12], homotopy analysis method (HAM) [13, 14], variational iteration method (VIM) [15–23], and wavelets method [24–29].
The variational iteration method (VIM) proposed by He [15–23] has been shown to be very efficient for handling a wide class of physical problems [16–18, 30–41]. If the exact solution of the nonlinear PDEs exists, the VIM gives rapidly convergent successive approximations; otherwise, a few approximations can be used for numerical purposes. In order to improve the efficiency of these algorithms, several modifications, such as variational iteration method using He’s Polynomials [42–48] or using Adomian’s Polynomials [49–54], have been developed and successfully applied to various engineering problems. However, since the variational iteration method provides the solution as a sequence of iterates, its successive iterations may be very complex, so that the resulting integrations in its iterative relation may be impossible to perform analytically.
In recent years, wavelets have found their way into many different fields of science and engineering. Various wavelets [24–29] have been used for studying problems with greater computational complexity and proved to be powerful tools to explore a new direction in solving differential equations. Unlike the variational iteration method that requires symbolic computations, the wavelets method converts the PDE into algebraic equations by the operational matrices, which can be solved by an iterative procedure. It is worthy to mention here that the method based on operational matrices of an orthogonal function for solving differential equations is computer oriented. The problem with this approach is that the algebraic equations may be singular and nonlinear.
Recently, some efficient modifications of ADM (using [55, 56]) and VIM or HAM [57] (using Legendre polynomials) are presented to approximate nonhomogeneous terms in nonlinear differential equations. Motivated and inspired by the ongoing research in these areas, we implement Legendre wavelets within the framework of VIM to facilitate the computational work of the method while still keeping the accuracy. The remainder of the paper is organized as follows. Section 2 introduces the VIM. In Section 3, we describe the basic formulation of Legendre wavelets and the operational matrix required for our subsequent development. In Section 4, we propose a new variational iteration method using Legendre wavelets (VIMLW). In order to demonstrate the validity and applicability of VIMLW, four examples are given in Section 5. Finally, a brief summary is presented.
2. Variational Iteration Method
This section introduces the basic ideas of variational iteration method (VIM). Here a description of the method [15–23] is given to handle the general nonlinear problem: where is a linear operator, is a nonlinear operator, and is a known analytic function. According to He’s VIM, we can construct a correction functional as follows: where is a general Lagrange multiplier which can be optimally identified via variational theory and is a restricted variation which means . Therefore, the Lagrange multiplier should be first determined via integration by parts. The successive approximation () of the solution will be readily obtained by using the obtained Lagrange multiplier and any selective function . The zeroth approximation may select any function that just meets, at least, the initial and boundary conditions. With determined, several approximations , , follow immediately. Consequently, the exact solution may be obtained as The VIM depends on the proper selection of the initial approximation . Finally, we approximate the solution of the initial value problem (1) by the thorder term . It has been validated that VIM is capable of effectively, easily, and accurately solving a large class of nonlinear problems.
3. Legendre Wavelets
3.1. Legendre Wavelets
Legendre wavelets have four arguments: is any positive integer, (), is the order for Legendre polynomials, and is the normalized time. They are defined on the interval [0, 1) as follows: where , . The coefficient is for orthonormality, the dilation parameter is , and the translation parameter . Here, are the wellknown Legendre polynomials of order defined on the interval .
A function defined over may be expanded by Legendre wavelet series as with in (6); denotes the inner product.
If the infinite series in (5) is truncated, then it can be written as where and are matrices given by
A twodimensional function defined over may be expanded by Legendre wavelet series as with Equation (10) can be written into the discrete form (in matrix form) by where is a matrix given by
The integration and derivative operation matrices of the Legendre wavelets have been derived in [58, 59].
The integration of the vector defined in (9) can be obtained as where is a matrix given by [58].
The derivative of the vector can be expressed by where is the operational matrix of derivative given by [59].
The integration of with respect to variable can be expressed as
Similarly, the integration of with respect to variable can be expressed as
The derivative of with respect to variable can be expressed as
Similarly, the derivative of with respect to variable can be expressed as
3.2. Block Pulse Functions
The block pulse functions (BPFs) form a complete set of orthogonal functions that are defined on the interval by for . It is also known that for arbitrary absolutely integrable function on can be expanded in block pulse functions: in which where are the coefficients of the block pulse function given by The elementary properties of BPFs are as follows.(1) Disjointness: the BPFs are disjoined with each other in the interval :
for .(2) Orthogonality: the BPFs are orthogonal with each other in the interval :
for .(3) Completeness: the BPFs set is complete when approaches infinity. This means that for every , when approaches to the infinity, Parseval’s identity holds: where
Definition 1. Let and be two matrices of , then .
Lemma 2. Assuming that and are two absolutely integrable functions, which can be expanded in block pulse function as , and respectively, then one has where .
Proof. According to the disjointness property of BPFS in (16), we have
Lemma 3. Let and be two absolutely integrable functions, which can be expanded in block pulse function as and , respectively, one has where .
3.3. Nonlinear Term Approximation
The Legendre wavelets can be expanded into set of block pulse functions as Taking the collocation points as follow, The msquare Legendre matrix is defined as
The operational matrix of product of Legendre wavelets can be obtained by using the properties of BPFs. Let and be two absolutely integrable functions, which can be expanded in Legendre wavelets as and , respectively.
From (31), we have and let , , .
By employing Lemma 3, we get where .
4. Variational Iteration Method Using Legendre Wavelets
In this section, we present a new modification of variational iteration method using Legendre wavelets (called VIMLW). This algorithm can be implemented for solving nonlinear PDEs effectively.
To deduce the basic relations of our proposed algorithm, consider the following forms of initial value problems: where and are linear operator and nonlinear operator, respectively, and is a known analytic function, subject to the initial condition . It should be noted here that contains the term , where is a positive integer.
According to the traditional VIM, we can construct the correction functional for (36) as The Lagrange multiplier of (37) is
In order to improve the performance of VIM, we introduce Legendre wavelets to approximate and the nonhomogeneous term as
Now for the nonlinear part, by nonlinear term approximation described in Section 3.3, we have where is matrix of order .
For the linear part, we have where is a matrix of order .
Then the iteration formula (37) can be constructed as
If is constant, we have
When is a function of , the Legendre wavelets are used to approximate as Substituting (44) into (42), we have Since we get According to the property of block pulse functions, we obtain where
Substituting (48) into (45), we have where .
Finally, we get the iteration formula as follows:
5. Numerical Examples
To demonstrate the effectiveness and good accuracy of the VIMLW, four different examples will be examined.
Example 4. Consider the regularized longwave (RLW) equation [39]:
with the initial condition and the exact solution is .
By assuming and from (52), we have
where , .
We utilize the methods presented in this paper to solve (52) with and . Table 1 shows the approximate solutions for (52) obtained for different points using the variational iteration and VIMLW method. Figure 1 presents the Exact solution and VIMLW approximate solution of Example 4. Note that only the fifthorder term of their solutions is used in evaluating the approximate solutions for Example 4. We can see that the approximate solution obtained with VIMLW gives almost the same results as that with VIM. It indicates that the approximate solution is quite close to the exact one.

Example 5. Consider the following equation [39]:
with the initial conditions , and the exact solution is .
By assuming and from (54), we have
where , .
We employ the methods presented in this paper to solve (54) with and . The numerical results are presented in Table 2 and shown in Figure 2. It is to be noted that only the fifthorder terms are used in evaluating the approximate solutions. The results obtained using the VIMLW are in good agreement with the results of VIM.

Example 6. We consider the following equation [40]:
with the initial conditions and the exact solution is , where .
By assuming , , we have
where , .
Table 3 shows the approximate solutions for (56) with and using the VIM and the VIMLW methods and the results are plotted in Figure 3. It is to be noted that only the fourthorder terms of VIM and VIMLW are used in evaluating the approximate solutions in Table 3. We observe that the approximate solution of (56) with VIMLW gives analogous results to that obtained by VIM, which shows that the approximate solution remains closed form to the exact one.

Example 7. Consider the following BurgersPoisson (BP) equation of the form [41]:
with the initial conditions , and the exact solution is .
By assuming , we have
where .
And
where .
Table 4 shows the approximate solutions to (58) with and with VIM and VIMLW, and Figure 4 presents the Exact solution and VIMLW approximate solution of Example 7. Only the fourthorder terms are used in evaluating the approximate solutions in Table 4. From Table 4 and Figure 4 the approximate solution of the given Example 7 by using VIMLW is in good agreement with the results of VIM and it clearly appears that the approximate solution remains closed form to exact solution.

6. Conclusion
A new modification of variational iteration method using Legendre wavelets is proposed and employed to solve a number of nonlinear partial differential equations. The proposed method can give approximations of higher accuracy and closed form solutions if existed. There are four important points to make here. First, unlike the VIM, the VIMLW can easily overcome the difficulty arising in the evaluation integration and the derivative of nonlinear terms and does not need symbolic computation. Second, by using the properties of BPFs, operational matrices of product of Legendre wavelets are derived and utilized to deal with nonlinear terms. Third, compared with Legendre wavelets method, the VIMLW only needs a few iterations instead of solving a system of nonlinear algebraic equations. Fourth and most important, VIMLW is computer oriented and can use existing fast algorithms to reduce the computation cost.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (Grant no. 41105063). The authors are very grateful to the reviewers for carefully reading the paper and for thier comments and suggestions which have improved the paper.
References
 G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, vol. 60 of Fundamental Theories of Physics, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1994. View at: Zentralblatt MATH  MathSciNet
 G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Adomian, “Solutions of nonlinear P. D. E,” Applied Mathematics Letters, vol. 11, no. 3, pp. 121–123, 1998. View at: Publisher Site  Google Scholar  MathSciNet
 Q. Esmaili, A. Ramiar, E. Alizadeh, and D. D. Ganji, “An approximation of the analytical solution of the JefferyHamel flow by decomposition method,” Physics Letters A, vol. 372, no. 19, pp. 3434–3439, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A.M. Wazwaz, “A new algorithm for calculating Adomian polynomials for nonlinear operators,” Applied Mathematics and Computation, vol. 111, no. 1, pp. 53–69, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Momani and Z. Odibat, “Analytical solution of a timefractional NavierStokes equation by Adomian decomposition method,” Applied Mathematics and Computation, vol. 177, no. 2, pp. 488–494, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 34, pp. 257–262, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “A coupling method of a homotopy technique and a perturbation technique for nonlinear problems,” International Journal of NonLinear Mechanics, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 J.H. He, “The homotopy perturbation method nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. T. MohyudDin and M. A. Noor, “Homotopy perturbation method for solving fourthorder boundary value problems,” Mathematical Problems in Engineering, vol. 2007, Article ID 98602, 15 pages, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. T. MohyudDin and M. A. Noor, “Homotopy perturbation method for solving partial differential equations,” Zeitschrift für Naturforschung A, vol. 64, no. 34, pp. 157–170, 2009. View at: Google Scholar  Zentralblatt MATH
 S. T. MohyudDin and M. A. Noor, “Homotopy perturbation method and Padé approximants for solving FlierlPetviashivili equation,” Applications and Applied Mathematics, vol. 3, no. 2, pp. 224–234, 2008. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 S. J. Liao, “An approximate solution technique not depending on small parameters: a special example,” International Journal of NonLinear Mechanics, vol. 30, no. 3, pp. 371–380, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. J. Liao, “Boundary element method for general nonlinear differential operators,” Engineering Analysis with Boundary Elements, vol. 20, no. 2, pp. 91–99, 1997. View at: Google Scholar
 J. H. He, “Variational iteration method—a kind of nonlinear analytical technique: some examples,” International Journal of NonLinear Mechanics, vol. 34, no. 4, pp. 699–708, 1999. View at: Publisher Site  Google Scholar
 J.H. He, “Variational iteration method for autonomous ordinary differential systems,” Applied Mathematics and Computation, vol. 114, no. 23, pp. 115–123, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “Variational principles for some nonlinear partial differential equations with variable coefficients,” Chaos, Solitons & Fractals, vol. 19, no. 4, pp. 847–851, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He and X.H. Wu, “Construction of solitary solution and compactonlike solution by variational iteration method,” Chaos, Solitons & Fractals, vol. 29, no. 1, pp. 108–113, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He and X.H. Wu, “Variational iteration method: new development and applications,” Computers & Mathematics with Applications, vol. 54, no. 78, pp. 881–894, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J.H. He, “Variational iteration method—some recent results and new interpretations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 3–17, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. H. He, G. C. Wu, and F. Austin, “The variational iterational method which should be follow,” Nonlinear Science Letters A, vol. 1, no. 1, pp. 1–30, 2010. View at: Google Scholar
 J.H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. H. He, “Asymptotic methods for solitary solutions and compactons,” Abstract and Applied Analysis, vol. 2012, Article ID 916793, 130 pages, 2012. View at: Publisher Site  Google Scholar
 G. Hariharan, K. Kannan, and K. R. Sharma, “Haar wavelet method for solving Fisher's equation,” Applied Mathematics and Computation, vol. 211, no. 2, pp. 284–292, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Ü. Lepik, “Numerical solution of evolution equations by the Haar wavelet method,” Applied Mathematics and Computation, vol. 185, no. 1, pp. 695–704, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. G. Venkatesh, S. K. Ayyaswamy, and S. Raja Balachandar, “The Legendre wavelet method for solving initial value problems of Bratutype,” Computers & Mathematics with Applications, vol. 63, no. 8, pp. 1287–1295, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. A. Yousefi, “Legendre wavelets method for solving differential equations of LaneEmden type,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1417–1422, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. K. Pandey, N. Kumar, A. Bhardwaj, and G. Dutta, “Solution of LaneEmden type equations using Legendre operational matrix of differentiation,” Applied Mathematics and Computation, vol. 218, no. 14, pp. 7629–7637, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Yin, J. Song, F. Lu, and H. Leng, “A coupled method of Laplace transform and legendre wavelets for LaneEmdentype differential equations,” Journal of Applied Mathematics, vol. 2012, Article ID 163821, 16 pages, 2012. View at: Publisher Site  Google Scholar
 L. M. B. Assas, “Variational iteration method for solving coupledKdV equations,” Chaos, Solitons and Fractals, vol. 38, no. 4, pp. 1225–1228, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. M. Odibat, “Construction of solitary solutions for nonlinear dispersive equations by variational iteration method,” Physics Letters A, vol. 372, no. 22, pp. 4045–4052, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. M. Abulwafa, M. A. Abdou, and A. A. Mahmoud, “The solution of nonlinear coagulation problem with mass loss,” Chaos, Solitons & Fractals, vol. 29, no. 2, pp. 313–330, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. M. Abulwafa, M. A. Abdou, and A. A. Mahmoud, “Nonlinear fluid flows in pipelike domain problem using variationaliteration method,” Chaos, Solitons & Fractals, vol. 32, no. 4, pp. 1384–1397, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. Bildik and A. Konuralp, “The use of variational iteration method, differential transform method and adomian decomposition method for solving different types of nonlinear partial differential equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 7, no. 1, pp. 65–70, 2006. View at: Google Scholar
 S. Momani and S. Abuasad, “Application of He's variational iteration method to Helmholtz equation,” Chaos, Solitons & Fractals, vol. 27, no. 5, pp. 1119–1123, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Momani and Z. Odibat, “Numerical comparison of methods for solving linear differential equations of fractional order,” Chaos, Solitons & Fractals, vol. 31, no. 5, pp. 1248–1255, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. T. MohyudDin, M. A. Noor, K. I. Noor, and M. M. Hosseini, “Solution of singular equation by He's variational iteration method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 11, no. 2, pp. 81–86, 2010. View at: Google Scholar
 N. H. Sweilam and M. M. Khader, “Variational iteration method for one dimensional nonlinear thermoelasticity,” Chaos, Solitons & Fractals, vol. 32, no. 1, pp. 145–149, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 D. D. Ganji, H. Tari, and M. Bakhshi Jooybari, “Variational iteration method and homotopy perturbation method for nonlinear evolution equations,” Computers & Mathematics with Applications, vol. 54, no. 78, pp. 1018–1027, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Z. Odibat and S. Momani, “Numerical methods for nonlinear partial differential equations of fractional order,” Applied Mathematical Modelling, vol. 32, no. 1, pp. 28–39, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. Hizel and S. K. Küçükarslan, “A numerical analysis of the BurgersPoisson (BP) equation using variational iteration method,” in Proceedings of the 3rd WSEAS International Conference on Applied and Theoretical Mechanics, Tenerife, Spain, December 2007. View at: Google Scholar
 S. T. MohyudDin and M. A. Noor, “Modified variational iteration method for solving Fisher's equations,” Journal of Applied Mathematics and Computing, vol. 31, no. 12, pp. 295–308, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Noor and S. T. MohyudDin, “Variational iteration method for solving higherorder nonlinear boundary value problems using He's polynomials,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 9, no. 2, pp. 141–156, 2008. View at: Google Scholar
 M. A. Noor and S. T. MohyudDin, “Modified variational iteration method for heat and wavelike equations,” Acta Applicandae Mathematicae, vol. 104, no. 3, pp. 257–269, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Noor and S. T. MohyudDin, “Modified variational iteration method for solving Helmholtz equations,” Computational Mathematics and Modeling, vol. 20, no. 1, pp. 40–50, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Noor and S. T. MohyudDin, “Variational iteration method for fifthorder boundary value problems using He's polynomials,” Mathematical Problems in Engineering, vol. 2008, Article ID 954794, 12 pages, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Noor and S. T. MohyudDin, “Modified variational iteration method for solving fourthorder boundary value problems,” Journal of Applied Mathematics and Computing, vol. 29, no. 12, pp. 81–94, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Noor and S. T. MohyudDin, “Modified variational iteration method for Goursat and Laplace problems,” World Applied Sciences Journal, vol. 4, no. 4, pp. 487–498, 2008. View at: Google Scholar
 S. Abbasbandy, “A new application of He's variational iteration method for quadratic Riccati differential equation by using Adomian's polynomials,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 59–63, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Abbasbandy, “Numerical solution of nonlinear KleinGordon equations by variational iteration method,” International Journal for Numerical Methods in Engineering, vol. 70, no. 7, pp. 876–881, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 S. T. MohyudDin and M. A. Noor, “Solving Schrödinger equations by modified variational iteration method,” World Applied Sciences Journal, vol. 5, no. 3, pp. 352–357, 2008. View at: Google Scholar
 S. T. MohyudDin, M. A. Noor, and K. I. Noor, “Modified variational iteration method for solving Sine Gordon equations,” World Applied Sciences Journal, vol. 5, no. 3, pp. 352–357, 2008. View at: Google Scholar
 M. A. Noor and S. T. MohyudDin, “Solution of singular and nonsingular initial and boundary value problems by modified variational iteration method,” Mathematical Problems in Engineering, vol. 2008, Article ID 917407, 23 pages, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Noor and S. T. MohyudDin, “Variational iteration decomposition method for solving eighthorder boundary value problems,” Differential Equations and Nonlinear Mechanics, vol. 2007, Article ID 19529, 16 pages, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. M. Hosseini, “Adomian decomposition method with Chebyshev polynomials,” Applied Mathematics and Computation, vol. 175, no. 2, pp. 1685–1693, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. C. Tien and C. K. Chen, “Adomian decomposition method by Legendre polynomials,” Chaos, Solitons and Fractals, vol. 39, no. 5, pp. 2093–2101, 2009. View at: Publisher Site  Google Scholar
 Z. Odibat, “On Legendre polynomial approximation with the VIM or HAM for numerical treatment of nonlinear fractional differential equations,” Journal of Computational and Applied Mathematics, vol. 235, no. 9, pp. 2956–2968, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Razzaghi and S. Yousefi, “The Legendre wavelets operational matrix of integration,” International Journal of Systems Science, vol. 32, no. 4, pp. 495–502, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 F. Mohammadi and M. M. Hosseini, “A new Legendre wavelet operational matrix of derivative and its applications in solving the singular ordinary differential equations,” Journal of the Franklin Institute, vol. 348, no. 8, pp. 1787–1796, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
Copyright © 2013 Fukang Yin et al. 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.