`Journal of Applied MathematicsVolume 2012, Article ID 163821, 16 pageshttp://dx.doi.org/10.1155/2012/163821`
Research Article

## A Coupled Method of Laplace Transform and Legendre Wavelets for Lane-Emden-Type Differential Equations

1College of Computer, National University of Defense Technology, Changsha 410073, China
2China Aerodynamics Research and Development Center, Sichuan, Mianyang 621000, China

Received 4 September 2012; Revised 11 October 2012; Accepted 15 October 2012

Copyright © 2012 Fukang Yin et al. 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

A coupled method of Laplace transform and Legendre wavelets is presented to obtain exact solutions of Lane-Emden-type equations. By employing properties of Laplace transform, a new operator is first introduced and then its Legendre wavelets operational matrix is derived to convert the Lane-Emden equations into a system of algebraic equations. Block pulse functions are used to calculate the Legendre wavelets coefficient matrices of the nonlinear terms. The results show that the proposed method is very effective and easy to implement.

#### 1. Introduction

The Lane-Emden-type equation has attracted much attention from mathematicians and physicists, since it is widely used to investigate the theory of stellar structure, the thermal behavior of spherical cloud gas, and theory of thermionic currents [13]. One of the general forms of Lane-Emden type equations is subject to conditions where and are constants and is a nonlinear function of and .

The solution of the Lane-Emden equations is numerically challenging due to the singularity behavior at the origin and nonlinearities. Therefore, much attention has been paid to searching for the better and more efficient methods for determining a solution, approximate or exact, analytical or numerical, to the Lane-Emden equations. The existing methods fall into two groups: the analytical methods and the numerical ones. The analytical methods express the exact solution of the equation in the form of elementary functions and convergent function series, such as the Adomian decomposition method (ADM) [47], homotopy perturbation method (HPM) [810], variational iteration method (VIM) [1113], homotopy analysis method (HAM) [14, 15], power series solution [16], and differential transformation method [17, 18]. Unlike the analytical methods, the numerical ones approximate the exact solution on a finite set of distinct points, such as the Legendre Tau method [19] and sinc-collocation method [20], Lagrangian approach [21], successive linearization method [22], wavelets and collocation method [2326], and spectral method [27].

The Laplace transform is a wonderful tool for solving linear differential equations and has enjoyed much success in this realm. However, it is totally incapable of handling nonlinear equations because of the difficulties caused by nonlinear terms. Since Laplace Adomian decomposition method (LADM) was proposed by Khuri [28] and then developed by Khan [29] and Khan and Gondal [30], the couple methods that based on Laplace transform and other methods have received considerable attention in the literature. What is more, the homotopy perturbation method [31] and the variational iteration method [32] are combined with the well-known Laplace transform to develop a highly effective technique for handling many nonlinear problems. For example, the coupled methods [3338] based on the homotopy perturbation method and Laplace transform have been proved to be very effective for the solution of nonlinear problems.

Wavelets theory, as a relatively new and emerging area in mathematical research, has received considerable attention in dealing with various problems of dynamic systems. The fundamental idea of this technique is that it reduces these problems to those of solving a system of algebraic equations, thus greatly simplifies the problem and reduces the computation cost [6]. Moreover, wavelets establish a connection with fast numerical algorithms. Yousefi [25] has obtained the numerical solutions of the Lane-Emden equations (1.1) by converting it into an integral equation and then using Legendre wavelets and Gaussian integration method. Furthermore, Pandey et al. [26] obtained the numerical solutions of Lane-Emden equations by using the operational matrix of derivative of the Legendre polynomials.

Motivated and inspired by the ongoing research in these areas, we propose a coupled method of Laplace transform and Legendre wavelets to establish exact solutions of (1.1). The advantage of our method is its capability of combining the two powerful methods to obtain exact solutions of nonlinear equations. The remainder of the paper is organized as follows. In Section 2, we describe some preliminaries about Legendre wavelets and Block function pulses. The proposed method is detailed in Section 3. Four examples are given in Section 4 to demonstrate the validity and applicability of the proposed method. Finally the concluding remarks are given in Section 5.

#### 2. Preliminaries

##### 2.1. Legendre Wavelets

Legendre wavelets have four arguments: , , is any positive integer, is the order for Legendre polynomials, and is the normalized time. They are defined on the interval as follows: where . The coefficient is for orthonormality, the dilation parameter , and the translation parameter . Here, are Legendre polynomials of order defined on the interval .

A function defined over , may be expanded by Legendre wavelet series as with and denoting the inner product.

If the infinite series in (2.2) is truncated, then it can be written as where and are matrices defined as

