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

Recommended by David Chelidze 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.


Introduction
In this paper, we general fifth-order boundary value problem y v x g x y q x 1.1 with boundary conditions where y x and f x, y are real and as many times differentiable as required for x ∈ a, b ; and A i , i 1, 2, 3 and B i , i 1, 2 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.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  .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][14][15][16][17][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 p 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.

Variational iteration method
To illustrate the basic concept of the technique, we consider the following general differential equation: where L is a linear operator, N a nonlinear operator, and g x 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 nth approximation, u n is considered as a restricted variation, that is, δ u n 0; 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 u n 1 , n ≥ 0 of the solution u will be readily obtained upon using the determined Lagrange multiplier and any selective function u 0 , consequently, the solution is given by For the convergence and error estimates of variational iteration method, see Ramos 35 .

Homotopy perturbation method
To explain the homotopy perturbation method, we consider a general equation of the type where L is an integral or differential operator.We define a convex homotopy H u, p by where F u is a functional operator with known solutions v 0 , which can be obtained easily.It is clear that for This shows that H u, p continuously traces an implicitly defined curve from a starting point H v 0 , 0 to a solution function H f, 1 .The embedding parameter monotonically increases from zero to unit as the trivial problem F u 0 is continuously deforms the original problem L u 0. 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 p → 1, 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 L u ; 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.

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 L is a linear operator, N a nonlinear operator, and g x 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 nth approximation, u n is considered as a restricted variation, that is, δ u n 0; 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.

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 .

5.2
The exact solution for this problem is y x e x .

5.3
The correct functional for the boundary value problem 5.1 and 5.2 is given as λ s d 5 y n dx 5 − e −x y 2 n x ds.

5.5
where A y 0 , B y iv 0 .

5.6
Applying the variational iteration method using He's polynomials,

5.7
Comparing the coefficient of like powers of p, consequently, we obtain the following approximants:

5.11
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 x 0.0, 0.1, 0.2, . . ., 1.0.The table clearly indicates the improvements as compared with B-spline method.Higher accuracy can be obtained by evaluating more components of y x .
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.

5.14
The exact solution of the problem is

5.15
The correct functional for the boundary value problem 5.13 and 5.14 is given as Making the correct functional stationary, using λ 1/4! s − x 4 , as the Lagrange multiplier 13-18, 36 , we get the following iterative formula:

5.17
where A y 0 , B y iv 0 .

5.18
Applying the variational iteration method using He's polynomials,

5.19
Comparing the coefficient of like powers of p, consequently, we obtain the following approximants:  The series solution is given as x 14 O x 15 .

5.21
Imposing the boundary conditions at x 1 and using y 1 y 1 −e leads to the following system:

5.22
The solution of the above system gives A −2.99999988, B −8.00000054.

5.24
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 x 0.0, 0.1, 0.2, . . ., 1.0.The table clearly indicates the improvements as compared with B-spline method.Higher accuracy can be obtained by evaluating more components of y x .
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.

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.

Example 5 .1 see 4 - 7 .
Consider the following nonlinear boundary value problem of fifthorder,
* Error exact solution − series solution.
* Error exact solution − series solution.