`Mathematical Problems in EngineeringVolume 2008, Article ID 954794, 11 pageshttp://dx.doi.org/10.1155/2008/954794`
Research Article

## Variational Iteration Method for Fifth-Order Boundary Value Problems Using He's Polynomials

Department of Mathematics, COMSATS Institute of Information Technology, Islamabad 44000, Pakistan

Received 8 January 2008; Revised 31 January 2008; Accepted 5 March 2008

Copyright © 2008 Muhammad Aslam Noor and Syed Tauseef Mohyud-Din. 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 apply the variational iteration method using He's polynomials (VIMHP) for solving the fifth-order boundary value problems. The proposed method is an elegant combination of variational iteration and the homotopy perturbation methods and is mainly due to Ghorbani (2007). The suggested algorithm is quite efficient and is practically well suited for use in these problems. The proposed iterative scheme finds the solution without any discritization, linearization, or restrictive assumptions. Several examples are given to verify the reliability and efficiency of the method. The fact that the proposed technique solves nonlinear problems without using Adomian's polynomials can be considered as a clear advantage of this algorithm over the decomposition method.

#### 1. Introduction

In this paper, we general fifth-order boundary value problem with boundary conditions where and are real and as many times differentiable as required for and , and , are real finite constants. This type of boundary value problems arises in the mathematical modeling of the viscoelastic flows and other branches of mathematical, physical, and engineering sciences, see [17] and the references therein. Several numerical methods including spectral Galerkin and collocation, decomposition, and sixth-order B-spline have been developed for solving fifth-order boundary value problems, see [1, 7] and the references therein. The use of spline function in the context of fifth-order boundary value problems was studied by Fyfe [3], who used the quintic polynomial spline functions to develop consistency relation connecting the values of solution with fifth-order derivative at the respective nodal points. He [818] developed the variational iteration method and homotopy perturbation method for solving linear, nonlinear, initial, and boundary value problems. It is worth mentioning that the origin of variational iteration method can be traced back to Inokuti et al. [19], but the real potential of this technique was explored by He [1318]. Moreover, He realized the physical significance of the variational iteration method, its compatibility with the physical problems and applied this promising technique to a wide class of linear and nonlinear, ordinary, partial, deterministic, or stochastic differential equation; see [1318]. The homotopy perturbation method [812, 17] was also developed by He by merging two techniques: the standard homotopy and the perturbation. The homotopy perturbation method was formulated by taking the full advantage of the standard homotopy and perturbation methods. The variational iteration method and homotopy perturbation method have been applied to a wide class of functional equations; see [4, 5, 836] and the references therein. In these methods, the solution is given in an infinite series usually converging to an accurate solution, see [4, 5, 822, 2536] and the references therein. In a later work, Ghorbani [23, 24] splits the nonlinear term into a series of polynomials calling them as He’s polynomials. Recently, Noor and Mohyud-Din used homotopy perturbation, variational iteration, and the iterative methods [46] for solving the fifth-order boundary value problems. The results are very encouraging and reveal the complete reliability of the new algorithm.

Inspired and motivated by the ongoing research in this area, we use the variational iteration method coupled with He’s polynomials for solving the fifth-order boundary value problems in this paper. It is worth mentioning that the proposed method is an elegant combination of variational iteration and the homotopy perturbation methods and is mainly due to Ghorbani [23, 24]. The use of He’s polynomials in the nonlinear term was first introduced by Ghorbani, see [23, 24]. The proposed algorithm provides the solution in a rapid convergent series which may lead the solution in a closed form. In this technique, the correction functional is developed [1322] and the Lagrange multipliers are calculated optimally via variational theory. The use of Lagrange multipliers reduces the successive application of the integral operator and the cumbersome of huge computational work, while still maintaining a very high level of accuracy. Finally, He’s polynomials are introduced in the correction functional and the comparison of like powers of gives solutions of various orders. The proposed iterative scheme takes full advantage of variational iteration method and the homotopy perturbation method. It is worth mentioning that the suggested method is applied without any discritization, restrictive assumption, or transformation and is free from round off errors. Unlike the method of separation of variables that requires initial and boundary conditions, the method provides an analytical solution by using the initial conditions only. The proposed method works efficiently and the results so far are very encouraging and reliable. The fact that the proposed VIMHP solves nonlinear problems without using Adomian’s polynomials can be considered as a clear advantage of this method over the decomposition method. Several examples are given to verify the reliability and efficiency of the algorithm.