##### 2.2. Block Pulse Functions

Block pulse functions (BPFs) form a complete set of orthogonal functions that are defined on the interval as for . It is also known that for arbitrary absolutely integrable function on can be expanded in block pulse functions: in which where are the coefficients of the block-pulse function, given by The elementary properties of BPFs are as follows.(1)Disjointness: the BPFs are disjoined with each other in the interval : for .(2)Orthogonality: the BPFs are orthogonal with each other in the interval : for .(3)Completeness: the BPFs set is complete when approaches infinity. This means that for every , when approaches to the infinity, Parseval’s identity holds: where

Definition 2.1. Let and are two vectors of ; then .

Lemma 2.2. Assuming and are two absolutely integrable functions, which can be expanded in block pulse function as and respectively, then one has where .

Proof. By using the property of BPFs in (2.11), we have and this completes the proof.

##### 2.3. Nonlinear Term Approximately

The Legendre wavelets can be expanded into -set of block-pulse functions as

Taking the collocation points as in the following,

The -square Legendre matrix is defined as

The operational matrix of product of Legendre wavelets can be obtained by using the properties of BPFs. Let and   be two absolutely integrable functions, which can be expanded in Legendre wavelets as and , respectively.

From (2.17), we have

By employing Lemma 2.2 and (2.17), we get

#### 3. Laplace Legendre Wavelets Method (LLWM)

In this section, we will briefly demonstrate the utilization of the LLWM for solving the Lane-Emden equations given in (1.1).

By multiplying and then applying the Laplace transform to both sides of (1.1), we obtain where is operator of Laplace transform and .

Equation (3.1) can be rewritten as

By integrating both sides of (3.2) from 0 to with respect to , we have

Taking the inverse Laplace transform to (3.3), we get

By using the initial conditions from (1.2), we have

For convenience, we define an operator . Therefore, (3.5) can be represented as

Now, we will show how to derive the Legendre wavelets operational matrix of operator . First of all, three corollaries are given.

Corollary 3.1. Let be the one-dimensional Legendre wavelets vector defined in (2.6), then one has where is matrix

Proof. Let ; then by expanding by , we obtain where is matrix. For example, when , we have
From (3.9), we obtain and this completes the proof.

Corollary 3.2. Let be the one-dimensional Legendre wavelets vector defined in (2.6), then one has

Proof. From the definition of operator , we can know
By further analysis, we obtain where matrix is give in (3.8) and matrix is defined as
So we finally have Let , we get (3.12), and this completes the proof.

Corollary 3.3. Let be the one dimension Legendre wavelets vector defined in (2.6); then one has where :

In order to use Legendre wavelets method, we approximate ,, and as

Substituting (3.19) into (3.6) and then using the Corollary 3.3, we have

Finally, we can get Equation (3.21) is a nonlinear equation which can be solved for the elements of in (3.19) by the Newton iterative method.

#### 4. Numerical Examples

In this section, four different examples are examined to demonstrate the effectiveness and high accuracy of the LLWM.

Example 4.1. Consider the lane-Emden equation given in [6] subject to the initial conditions with the exact solution .
Applying the method proposed in Section 3, we have

When and , the numerical solution and the absolute error of Example 4.1 are plotted in Figure 1. We see that a good approximation is obtained by using a few terms of Legendre wavelets.

Figure 1: Numerical solution and absolute error of Example 4.1.

Example 4.2. Consider the lane-Emden equation of index : subject to the initial conditions When   and  , the exact solutions are and , respectively.
Applying the method developed in Section 3, we have
From [18], we can get the four terms approximate solution of (4.4) as follows:
In the case of , when and , the results of Example 4.2 are plotted in Figure 2. We observe that the accuracy of LLWM solution is very high. For the case of , which has strong nonlinearity, we plot the results when and in Figure 3. It can be noted that the LLWM solution is very close to the exact solution. In addition, we give the numerical solutions of difference in Figure 4.

Figure 2: Numerical solution and absolute error of Example 4.2 when .
Figure 3: Numerical solution and absolute error of Example 4.2 when .
Figure 4: Numerical solutions of Example 4.2 for difference .

Example 4.3. Consider the Lane-Emden equation given in [7] by subject to the initial conditions with the exact solution .
Applying the method developed in Section 3, we have
When and , the numerical solution and the absolute error of Example 4.3 are plotted in Figure 5. We observe that the accuracy of LLWM is very high and only needs a few terms of Legendre wavelets.

Figure 5: Numerical solution and absolute error of Example 4.3.

