SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Degenerate and Singular Differential Operators with Applications to Boundary Value Problems.

Open Access Research Article

Comparison between the Variational Iteration Method and the Homotopy Perturbation Method for the Sturm-Liouville Differential Equation

A Neamaty1 and R Darzi2*

Author Affiliations

1 Department of Mathematics, University of Mazandaran, Babolsar 47416-95447, Iran

2 Department of Mathematics, Islamic Azad University Neka Branch, Neka 48411-86114, Iran

For all author emails, please log on.

Boundary Value Problems 2010, 2010:317369  doi:10.1155/2010/317369


The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2010/1/317369


Received:28 October 2009
Accepted:10 April 2010
Published:17 May 2010

© 2010 The Author(s)

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.

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) [14] and homotopy perturbation method (HPM) [58], 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 [1620]. 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

(11)

where

(12)

and is known analytic function representing the nonhomogeneous term. Associated with the differential equation (1.1) are the following separated homogeneous boundary conditions:

(13)

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 [14], we consider the following nonlinear differential equation:

(21)

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 [14, 9]:

(22)

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 [14], 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 [58], we consider the following nonlinear differential equation:

(31)

with boundary conditions

(32)

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:

(33)

By the homotopy technique, we construct a homotopy as follows:

(34)

which satisfies

(35)

or

(36)

where is an embedding parameter, and is an initial approximation of (3.1) which satisfies the boundary conditions. Obviously, from (3.5), we have

(37)

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 :

(38)

Setting results in the approximate solution  (3.2):

(39)

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

(41)

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:

(42)

By the homotopy perturbation technique proposed by He [58], we can construct a homotopy

(43)

or

(44)

One may now try to obtain a solution of (4.2) in the form

(45)

where the for are functions yet to be determined. Substituting (4.5) into (4.4) yields

(46)

Collecting terms of the same powers of yields

(47)

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

(51)

with initial approximation

(52)

where and are constants. To solve (5.1) using the VIM, we have correction functional

(53)

where is Lagrange multiplier. Making the above correction functional stationary, we can obtain the following stationary conditions:

(54)

The Lagrange multiplier can, therefore, be identified as

(55)

Substituting (5.5) for correction functional (5.3), we have the following iteration formula:

(56)

Using the iteration formula (5.6) and initial approximation (5.2), we get

(57)

In the same way, we obtain

(58)

which means that

(59)

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

(510)

or

(511)

We consider as

(512)

Substituting (5.12) into (5.11), collecting terms of the same power, and using initial approximation, we have the following set of linear equations:

(513)

Solving the above equations, we have

(514)

Continuing in this manner, we can obtain

(515)

which is exactly the same as that obtained by VIM.

Example 5.2.

As another example, we consider Sturm-Liouville problem

(516)

with initial conditions

(517)

where and are constants. To solve (5.16) by means of variational method, we construct a correction functional

(518)

where is the Lagrange multiplier and denotes restricted variation that is . Then, we have

(519)

Calculus of variations and integration by parts give the stationary conditions

(520)

for which the Lagrange multiplier should satisfy. The Lagrange multiplier can, therefore, be identified as

(521)

Substituting (5.21) into correction functional (5.18) results in the following iteration formula:

(522)

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:

(523)

In order to solve system (5.16)-(5.17) using HPM, after applying HPM and rearranging based on powers of -terms, we have

(524)

Solving the above equations, we get

(525)

Example 5.3.

Finally, we consider eigenvalue Sturm-Liouville problem

(526)

along with the Dirichlet boundary conditions

(527)

To solve (5.26) by means of variational method, we construct a correction functional for (5.26) that reads as

(528)

where is Lagrange multiplier. Following the discussion presented in the previous example, we obtain the following iteration formula:

(529)

Let us begin with an initial approximation where and are constants to be determined. Substituting the proposed initial iterate in (5.29) gives

(530)

In the same way, we obtain

(531)

So, we can derive that

(532)

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

(533)

Now, we choose . Solving the above sets of equations yields

(534)

Hence, from (4.4) we get

(535)

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

(536)

So, there are two infinite sequences of eigenvalues :

(537)

Thus, corresponding linearly nontrivial solutions are

(538)

Since and are of class , that is, are continuous real-valued functions of , using the definition of inner product on , that is,

(539)

and the norm induced by inner product

(540)

we get the normalization constants as

(541)

Consequently, we obtain