#### 2. Variational Iteration Method

To illustrate the basic concept of the technique, we consider the following general differential equation: where is a linear operator, a nonlinear operator, and is the forcing term. According to variational iteration method [5, 1322, 25, 3236], we can construct a correct functional as follows: where is a Lagrange multiplier [1319], which can be identified optimally via variational iteration method. The subscripts n denote the nth approximation, is considered as a restricted variation, that is, and (2.2) is called as a correct functional. The solution of the linear problems can be solved in a single iteration step due to the exact identification of the Lagrange multiplier. The principles of variational iteration method and its applicability for various kinds of differential equations are given in [1319]. In this method, it is required first to determine the Lagrange multiplier optimally. The successive approximation , of the solution u will be readily obtained upon using the determined Lagrange multiplier and any selective function consequently, the solution is given by For the convergence and error estimates of variational iteration method, see Ramos [35].

#### 3. Homotopy Perturbation Method

To explain the homotopy perturbation method, we consider a general equation of the type where is an integral or differential operator. We define a convex homotopy by where is a functional operator with known solutions , which can be obtained easily. It is clear that for we have This shows that continuously traces an implicitly defined curve from a starting point to a solution function . The embedding parameter monotonically increases from zero to unit as the trivial problem is continuously deforms the original problem . The embedding parameter p (0, 1] can be considered as an expanding parameter [3, 811, 23, 24]. The homotopy perturbation method uses the homotopy parameter p as an expanding parameter [4, 812, 23, 24, 2731] to obtain If , then (3.5) corresponds to (3.2) and becomes the approximate solution of the form It is well known that series (3.6) is convergent for most of the cases and also the rate of convergence depends upon ; see [4, 812, 23, 24, 2731]. We assume that (3.6) has a unique solution. The comparisons of like powers of p give solutions of various orders.

#### 4. Variational Iteration Method Using He’s Polynomials (VIMHP)

To illustrate the basic concept of the variational homotopy perturbation method, we consider the following general differential equation: where is a linear operator, a nonlinear operator, and is the forcing term. According to variational iteration method [5, 6, 1322, 25, 26, 3236], we can construct a correct functional as follows: where is a Lagrange multiplier [1319], which can be identified optimally via variational iteration method. The subscripts n denote the nth approximation, is considered as a restricted variation, that is, and (4.2) is called as a correct functional. Now, we apply the homotopy perturbation method

which is the coupling of variational iteration and He’s polynomials. The comparison of like powers of p gives solutions of various orders.

#### 5. Numerical Applications

In this section, we apply the variational iteration method using He’s polynomials (VIMHP) for solving the fifth-order boundary value problems. We develop the correct functional and calculate the Lagrange multipliers optimally via variational theory, which reduces the successive application of the integral operator. The selection of initial value is done carefully because the approximants are heavily dependant upon initial value. He’s polynomials are introduced in the correct functional and finally, the comparison of like powers of p gives solutions of various orders. Numerical results are very encouraging. For the sake of comparison, we take the same examples as used in [47].

Example 5.1 (See [47]). Consider the following nonlinear boundary value problem of fifth-order, with boundary conditions The exact solution for this problem is The correct functional for the boundary value problem (5.1) and (5.2) is given as Making the correct functional stationary, using as the Lagrange multiplier [1318, 36], we get the following iterative formula: where

Applying the variational iteration method using He’s polynomials, Comparing the coefficient of like powers of p, consequently, we obtain the following approximants:

The series solution is given as Imposing the boundary conditions at and using leads to the following system: The solution of the above system gives Consequently, the series solution is given as which is in full agreement with [47].

Table 1 shows the exact values and the errors obtained by using the homotopy perturbation method (HPM) [4], variational iteration method (VIM) [5], decomposition method (ADM) [7], the sixth degree B-spline method [1], iterative method (ITM) [6], and the variational iteration method using He’s polynomials (VIMHP) for The table clearly indicates the improvements as compared with B-spline method. Higher accuracy can be obtained by evaluating more components of .