Example 4.4. Consider the following nonlinear Lane-Emden differential equation [7, 8, 13]: subject to initial conditions The exact solution is .
Here, we first expand by Taylor series and then have
Considering only the first five terms we can write
We let and then approximate and by the Legendre wavelets as
By applying the LLWM, we have
The analytical methods, such as the ADM [7], HPM [8], and VIM [13], can get the exact solution of Example 4.4. From [7, 8, 13], we can know that the four-term solution of the above methods is .

When and , the numerical solutions of Example 4.4 are plotted in Figure 6. We observe that the LLWM solution is more accurate than the four-term solution of HPM or ADM or VIM. The LLWM can only get the approximate solution because of the error caused by expanding the nonlinear term. However, the LLWM can avoid the symbolic computation and only needs a few terms of Legendre wavelets. Since the exact solution is an analytic function, the higher accuracy can be obtained by taking .

Figure 6: Numerical solution of Example 4.4.

#### 5. Conclusion

In this paper, we have successfully developed a coupled method of Laplace transform and Legendre wavelets (LLWM) for solution of Lane-Emden equations. The advantage of our method is that only small size operational matrix is required to provide the solution at high accuracy. It can be clearly seen in the paper that the proposed method works well even in the case of high nonlinearity. Compared to the analytical methods, such as the ADM, VIM, and HPM, the LLWM can only get the approximate solution because of the error caused by expanding highly nonlinear terms. However, the developed vector-matrix form makes LLWM a promising tool for Lane-Emden-type equations, because LLWM is computer-oriented and can use many existing fast algorithms to reduce the computational cost.

#### Acknowledgments

This work is supported by National Natural Science Foundation of China (Grant no. 41105063). The authors are very grateful to reviewers for carefully reading the paper and for his (her) comments and suggestions which have improved the paper.

#### References

