Homotopy perturbation method with Laplace Transform (LT-HPM) for solving Lane–Emden type differential equations (LETDEs)

In this communication, we describe the Homotopy Perturbation Method with Laplace Transform (LT-HPM), which is used to solve the Lane–Emden type differential equations. It’s very difficult to solve numerically the Lane–Emden types of the differential equation. Here we implemented this method for two linear homogeneous, two linear nonhomogeneous, and four nonlinear homogeneous Lane–Emden type differential equations and use their appropriate comparisons with exact solutions. In the current study, some examples are better than other existing methods with their nearer results in the form of power series. The Laplace transform used to accelerate the convergence of power series and the results are shown in the tables and graphs which have good agreement with the other existing method in the literature. The results show that LT-HPM is very effective and easy to implement.


Boundary conditions are
Thus, the solutions describe the gallop of pressure and density (with radius), which is known as polytropes of n.
The Lane-Emden equation has been useful to model some phenomena in astrophysics and mathematical physics such as the principle of stellar structure, the thermal nature of the spherical bunch of the gas, isothermal gas spheres (IGSs), and the principle of thermionic currents (Wazwaz 2011). According to the extensive study of many physicists, these equations have been applicable in the case of astrophysics such as kinetics of combustion and the Landau-Ginzburg major phenomenon (Dixon and Tuszynski 1990;Fermi 1927;Fowler 1930;Frank-Kamenetskii 1921). The numerical solutions of the Lane-Emden equations [LEes] are very difficult due to the expressive nature of the nonlinearities term. Therefore, much attention has been applied to the better and more powerful methods for establishing a solution, approximate or exact, analytical or numerical, to the Lane-Emden equations (LEes).
Recently many analytical techniques have been used for the solution of Lane-Emden type equation, for example, Hosseini and Abbasbandy (2015) described the hybrid Spectral Adomain Decomposition Method for solving Lane-Emden type of differential equations by combining the spectral method and Adomain Decomposition method. Since the description of this method is very long and difficult for solving these types of numerical problems. Our method is suitable and best to determine these results, Modified Laplace decomposition method for Lane-Emden type differential equations by Yin et al. (2013), A new algorithm for solving singular IVPs of Lane-Emden type differential equation by Motsa and Sibanda (2010). The author has solved the Lane-Emden type equations by Successive Linearization Method Since it is a very complicated method to get the solution of these types of problem in terms of an exact solution. Various methods for Lane-Emden equations have described by some authors in Rafiq et al. (2009), Baranwal et al. (2012), Liao (2003), Shawagfeh (1993), Wazwaz (2001), A new method for solving singular IVPs in the second order ordinary differential equations by Wazwaz (2001), Nouh (2004), Romas (2003) and other researchers have been studied several methods to attempt nonlinear problems. These methods have also been successfully applied to, analytical solution of convection-diffusion problem by combining Laplace transform method and homotopy perturbation method by Gupta et al. (2015), Mandelzweig and Quasi (2001), Singh et al. (2012), Nazari-Golshan et al. (2013, Explicit solution of Helmholtz equation and sixth-order KdV equation by Rafei and Ganji (2006), Ganji and Rajabi (2006), Jang (2016), Exact solutions of some coupled nonlinear partial differential equations(NPDE) using the homotopy perturbation method by Sweilam and Khader (2009), Homotopy perturbation method for solving viral dynamical model (VDM) by Merdan and Khaniyev (2010), nonlinear population dynamics models(NPDMs) by Chowdhury and Hashim (2007, Hashim andChowdhury (2007), A new dispersion-relation preserving method for integrating the classical boussinesq equation by Jang (2017), The modified homotopy perturbation method (MHPM) for solving strongly nonlinear oscillators by Momani et al. (2009), Pandit (2014 and pure nonlinear differential by Cveticanin (2006), Inverse problem (IP) of diffusion equation by He's homotopy perturbation method by Shakeri and Dehghan (2007), A Higher order Numerical Scheme for singularly perturbed Burger-Huxley equation by Jiwari and Mittal (2011).
The Laplace transform is a superb technique for solving linear and nonlinear Lane-Emden type differential equation and has enjoyed much success in the field of science and engineering. On the other hand, Laplace Transform (LT) has played an important role in mathematics (Spiegel and Teoríay 1988), not only for its theoretical interest but also because such method allows solving, in a simpler fashion, many problems in the realm of science, in comparison with other mathematical techniques. It is totally difficult to solve nonlinear equations because of the problems caused by nonlinear terms. Homotopy perturbation technique by He (1999), the homotopy perturbation method using Laplace Transform by Madani et al. (2011;Abbasbandy (2006); Gupta and Gupta 2011), a numerical solution of two-point boundary value problems using Galerkin-Finite element method by Sharma et al. (2012) have solved nonlinear problems. The Homotopy perturbation methods with Laplace transform (LT-HPM) and other methods have external significant thought in the literature. Moreover, The Homotopy Perturbation Method (HPM) by He (1999aHe ( , b, 2003 and the Variational Iteration Method (VIM) (Khuri and Sayfy 2012) are combined with the Laplace transform (LT) to develop a more effective technique for handling many nonlinear problems. A comparative study of model of matrix and finite elements methods for two-point boundary value problems is given by Sharma et al. (2012).
In the present paper, a Homotopy Perturbation Method (HPM) with Laplace Transform (LT) to solve the general type of Lane-Emden differential equations is proposed; the paper is organized as follows: The Homotopy Perturbation method is given in "Preliminaries" section. The Lane-Emden Equations (LEes) is given in "The Land-Emden equation" section. Homotopy Perturbation Method with Laplace Transform (LT-HPM) is given in "Homotopy perturbation method" section. Some examples of a different kind are given in "Results and discussion" section. Finally, the conclusion is explained in "Conclusions" section.

Homotopy perturbation method
Consider the nonlinear differential equation with the boundary conditions of where A, B, ζ(r) and χ are a general differential operator, a boundary operator, a known analytic function and the boundary of the domain Ω, respectively and ∂Γ ∂n denotes the differentiation of Γ with respect to n.
We can distribute the operator A into a linear part K and a nonlinear part N. Equation (1) may possibly be written as: (3) By the Homotopy technique, we construct a Homotopy v(r, p) : Ω × [0, 1] → R which satisfies: p∈ [0, 1] is an embedding parameter, while Γ 0 is an initial approximation which satisfies the boundary conditions. Obviously, from above equations, we will have Now changing the process of p from 0 to 1 is even-handed that of v(r, p) from Γ 0 to Γ (r) . According to the concept of topology, this is called deformation, whereas K (v) − K (Γ 0 ) and A(v) − ζ(r) are called homotopy. If we consider the embedding parameter p is a minor parameter, applying the classical perturbation technique, we can assume that the solution of Eqs. (4) and (5) can be defined as a power series in p: The coupling of the perturbation method and the homotopy method is said to be HPM. The series (6) is convergent for most cases. However, the convergent rate depends on the nonlinear operator A(v). Moreover, He (1999) made the following suggestions: (1) The second derivative of N (v) with respect to v must be minor because the parameter may be comparatively large, i.e. p → 1.
(2) The norm of K −1 ∂N ∂v must be lesser than one so that the series converges.

Laplace transform method
Definition The Laplace transform of a function ξ(τ), is defined by (Whenever integral on RHS exists) where, τ ≥ 0, s is real and L is the Laplace transform operator.

The Lane-Emden equation
Lane-Emden type differential equations are singular initial value problems (IVPs) describing the second order homogeneous and nonhomogeneous linear and nonlinear differential equations which have been applicable in the many fields. The mathematical representation of Lane-Emden equation is: subject to conditions, where A and B are constants ζ (Γ ) is a real-valued continuous function.
These types of equations generally occur in the principle of stellar structure, the thermal behaviour of a spherical bunch of gas, isothermal gas spheres (IGSs) and the principle of thermionic currents (Richardson 1921;Chandrasekhar 1967;Davis 1962).
On the other hand, A nonlinear class of singular initial value problems of Lane-Emden type has the following form: The solution of the Lane-Emden type differential equation is numerically challenging because of the singularity behaviour at the origin. The solutions of the Lane-Emden equation were given by Wazwaz (2001), Shawagfeh (1993), Mishra (2014), the homotopy perturbation method (HPM) of above type by Davis (1962), Yildrim and Ozis (2007), He (2003He ( ), 2006, Ramos (2008), Exact solution of Generalized Lane-Emden equation is given by Goenner and Havas (2000). The major advantage of this method is its capability of combining the two powerful methods to obtain exact solutions of nonlinear equations. Therefore, Homotopy Perturbation Method using Laplace transform (LT-HPM) accelerates the rapid convergence of the series solution. In this paper, we will apply the (LT-HPM) to obtain exact or approximate analytical solutions of the Lane-Emden type equations.

Homotopy perturbation method with laplace transform (LT-HPM)
In this section, we will briefly discuss the use of the LT-HPM for the solution of Lane-Emden equation given in "The Land-Emden equation" section, consider the following: Multiplying τ and then taking the Laplace transform on both sides of (11) we get: where L is the operator of Laplace transform and L ′ (Γ ) = dL(Γ ) ds by integrating both sides of (13) with respect to s, we have taking inverse Laplace transform on both sides of (14), we get By using initial condition (9), we have and N [Γ (τ )] denote the linear term and the nonlinear term respectively. The Homotopy perturbation method and the He's polynomials can be used to handle Eq. (12) and to address the nonlinear term. LT-HPM defines a solution by an infinite series of components given by: where the terms Γ n (τ ) are to recursively calculate and the nonlinear term ζ (Γ ) can be given as where N (Γ ) is a non-linear term and H n (Γ ) is He's polynomial.
For some He's polynomial H n (Mishra 2012) that are given by Substituting the value of (19) and (20) in (18), we get which is the coupling of the Laplace transformation and the Homotopy Perturbation Method (LT-HPM) using He's polynomials by Mishra (2012Mishra ( , 2014. Comparing the coefficient of like powers of p, the following approximations are obtain

Results and discussion
In this section, we will apply the method presented in this paper to solve singular IVPs of Lane-Emden-type.

Fig. 1 Comparison between LT-HPM, Exact and ADM solution
The closed form of the series (28) is y = exp(τ 2 ) which gives an approximate solution of the problem.
In this problem, the Table 1 shows the comparison of values of LT-HPM with exact and ADM in the terms of different values of τ (Fig. 1).
The graphical comparison of LT-HPM, exact solution and ADM solution are given as follows: Example 2 Consider the linear, homogeneous Lane-Emden differential equation subject to the initial condition Taking Laplace transform on both sides, we get putting n = 0 in Eq. (31), we get Integrating above equation with respect to s, we get Taking inverse Laplace Transform on both sides, we get which is the exact solution.
The Table 2 shows the comparison of values of LT-HPM within the terms of different values of τ.
Thus we get the solution of this series as follows: The closed form of the series (36) is Γ (τ ) = τ 2 + τ 3 which gives an exact solution of the problem.
The comparison with exact solution is given by Table 3. The Table 3 shows the comparison of values of LT-HPM with exact in the terms of different values of τ.
The closed form of the series (40), is Γ (τ ) = τ 2 which gives an exact solution of the problem.
The comparison with exact solution is given by Table 4. The Table 4 shows the comparison of values of LT-HPM with exact in the terms of different values of τ.
Example 5 Consider the linear homogeneous differential equation with the initial conditions Applying the Laplace transform on both sides, we get By integrating both sides with respect to s, we get Applying inverse Laplace Transformation on both sides, we get Applying HPM on both sides, we get ds.
Which is given by using the value of Eq. (50) in Eq. (49) and equating the coefficient of like the power of p, we get Thus we get the solution of this series as follows: The closed form of the series (52) is Γ (τ ) = −2In(1 + τ 2 ) which gives an exact solution of the problem.
The closed form of the series (57) is Γ (τ ) = exp(τ 2 ) which gives an exact solution of the problem.
The Table 5 shows the comparison of values of LT-HPM with exact in the terms of different values of τ.
Example 8 Consider the following Lane-Emden type differential equation: subject to the initial condition Applying Laplace Transform (LT) on both sides, we get Integrating the above equation with respect to s, we get Applying ILT on both sides, we get Applying HPM on both sides, we get Here ζ (Γ ) = sin(Γ ) is a non-linear term and H n (Γ ) is He's polynomial. Γ (s) = − s −2 ds + s −2 ((L(τ sin(Γ ))))ds.