Mathematical Problems in Engineering

Volume 2008, Article ID 954794, 11 pages

http://dx.doi.org/10.1155/2008/954794

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

Academic Editor: David Chelidze

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 [1–7] 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 [8–18] 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 [13–18]. 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 [13–18]. The homotopy perturbation method [8–12, 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, 8–36] and the references therein. In these methods, the solution is given in an infinite series usually converging to an accurate solution, see [4, 5, 8–22, 25–36] 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 [4–6] 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 [13–22] 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, 13–22, 25, 32–36], we
can construct a correct functional as follows: where is a Lagrange multiplier [13–19], which can be
identified optimally via variational iteration method. The subscripts *n* denote the *n*th 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 [13–19]. 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, 8–11, 23, 24]. The homotopy perturbation method uses the
homotopy parameter *p* as an expanding
parameter [4, 8–12, 23, 24, 27–31] 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, 8–12, 23, 24, 27–31]. 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, 13–22, 25, 26, 32–36], we can
construct a correct functional as follows: where is a Lagrange multiplier [13–19], which can
be identified optimally via variational iteration method. The subscripts *n* denote the *n*th 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 [4–7].

*Example 5.1 (See [4–7]). * 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 [13–18, 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 [4–7].

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 .

*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 [4–7]). * 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 [13–18, 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 [4–7].

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 .

*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

- 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - D. J. Fyfe, “Linear dependence relations connecting equal interval
*N*th degree splines and their derivatives,”*Journal of the Institute of Mathematics and Its Applications*, vol. 7, no. 3, pp. 398–406, 1971. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He, “Homotopy perturbation method for solving boundary value problems,”
*Physics Letters A*, vol. 350, no. 1-2, pp. 87–88, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - J.-H. He, “Comparison of homotopy perturbation method and homotopy analysis method,”
*Applied Mathematics and Computation*, vol. 156, no. 2, pp. 527–539, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Google Scholar - J.-H. He, “The homotopy perturbation method for nonlinear oscillators with discontinuities,”
*Applied Mathematics and Computation*, vol. 151, no. 1, pp. 287–292, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar - J.-H. He, “Variational iteration method for autonomous ordinary differential systems,”
*Applied Mathematics and Computation*, vol. 114, no. 2-3, pp. 115–123, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He, “The variational iteration method for eighth-order initial-boundary value problems,”
*Physica Scripta*, vol. 76, no. 6, pp. 680–682, 2007. View at Publisher · View at Google Scholar - 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. View at Google Scholar - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Publisher · View at Google Scholar - 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. View at Google Scholar - A. Ghorbani, “Beyond Adomian's polynomials: He polynomials,”
*Chaos, Solitons & Fractals*. In press. View at Publisher · 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. - 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. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar - 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. View at Google Scholar · View at MathSciNet - 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. View at Google Scholar - 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. - 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. View at Publisher · View at Google Scholar - 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. View at Publisher · View at Google Scholar - 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. - J. I. Ramos, “On the variational iteration method and other iterative techniques for nonlinear differential equations,”
*Applied Mathematics and Computation*. In press. - L. Xu, “The variational iteration method for fourth-order boundary value problems,”
*Chaos, Solitons & Fractals*. In press.