1. S. Chandrasekhar, An Introduction to the Study of Stellar Structure, Dover, New York, NY, USA, 1967.
2. H. T. Davis, Introduction to Nonlinear Differential and Integral Equations, Dover, New York, NY, USA, 1962.
3. O. U. Richardson, The Emission of Electricity From Hot Bodies, Zongmans Green and Company, London, UK, 1921.
4. N. T. Shawagfeh, “Nonperturbative approximate solution for Lane-Emden equation,” Journal of Mathematical Physics, vol. 34, no. 9, pp. 4364–4369, 1993.
5. A.-M. Wazwaz, “A new algorithm for solving differential equations of Lane-Emden type,” Applied Mathematics and Computation, vol. 118, no. 2-3, pp. 287–310, 2001.
6. A.-M. Wazwaz, “A new method for solving singular initial value problems in the second-order ordinary differential equations,” Applied Mathematics and Computation, vol. 128, no. 1, pp. 45–57, 2002.
7. F. Olga and Š. Zdeněk, “A domian decomposition method for certain singular initial value problems,” Journal of Applied Mathematics, vol. 3, no. 2, pp. 91–98, 2010.
8. J. I. Ramos, “Series approach to the Lane-Emden equation and comparison with the homotopy perturbation method,” Chaos, Solitons and Fractals, vol. 38, no. 2, pp. 400–408, 2008.
9. A. Yildirim and T. Öziş, “Solutions of singular IVP's of Lane-Emden type by homotopy pertutbation method,” Physics Letters A, vol. 369, pp. 70–76, 2007.
10. M. S. H. Chowdhury and I. Hashim, “Solutions of a class of singular second-order IVPs by homotopy-perturbation method,” Physics Letters A, vol. 365, no. 5-6, pp. 439–447, 2007.
11. J.-H. He, “Variational approach to the Lane-Emden equation,” Applied Mathematics and Computation, vol. 143, no. 2-3, pp. 539–541, 2003.
12. M. Dehghan and F. Shakeri, “Approximate solution of a differential equation arising in astrophysics using the variational iteration method,” New Astronomy, vol. 13, no. 1, pp. 53–59, 2008.
13. A. Yildirım and T. Öziş, “Solutions of singular IVPs of Lane-Emden type by the variational iteration method,” Nonlinear Analysis, vol. 70, no. 6, pp. 2480–2484, 2009.
14. S. Liao, “A new analytic algorithm of Lane-Emden type equations,” Applied Mathematics and Computation, vol. 142, no. 1, pp. 1–16, 2003.
15. A. S. Bataineh, M. S. M. Noorani, and I. Hashim, “Homotopy analysis method for singular IVPs of Emden-Fowler type,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1121–1131, 2009.
16. C. Mohan and A. R. Al-Bayaty, “Power-series solutions of the Lane-Emden equation,” Astrophysics and Space Science, vol. 73, no. 1, pp. 227–239, 1980.
17. V. S. Ertürk, “Differential transformation method for solving differential equations of Lane-Emden type,” Mathematical & Computational Applications, vol. 12, no. 3, pp. 135–139, 2007.
18. S. Mukherjee, B. Roy, and P. K. Chaterjee, “Solution of Lane-Emden equation by differential transform method,” International Journal of Nonlinear Science, vol. 12, no. 4, pp. 478–484, 2011.
19. K. Parand and M. Razzaghi, “Rational legendre approximation for solving some physical problems on semi-infinite intervals,” Physica Scripta, vol. 69, no. 5, pp. 353–357, 2004.
20. K. Parand and A. Pirkhedri, “Sinc-Collocation method for solving astrophysics equations,” New Astronomy, vol. 15, no. 6, pp. 533–537, 2010.
21. K. Parand, A. R. Rezaei, and A. Taghavi, “Lagrangian method for solving Lane-Emden type equation arising in astrophysics on semi-infinite domains,” Acta Astronautica, vol. 67, no. 7-8, pp. 673–680, 2010.
22. S. S. Motsa and P. Sibanda, “A new algorithm for solving singular IVPs of Lane-Emden type,” in Proceedings of the 4th International Conference on Applied Mathematics, Simulation, Modelling (ASM '10), pp. 176–180, NAUN International Conferences, Corfu Island, Greece, July 2010.
23. S. Karimi Vanani and A. Aminataei, “On the numerical solution of differential equations of Lane-Emden type,” Computers & Mathematics with Applications, vol. 59, no. 8, pp. 2815–2820, 2010.
24. C. Yang and J. Hou, “A numerical method for Lane-Emden equations using hybrid functions and the collocation method,” Journal of Applied Mathematics, vol. 2012, Article ID 316534, 9 pages, 2012.
25. S. A. Yousefi, “Legendre wavelets method for solving differential equations of Lane-Emden type,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1417–1422, 2006.
26. R. K. Pandey, N. Kumar, A. Bhardwaj, and G. Dutta, “Solution of Lane-Emden type equations using legendre operational matrix of differentiation,” Applied Mathematics and Computation, vol. 218, no. 14, pp. 7629–7637, 2012.
27. J. P. Boyd, “Chebyshev spectral methods and the Lane-Emden problem,” Numerical Mathematics, vol. 4, no. 2, pp. 142–157, 2011.
28. S. A. Khuri, “A Laplace decomposition algorithm applied to a class of nonlinear differential equations,” Journal of Applied Mathematics, vol. 1, no. 4, pp. 141–155, 2001.
29. Y. Khan, “An effective modification of the laplace decomposition method for nonlinear equations,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 10, no. 11-12, pp. 1373–1376, 2009.
30. M. Khan and M. A. Gondal, “A new analytical solution of foam drainage equation by laplace decomposition method,” Journal of Advanced Research in Scientific Computing, vol. 2, no. 2, pp. 53–64, 2010.
31. S. Abbasbandy, “Application of He's homotopy perturbation method for Laplace transform,” Chaos, Solitons and Fractals, vol. 30, no. 5, pp. 1206–1212, 2006.
32. S. A. Khuri and A. Sayfy, “A laplace variational iteration strategy for the solution of differential equations,” Applied Mathematics Letters, vol. 25, no. 12, pp. 2298–2305, 2012.
33. M. Madani and M. Fathizadeh, “Homotopy perturbation algorithm using laplace transformation,” Nonlinear Science Letters A, vol. 1, pp. 263–267, 2010.
34. M. Khan, M. A. Gondal, and S. Karimi Vanani, “On the coupling of homotopy perturbation and laplace transformation for system of partial differential equations,” Applied Mathematical Sciences, vol. 6, no. 9–12, pp. 467–478, 2012.
35. M. Madani, M. Fathizadeh, Y. Khan, and A. Yildirim, “On the coupling of the homotopy perturbation method and laplace transformation,” Mathematical and Computer Modelling, vol. 53, no. 9-10, pp. 1937–1945, 2011.
36. Y. Khan and Q. Wu, “Homotopy perturbation transform method for nonlinear equations using He's polynomials,” Computers & Mathematics with Applications, vol. 61, no. 8, pp. 1963–1967, 2011.
37. V. G. Gupta and S. Gupta, “Homotopy perturbation transform method for solving initial boundary value problems of variable coefficients,” International Journal of Nonlinear Science, vol. 12, no. 3, pp. 270–277, 2011.
38. J. Singh, D. Kumar, and S. Rathore, “Application of homotopy perturbation transform method for solving linear and nonlinear Klein-Gordon equations,” Journal of Information and Computing Science, vol. 7, pp. 131–139, 2012.