Table 1: Error estimates.

Remark 5.2. The numerical results clearly indicate that the results obtained by HPM, VIM, ADM, ITM, and the proposed VIMHP are the same. Moreover, it shows the improvements as compare to B-spline method.

Example 5.3 (See [47]). Consider the following linear boundary value problem of fifth-order: with boundary conditions The exact solution of the problem is The correct functional for the boundary value problem (5.13) and (5.14) is given as Making the correct functional stationary, using as the Lagrange multiplier [1318, 36], we get the following iterative formula: where Applying the variational iteration method using He’s polynomials, Comparing the coefficient of like powers of p, consequently, we obtain the following approximants: The series solution is given as Imposing the boundary conditions at and using leads to the following system: The solution of the above system gives
The series solution is given by which is in full agreement with [47].

Table 2 shows the exact values and the errors obtained by using the homotopy perturbation method (HPM) [4], variational iteration method (VIM) [5], decomposition method (ADM) [7], the sixth degree B-spline method [1], iterative method (ITM) [6], and the variational iteration method using He’s polynomials (VIMHP) for The table clearly indicates the improvements as compared with B-spline method. Higher accuracy can be obtained by evaluating more components of .

Table 2: Error estimates.

Remark 5.4. The numerical results clearly indicate that the results obtained by HPM, VIM, ADM, ITM, and the proposed VIMHP are the same. Moreover, it shows the improvements as compare to B-spline method.

#### 6. Conclusions

In this paper, we applied the variational iteration method using He’s polynomials (VIMHP) for finding the solution of boundary value problems of fifth-order. The method is applied in a direct way without using linearization, transformation, discretization, or restrictive assumptions. It may be concluded that VIMHP is very powerful and efficient in finding the analytical solutions for a wide class of boundary value problems. The method gives more realistic series solutions that converge very rapidly in physical problems. It is worth mentioning that the method is capable of reducing the volume of the computational work as compared to the classical methods while still maintaining the high accuracy of the numerical result, the size reduction amounts to the improvement of performance of approach. The fact that the VIMHP solves nonlinear problems without using Adomian’s polynomials is a clear advantage of this technique over the decomposition method.

#### Acknowledgments

The authors are highly grateful to both the Referees and Professor David Chelidze for their constructive comments. They would like to thank Dr. S. M. Junaid Zaidi, Rector CIIT, for providing excellent research environment and facilities.

#### References