(542)

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

  1. He, J-H: A new approach to nonlinear partial differential equations. Communications in Nonlinear Science and Numerical Simulation. 2(4), 230–235 (1997). Publisher Full Text OpenURL

  2. He, J-H: Approximate analytical solution for seepage flow with fractional derivatives in porous media. Computer Methods in Applied Mechanics and Engineering. 167(1-2), 57–68 (1998). Publisher Full Text OpenURL

  3. He, J-H: Approximate solution of nonlinear differential equations with convolution product nonlinearities. Computer Methods in Applied Mechanics and Engineering. 167(1-2), 69–73 (1998). Publisher Full Text OpenURL

  4. He, J-H: Variational iteration method—a kind of non-linear analytical technique: some examples. International Journal of Non-Linear Mechanics. 34(4), 699–708 (1999). Publisher Full Text OpenURL

  5. He, J-H: Homotopy perturbation technique. Computer Methods in Applied Mechanics and Engineering. 178(3-4), 257–262 (1999). Publisher Full Text OpenURL

  6. He, J-H: Homotopy perturbation method: a new nonlinear analytical technique. Applied Mathematics and Computation. 135(1), 73–79 (2003). Publisher Full Text OpenURL

  7. He, J-H: A coupling method of a homotopy technique and a perturbation technique for non-linear problems. International Journal of Non-Linear Mechanics. 35(1), 37–43 (2000). Publisher Full Text OpenURL

  8. He, J-H: Homotopy perturbation method for solving boundary value problems. Physics Letters A. 350(1-2), 87–88 (2006). Publisher Full Text OpenURL

  9. He, J-H: Variational iteration method for autonomous ordinary differential systems. Applied Mathematics and Computation. 114(2-3), 115–123 (2000). Publisher Full Text OpenURL

  10. He, J-H: Variational iteration method for delay differential equations. Communications in Nonlinear Science and Numerical Simulation. 2(4), 235–236 (1997). Publisher Full Text OpenURL

  11. Abdou, MA, Soliman, AA: New applications of variational iteration method. Physica D. 211(1-2), 1–8 (2005). Publisher Full Text OpenURL

  12. Abdou, MA, Soliman, AA: Variational iteration method for solving Burger's and coupled Burger's equations. Journal of Computational and Applied Mathematics. 181(2), 245–251 (2005). Publisher Full Text OpenURL

  13. Momani, S, Abuasad, S: Application of He's variational iteration method to Helmholtz equation. Chaos, Solitons & Fractals. 27(5), 1119–1123 (2006). PubMed Abstract | Publisher Full Text OpenURL

  14. El-Shahed, M: Application of He's homotopy perturbation method to Volterra's integro-differential equation. International Journal of Nonlinear Sciences and Numerical Simulation. 6(2), 163–168 (2005). Publisher Full Text OpenURL

  15. He, J-H: Application of He's homotopy perturbation method to nonlinear wave equations. Chaos Solitons & Fractals. 167(1-2), 69–73 (1998). PubMed Abstract OpenURL

  16. Abbasbandy, S: Homotopy perturbation method for quadratic Riccati differential equation and comparison with Adomian's decomposition method. Applied Mathematics and Computation. 172(1), 485–490 (2006). Publisher Full Text OpenURL

  17. 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. 7(4), 411–418 (2006). Publisher Full Text OpenURL

  18. He, J-H: Homotopy perturbation method for bifurcation of nonlinear problems. International Journal of Nonlinear Sciences and Numerical Simulation. 6(2), 207–208 (2005). Publisher Full Text OpenURL

  19. He, J-H: The homotopy perturbation method nonlinear oscillators with discontinuities. Applied Mathematics and Computation. 151(1), 287–292 (2004). Publisher Full Text OpenURL

  20. 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. 7(3), 321–328 (2006). Publisher Full Text OpenURL

  21. He, J-H: Variational principles for some nonlinear partial differential equations with variable coefficients. Chaos, Solitons & Fractals. 19(4), 847–851 (2004). PubMed Abstract | Publisher Full Text OpenURL

  22. 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. 14(1), 23–28 (1997). Publisher Full Text OpenURL

  23. He, J-H, Wan, Y-Q, Guo, Q: An iteration formulation for normalized diode characteristics. International Journal of Circuit Theory and Applications. 32(6), 629–632 (2004). Publisher Full Text OpenURL

  24. He, J-H, Wu, G-C, Austin, F: The variational iteration method which should be followed. Nonlinear Science Letters A. 1(1), 1–30 (2010)