- Research Article
- Open access
- Published:
Comparison between the Variational Iteration Method and the Homotopy Perturbation Method for the Sturm-Liouville Differential Equation
Boundary Value Problems volume 2010, Article number: 317369 (2010)
Abstract
We applied the variational iteration method and the homotopy perturbation method to solve Sturm-Liouville eigenvalue and boundary value problems. The main advantage of these methods is the flexibility to give approximate and exact solutions to both linear and nonlinear problems without linearization or discretization. The results show that both methods are simple and effective.
1. Introduction
The variational iteration method (VIM) [1–4] and homotopy perturbation method (HPM) [5–8], proposed by He, are powerful analytical methods for various kinds of linear and nonlinear problems. For example, the variational iteration method has been applied to autonomous ordinary differential equation [9] and delay differential equation [10]. Abdou and Soliman applied this method to Schrodinger-KDV, generalized KDV, and Shallow water equations [11], Burger's equations, and coupled Burger's equations [12]. Furthermore, Momani and Abuasad [13] used VIM for Helmoltz partial equation. Also homotopy perturbation method was successfully applied to Voltra's integrodifferential equation [14], boundary value problem [8], nonlinear wave equations [15], and so forth; see [16–20]. In this paper, we exert these methods for linear Sturm-Liouville eigenvalue and boundary value problems (BVPs). A linear Sturm-Liouville operator has the form
where
and is known analytic function representing the nonhomogeneous term. Associated with the differential equation (1.1) are the following separated homogeneous boundary conditions:
where and are arbitrary constants. For simplicity, we will assume that and are continuous. The values of for which BVP has a nontrivial solution are called eigenvalues of , and a nontrivial solution corresponding to an eigenvalue is called an eigenfunction.
The paper is organized as follows: in Sections 2 and 3, an analysis of the variational iteration and homotopy perturbation methods will be given. In Section 4, we apply HPM to solve Sturm-Liouville problems. We present 3 examples to show the efficiency and simplicity of the proposed methods in Section 5. Finally, we give our conclusions in Section 6.
2. He's Variational Iteration Method
To illustrate the basic concept of He's variational iteration method [1–4], we consider the following nonlinear differential equation:
where is a linear operator, is a nonlinear operator, and is a nonhomogeneous term. He has modified the general Lagrange multiplier method into an iteration method which is called correction functional as follows [1–4, 9]:
where is a general Lagrange multiplier, which can be identified optimally via the variational theory [21]. The subscript denotes the th approximation, and is considered as a restricted variation [1–4], that is, . Employing the restricted variation in (2.2) makes it easy to compute the Lagrange multiplier; see [22, 23]. It is shown that this method is very effective and easy and can solve a large class of nonlinear problems. For linear problems, its exact solution can be obtained only one iteration because can be exactly identified.
3. Homotopy Perturbation Method
In this section, we will present a review of the homotopy perturbation method. To clarify the basic idea of the HPM [5–8], we consider the following nonlinear differential equation:
with boundary conditions
where is a general differential operator, is a boundary operator, is a known analytic function, and is the boundary of the domain . The operator can, generally speaking, be divided into parts and while is nonlinear. Equation (3.1), therefore, can be rewritten as follows:
By the homotopy technique, we construct a homotopy as follows:
which satisfies
or
where is an embedding parameter, and is an initial approximation of (3.1) which satisfies the boundary conditions. Obviously, from (3.5), we have
The changing process of from zero to unity is just that of from to . In topology, this is called deformation and , and are called homotopic. According to HPM, we can assume that the solution of (3.5) can be written as a power series in :
Setting results in the approximate solution  (3.2):
The coupling of the perturbation method and the homotopy method is called the homotopy perturbation method which has eliminated limitations of the traditional perturbation method. On the other hand, the proposed technique can take full advantage of the traditional perturbations techniques.
4. Applying HPM to Solve Sturm-Liouville Problem
To solve (1.1), by means of homotopy perturbation method, we choose linear operator
with the property , where is constant of integration and suggests that we define a nonlinear operator as . Also is known analytic function representing the nonhomogeneous term. Therefore, (1.1) can be rewritten as follows:
By the homotopy perturbation technique proposed by He [5–8], we can construct a homotopy
or
One may now try to obtain a solution of (4.2) in the form
where the for are functions yet to be determined. Substituting (4.5) into (4.4) yields
Collecting terms of the same powers of yields
The initial approximation or can be freely chosen.
5. The Applications
To incorporate our discussion above, three special cases of the Sturm-Liouville equation (1.1) will be studied.
Example 5.1.
Consider the Sturm-Liouville equation
with initial approximation
where and are constants. To solve (5.1) using the VIM, we have correction functional
where is Lagrange multiplier. Making the above correction functional stationary, we can obtain the following stationary conditions:
The Lagrange multiplier can, therefore, be identified as
Substituting (5.5) for correction functional (5.3), we have the following iteration formula:
Using the iteration formula (5.6) and initial approximation (5.2), we get
In the same way, we obtain
which means that
is the exact solution of (5.1).
In order to solve (5.1) using the HPM according to (4.4), we can readily construct a homotopy which satisfies
or
We consider as
Substituting (5.12) into (5.11), collecting terms of the same power, and using initial approximation, we have the following set of linear equations:
Solving the above equations, we have
Continuing in this manner, we can obtain
which is exactly the same as that obtained by VIM.
Example 5.2.
As another example, we consider Sturm-Liouville problem
with initial conditions
where and are constants. To solve (5.16) by means of variational method, we construct a correction functional
where is the Lagrange multiplier and denotes restricted variation that is . Then, we have
Calculus of variations and integration by parts give the stationary conditions
for which the Lagrange multiplier should satisfy. The Lagrange multiplier can, therefore, be identified as
Substituting (5.21) into correction functional (5.18) results in the following iteration formula:
According to initial conditions (5.17), it is natural to choose initial approximation Using the above variational formula (5.22), we can obtain the following result:
In order to solve system (5.16)-(5.17) using HPM, after applying HPM and rearranging based on powers of -terms, we have
Solving the above equations, we get
Example 5.3.
Finally, we consider eigenvalue Sturm-Liouville problem
along with the Dirichlet boundary conditions
To solve (5.26) by means of variational method, we construct a correction functional for (5.26) that reads as
where is Lagrange multiplier. Following the discussion presented in the previous example, we obtain the following iteration formula:
Let us begin with an initial approximation where and are constants to be determined. Substituting the proposed initial iterate in (5.29) gives
In the same way, we obtain
So, we can derive that
is the exact solution of (5.26).
In order to solve (5.26) using HPM, similar to previous examples, after applying HPM and rearranging based on powers of -terms, we have
Now, we choose . Solving the above sets of equations yields
Hence, from (4.4) we get
which is exactly the same as that obtained by VIM. Now, we use the boundary condition (5.27) to obtain eigenvalue and eigenfunctions of (5.26). Imposing the boundary conditions in (5.35) yields
So, there are two infinite sequences of eigenvalues :
Thus, corresponding linearly nontrivial solutions are
Since and are of class , that is, are continuous real-valued functions of , using the definition of inner product on , that is,
and the norm induced by inner product
we get the normalization constants as
Consequently, we obtain
where and are normalized eigenfunctions, that is, and .
6. Conclusion
In this work, we proposed variational method and compared with homotopy perturbation method to solve ordinary Sturm-Liouville differential equation. The variational iteration algorithm used in this paper is the variational iteration algorithm-I; there are also variational iteration algorithm-II and variational iteration algorithm-III [24], which can also be used for the present paper. It may be concluded that the two methods are powerful and efficient techniques to find exact as well as approximate solutions for wide classes of ordinary differential equations.
References
He J-H: A new approach to nonlinear partial differential equations. Communications in Nonlinear Science and Numerical Simulation 1997,2(4):230-235. 10.1016/S1007-5704(97)90007-1
He J-H: Approximate analytical solution for seepage flow with fractional derivatives in porous media. Computer Methods in Applied Mechanics and Engineering 1998,167(1-2):57-68. 10.1016/S0045-7825(98)00108-X
He J-H: Approximate solution of nonlinear differential equations with convolution product nonlinearities. Computer Methods in Applied Mechanics and Engineering 1998,167(1-2):69-73. 10.1016/S0045-7825(98)00109-1
He J-H: Variational iteration method—a kind of non-linear analytical technique: some examples. International Journal of Non-Linear Mechanics 1999,34(4):699-708. 10.1016/S0020-7462(98)00048-1
He J-H: Homotopy perturbation technique. Computer Methods in Applied Mechanics and Engineering 1999,178(3-4):257-262. 10.1016/S0045-7825(99)00018-3
He J-H: Homotopy perturbation method: a new nonlinear analytical technique. Applied Mathematics and Computation 2003,135(1):73-79. 10.1016/S0096-3003(01)00312-5
He J-H: A coupling method of a homotopy technique and a perturbation technique for non-linear problems. International Journal of Non-Linear Mechanics 2000,35(1):37-43. 10.1016/S0020-7462(98)00085-7
He J-H: Homotopy perturbation method for solving boundary value problems. Physics Letters A 2006,350(1-2):87-88. 10.1016/j.physleta.2005.10.005
He J-H: Variational iteration method for autonomous ordinary differential systems. Applied Mathematics and Computation 2000,114(2-3):115-123. 10.1016/S0096-3003(99)00104-6
He J-H: Variational iteration method for delay differential equations. Communications in Nonlinear Science and Numerical Simulation 1997,2(4):235-236. 10.1016/S1007-5704(97)90008-3
Abdou MA, Soliman AA: New applications of variational iteration method. Physica D 2005,211(1-2):1-8. 10.1016/j.physd.2005.08.002
Abdou MA, Soliman AA: Variational iteration method for solving Burger's and coupled Burger's equations. Journal of Computational and Applied Mathematics 2005,181(2):245-251. 10.1016/j.cam.2004.11.032
Momani S, Abuasad S: Application of He's variational iteration method to Helmholtz equation. Chaos, Solitons & Fractals 2006,27(5):1119-1123. 10.1016/j.chaos.2005.04.113
El-Shahed M: Application of He's homotopy perturbation method to Volterra's integro-differential equation. International Journal of Nonlinear Sciences and Numerical Simulation 2005,6(2):163-168. 10.1515/IJNSNS.2005.6.2.163
He J-H: Application of He's homotopy perturbation method to nonlinear wave equations. Chaos Solitons & Fractals 1998,167(1-2):69-73.
Abbasbandy S: Homotopy perturbation method for quadratic Riccati differential equation and comparison with Adomian's decomposition method. Applied Mathematics and Computation 2006,172(1):485-490. 10.1016/j.amc.2005.02.014
Ganji DD, Sadighi A: Application of He's homotopy-perturbation method to nonlinear coupled systems of reaction-diffusion equations. International Journal of Nonlinear Sciences and Numerical Simulation 2006,7(4):411-418. 10.1515/IJNSNS.2006.7.4.411
He J-H: Homotopy perturbation method for bifurcation of nonlinear problems. International Journal of Nonlinear Sciences and Numerical Simulation 2005,6(2):207-208. 10.1515/IJNSNS.2005.6.2.207
He J-H: The homotopy perturbation method nonlinear oscillators with discontinuities. Applied Mathematics and Computation 2004,151(1):287-292. 10.1016/S0096-3003(03)00341-2
Rafei M, Ganji DD: Explicit solutions of Helmholtz equation and fifth-order KdV equation using homotopy perturbation method. International Journal of Nonlinear Sciences and Numerical Simulation 2006,7(3):321-328. 10.1515/IJNSNS.2006.7.3.321
He J-H: Variational principles for some nonlinear partial differential equations with variable coefficients. Chaos, Solitons & Fractals 2004,19(4):847-851. 10.1016/S0960-0779(03)00265-0
He J-H: Semi-inverse method of establishing generalized variational principles for fluid mechanics with emphasis on turbomachinery aerodynamics. International Journal of Turbo and Jet Engines 1997,14(1):23-28. 10.1515/TJJ.1997.14.1.23
He J-H, Wan Y-Q, Guo Q: An iteration formulation for normalized diode characteristics. International Journal of Circuit Theory and Applications 2004,32(6):629-632. 10.1002/cta.300
He J-H, Wu G-C, Austin F: The variational iteration method which should be followed. Nonlinear Science Letters A 2010,1(1):1-30.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Neamaty, A., Darzi, R. Comparison between the Variational Iteration Method and the Homotopy Perturbation Method for the Sturm-Liouville Differential Equation. Bound Value Probl 2010, 317369 (2010). https://doi.org/10.1155/2010/317369
Received:
Accepted:
Published:
DOI: https://doi.org/10.1155/2010/317369