1. H. N. Çaglar, S. H. Çaglar, and E. H. Twizell, “The numerical solution of fifth-order boundary value problems with sixth-degree B-spline functions,” Applied Mathematics Letters, vol. 12, no. 5, pp. 25–30, 1999.
2. A. R. Davies, A. Karageorghis, and T. N. Phillips, “Spectral Galerkin methods for the primary two-point boundary value problem in modelling viscoelastic flows,” International Journal for Numerical Methods in Engineering, vol. 26, no. 3, pp. 647–662, 1988.
3. D. J. Fyfe, “Linear dependence relations connecting equal interval Nth degree splines and their derivatives,” Journal of the Institute of Mathematics and Its Applications, vol. 7, no. 3, pp. 398–406, 1971.
4. M. A. Noor and S. T. Mohyud-Din, “An efficient algorithm for solving fifth-order boundary value problems,” Mathematical and Computer Modelling, vol. 45, no. 7-8, pp. 954–964, 2007.
5. M. A. Noor and S. T. Mohyud-Din, “Variational iteration technique for solving higher order boundary value problems,” Applied Mathematics and Computation, vol. 189, no. 2, pp. 1929–1942, 2007.
6. M. A. Noor and S. T. Mohyud-Din, “Modified decomposition method for solving linear and nonlinear fifth-order boundary value problems,” International Journal of Applied Mathematics and Computer Science. In press.
7. A.-M. Wazwaz, “The numerical solution of fifth-order boundary value problems by the decomposition method,” Journal of Computational and Applied Mathematics, vol. 136, no. 1-2, pp. 259–270, 2001.
8. J.-H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 1-2, pp. 87–88, 2006.
9. J.-H. He, “Comparison of homotopy perturbation method and homotopy analysis method,” Applied Mathematics and Computation, vol. 156, no. 2, pp. 527–539, 2004.
10. J.-H. He, “Homotopy perturbation method for bifurcation of nonlinear problems,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 6, no. 2, pp. 207–208, 2005.
11. J.-H. He, “The homotopy perturbation method for nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004.
12. J.-H. He, “A coupling method of a homotopy technique and a perturbation technique for nonlinear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.
13. 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.
14. J.-H. He, “Variational iteration method—a kind of nonlinear analytical technique: some examples,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 699–708, 1999.
15. J.-H. He, “Variational iteration method for autonomous ordinary differential systems,” Applied Mathematics and Computation, vol. 114, no. 2-3, pp. 115–123, 2000.
16. 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.
17. J.-H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.
18. J.-H. He, “The variational iteration method for eighth-order initial-boundary value problems,” Physica Scripta, vol. 76, no. 6, pp. 680–682, 2007.
19. M. Inokuti, H. Sekine, and T. Mura, “General use of the Lagrange multiplier in nonlinear mathematical physics,” in Variational Method in the Mechanics of Solids, S. Nemat-Naseer, Ed., pp. 156–162, Pergamon Press, New York, NY, USA, 1978.
20. S. Abbasbandy, “Numerical solution of nonlinear Klein-Gordon equations by variational iteration method,” International Journal for Numerical Methods in Engineering, vol. 70, no. 7, pp. 876–881, 2007.
21. 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.
22. S. Abbasbandy, “Numerical method for nonlinear wave and diffusion equations by the variational iteration method,” International Journal for Numerical Methods in Engineering, vol. 73, no. 12, pp. 1836–1843, 2008.
23. A. Ghorbani and J. Saberi-Nadjafi, “He's homotopy perturbation method for calculating adomian polynomials,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 229–232, 2007.
24. A. Ghorbani, “Beyond Adomian's polynomials: He polynomials,” Chaos, Solitons & Fractals. In press.
25. 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.
26. S. T. Mohyud-Din, “A reliable algorithm for Blasius equation,” in Proceedings of the International Conference of Mathematical Sciences (ICMS '07), pp. 616–626, Bangi-Putrajaya, Malaysia, November 2007.
27. 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.
28. M. A. Noor and S. T. Mohyud-Din, “Homotopy perturbation method for solving sixth-order boundary value problems,” Computers and Mathematics with Applications. In press.
29. M. A. Noor and S. T. Mohyud-Din, “Homotopy method for solving eighth order boundary value problems,” Journal of Mathematical Analysis and Approximation Theory, vol. 1, no. 2, pp. 161–169, 2006.
30. M. A. Noor and S. T. Mohyud-Din, “Approximate solutions of Flieral Petviashivili equation and its variants,” International Journal of Mathematics and Computer Science, vol. 2, no. 4, pp. 345–360, 2007.
31. M. A. Noor and S. T. Mohyud-Din, “A reliable approach for solving linear and nonlinear sixth-order boundary value problems,” International Journal of Computational and Applied Mathematics. In press.
32. M. A. Noor and S. T. Mohyud-Din, “An efficient method for fourth-order boundary value problems,” Computers and Mathematics with Applications, vol. 54, no. 7-8, pp. 1101–1111, 2007.
33. M. A. Noor and S. T. Mohyud-Din, “Variational iteration technique for solving tenth-order boundary value problems,” Applied Mathematics and Computation, vol. 189, no. 2, pp. 1929–1942, 2007.
34. M. A. Noor and S. T. Mohyud-Din, “Variational iteration decomposition method for solving eighth-order boundary value problems,” Differential Equations and Nonlinear Mechanics. In press.
35. J. I. Ramos, “On the variational iteration method and other iterative techniques for nonlinear differential equations,” Applied Mathematics and Computation. In press.
36. L. Xu, “The variational iteration method for fourth-order boundary value problems,” Chaos, Solitons & Fractals. In press.