/ / Article

Research Article | Open Access

Volume 2013 |Article ID 157956 | 11 pages | https://doi.org/10.1155/2013/157956

# Couple of the Variational Iteration Method and Legendre Wavelets for Nonlinear Partial Differential Equations

Revised27 Dec 2012
Accepted27 Dec 2012
Published26 Feb 2013

#### 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 vector-matrix 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) , homotopy perturbation method (HPM) , homotopy analysis method (HAM) [13, 14], variational iteration method (VIM) , and wavelets method .

The variational iteration method (VIM) proposed by He  has been shown to be very efficient for handling a wide class of physical problems [1618, 3041]. 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  or using Adomian’s Polynomials , 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  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  (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  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 th-order 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 well-known 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 two-dimensional 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 .

The derivative of the vector can be expressed by where is the operational matrix of derivative given by .

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 m-square 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 long-wave (RLW) equation : 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 fifth-order 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.

 0.25 0.50 0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 VIM 0.20002 0.40004 0.60006 0.80007 0.16702 0.33405 0.50107 0.66810 VIMLW 0.20002 0.40004 0.60006 0.80006 0.16702 0.33405 0.50107 0.66809 Exact 0.20000 0.40000 0.60000 0.80000 0.16667 0.33333 0.50000 0.66667 0.75 1.00 0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 VIM 0.14461 0.28921 0.43382 0.57842 0.12989 0.25977 0.38966 0.51955 VIMLW 0.14460 0.28921 0.43382 0.57841 0.12989 0.25977 0.38966 0.51953 Exact 0.14286 0.28571 0.42857 0.57143 0.12500 0.25000 0.37500 0.50000

Example 5. Consider the following equation : 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 fifth-order terms are used in evaluating the approximate solutions. The results obtained using the VIMLW are in good agreement with the results of VIM.

 0.25 0.50 0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 VIM 0.05001 0.20002 0.45004 0.80007 0.04176 0.16702 0.37581 0.66810 VIMLW 0.05000 0.20002 0.45004 0.80021 0.04168 0.16699 0.37572 0.66831 Exact 0.05000 0.20000 0.45000 0.80000 0.04167 0.16667 0.37500 0.66667 0.75 1.00 0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 VIM 0.03615 0.14461 0.32536 0.57842 0.03247 0.12989 0.29225 0.51955 VIMLW 0.03577 0.14443 0.32491 0.57751 0.03128 0.12932 0.29083 0.51421 Exact 0.03571 0.14286 0.32143 0.57143 0.03125 0.12500 0.28125 0.50000

Example 6. We consider the following equation : 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 fourth-order 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.

 0.25 0.50 0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 VIM 0.06250 0.12500 0.18750 0.25000 0.12508 0.25015 0.37523 0.50030 VIMLW 0.06250 0.12500 0.18750 0.25000 0.12508 0.25015 0.37523 0.50030 Exact 0.06250 0.12500 0.18750 0.25000 0.12500 0.25000 0.37500 0.50000 0.75 1.00 0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 VIM 0.18882 0.37764 0.56646 0.75527 0.25992 0.51983 0.77975 1.03970 VIMLW 0.18882 0.37764 0.56646 0.75527 0.25992 0.51983 0.77975 1.03970 Exact 0.18750 0.37500 0.56250 0.75000 0.25000 0.50000 0.75000 1.00000

Example 7. Consider the following Burgers-Poisson (BP) equation of the form : 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 fourth-order 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.

 0.25 0.50 0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 VIM 0.00009 0.20011 0.40013 0.60015 −0.16488 0.00215 0.16917 0.33620 VIMLW 0.00009 0.20011 0.40013 0.59864 −0.16489 0.00215 0.16920 0.33283 Exact 0.00000 0.20000 0.40000 0.60000 −0.16667 0.00000 0.16667 0.33333 0.75 1.00 0.25 0.50 0.75 1.00 0.25 0.50 0.75 1.00 VIM −0.27697 −0.13237 0.01224 0.15694 −0.35057 −0.22068 −0.09079 0.03909 VIMLW −0.27701 −0.13237 0.01228 0.15213 −0.35070 −0.22068 −0.09078 0.03422 Exact −0.28571 −0.14286 0.00000 0.14286 −0.37500 −0.25000 −0.12500 0.00000

#### 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.

1. 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
2. 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.
3. 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
4. Q. Esmaili, A. Ramiar, E. Alizadeh, and D. D. Ganji, “An approximation of the analytical solution of the Jeffery-Hamel flow by decomposition method,” Physics Letters A, vol. 372, no. 19, pp. 3434–3439, 2008.
5. 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.
6. S. Momani and Z. Odibat, “Analytical solution of a time-fractional Navier-Stokes equation by Adomian decomposition method,” Applied Mathematics and Computation, vol. 177, no. 2, pp. 488–494, 2006.
7. J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.
8. J.-H. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site | Google Scholar | MathSciNet
9. J.-H. He, “The homotopy perturbation method nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004.
10. S. T. Mohyud-Din and M. A. Noor, “Homotopy perturbation method for solving fourth-order boundary value problems,” Mathematical Problems in Engineering, vol. 2007, Article ID 98602, 15 pages, 2007.
11. S. T. Mohyud-Din and M. A. Noor, “Homotopy perturbation method for solving partial differential equations,” Zeitschrift für Naturforschung A, vol. 64, no. 3-4, pp. 157–170, 2009. View at: Google Scholar | Zentralblatt MATH
12. S. T. Mohyud-Din and M. A. Noor, “Homotopy perturbation method and Padé approximants for solving Flierl-Petviashivili equation,” Applications and Applied Mathematics, vol. 3, no. 2, pp. 224–234, 2008.
13. S. J. Liao, “An approximate solution technique not depending on small parameters: a special example,” International Journal of Non-Linear Mechanics, vol. 30, no. 3, pp. 371–380, 1995.
14. 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
15. J. H. He, “Variational iteration method—a kind of non-linear analytical technique: some examples,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 699–708, 1999. View at: Publisher Site | Google Scholar
16. J.-H. He, “Variational iteration method for autonomous ordinary differential systems,” Applied Mathematics and Computation, vol. 114, no. 2-3, pp. 115–123, 2000.
17. 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.
18. J.-H. He and X.-H. Wu, “Construction of solitary solution and compacton-like solution by variational iteration method,” Chaos, Solitons & Fractals, vol. 29, no. 1, pp. 108–113, 2006.
19. J.-H. He and X.-H. Wu, “Variational iteration method: new development and applications,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 881–894, 2007.
20. 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.
21. 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
22. J.-H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.
23. 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
24. 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.
25. Ü. Lepik, “Numerical solution of evolution equations by the Haar wavelet method,” Applied Mathematics and Computation, vol. 185, no. 1, pp. 695–704, 2007.
26. S. G. Venkatesh, S. K. Ayyaswamy, and S. Raja Balachandar, “The Legendre wavelet method for solving initial value problems of Bratu-type,” Computers & Mathematics with Applications, vol. 63, no. 8, pp. 1287–1295, 2012.
27. S. A. Yousefi, “Legendre wavelets method for solving differential equations of Lane-Emden type,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1417–1422, 2006.
28. R. K. Pandey, N. Kumar, A. Bhardwaj, and G. Dutta, “Solution of Lane-Emden type equations using Legendre operational matrix of differentiation,” Applied Mathematics and Computation, vol. 218, no. 14, pp. 7629–7637, 2012.
29. F. Yin, J. Song, F. Lu, and H. Leng, “A coupled method of Laplace transform and legendre wavelets for Lane-Emden-type differential equations,” Journal of Applied Mathematics, vol. 2012, Article ID 163821, 16 pages, 2012. View at: Publisher Site | Google Scholar
30. L. M. B. Assas, “Variational iteration method for solving coupled-KdV equations,” Chaos, Solitons and Fractals, vol. 38, no. 4, pp. 1225–1228, 2008.
31. 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.
32. 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.
33. E. M. Abulwafa, M. A. Abdou, and A. A. Mahmoud, “Nonlinear fluid flows in pipe-like domain problem using variational-iteration method,” Chaos, Solitons & Fractals, vol. 32, no. 4, pp. 1384–1397, 2007.
34. 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
35. 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.
36. 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.
37. S. T. Mohyud-Din, 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
38. 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.
39. 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. 7-8, pp. 1018–1027, 2007.
40. 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.
41. E. Hizel and S. K. Küçükarslan, “A numerical analysis of the Burgers-Poisson (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
42. S. T. Mohyud-Din and M. A. Noor, “Modified variational iteration method for solving Fisher's equations,” Journal of Applied Mathematics and Computing, vol. 31, no. 1-2, pp. 295–308, 2009.
43. M. A. Noor and S. T. Mohyud-Din, “Variational iteration method for solving higher-order 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
44. M. A. Noor and S. T. Mohyud-Din, “Modified variational iteration method for heat and wave-like equations,” Acta Applicandae Mathematicae, vol. 104, no. 3, pp. 257–269, 2008.
45. M. A. Noor and S. T. Mohyud-Din, “Modified variational iteration method for solving Helmholtz equations,” Computational Mathematics and Modeling, vol. 20, no. 1, pp. 40–50, 2009.
46. M. A. Noor and S. T. Mohyud-Din, “Variational iteration method for fifth-order boundary value problems using He's polynomials,” Mathematical Problems in Engineering, vol. 2008, Article ID 954794, 12 pages, 2008.
47. M. A. Noor and S. T. Mohyud-Din, “Modified variational iteration method for solving fourth-order boundary value problems,” Journal of Applied Mathematics and Computing, vol. 29, no. 1-2, pp. 81–94, 2009.
48. M. A. Noor and S. T. Mohyud-Din, “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
49. 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.
50. S. Abbasbandy, “Numerical solution of non-linear Klein-Gordon 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
51. S. T. Mohyud-Din 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
52. S. T. Mohyud-Din, 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
53. M. A. Noor and S. T. Mohyud-Din, “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.
54. M. A. Noor and S. T. Mohyud-Din, “Variational iteration decomposition method for solving eighth-order boundary value problems,” Differential Equations and Nonlinear Mechanics, vol. 2007, Article ID 19529, 16 pages, 2007.
55. M. M. Hosseini, “Adomian decomposition method with Chebyshev polynomials,” Applied Mathematics and Computation, vol. 175, no. 2, pp. 1685–1693, 2006.
56. 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
57. 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.
58. 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.
59. 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.

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at help@hindawi.com to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. 