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.

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.

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.

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 .

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.