In this paper, we use the concept of homotopy, Laplace transform, and He’s polynomials, to propose the auxiliary Laplace homotopy parameter method (ALHPM). We construct a homotopy equation consisting on two auxiliary parameters for solving nonlinear differential equations, which switch nonlinear terms with He’s polynomials. The existence of two auxiliary parameters in the homotopy equation allows us to guarantee the convergence of the obtained series. Compared with numerical techniques, the method solves nonlinear problems without any discretization and is capable to reduce computational work. We use the method for different types of singular Emden–Fowler equations. The solutions, constructed in the form of a convergent series, are in excellent agreement with the existing solutions.

1. Introduction

Initial value problems with singularity and of type Lane–Emden differential equation,have been used to model a large category of phenomena in various science, such as mathematical physics and astrophysics. The first studies on these equations have been published by Lane in 1870 [1]. A further research was continued by Emden in [2]. Then, this equation was used in the modeling of some problems, such as the thermal behavior of spherical cloud of gas.

In astrophysics and in the study of a self-gravitating spherically symmetric polytropic fluid, this equation appears as its gravitational Poisson’s equation. There are several phenomena, such as astrophysics, aerodynamics, stellar structure, chemistry, biochemistry, and many others which can be modeled by the Lane–Emden equation [35]. Fowler [6, 7] generalized the Lane–Emden equation to the Emden–Fowler equation:for some given functions and . In this connection, we note that the following heat equation,where , and , appears in the modeling of the diffusion of heat perpendicular to surfaces of parallel planes. For and and for the steady state, equation (3) is the Emden–Fowler equation. Singularity behavior that occurs at is a difficult element in the analysis of this type of equations, and due to this problem, common methods need to be reviewed. For example, Wazwaz [8] used the Adomian decomposition method to solve these equations; however, the appropriate choices of operator were required to overcome the singularity behavior at the origin. In this paper and in our proposed method, taking the Laplace transform from both sides of the equation solves this difficulty. Indeed, in vast majority of cases in differential equations with variable coefficients, we cannot obtain an exact solution, so we must look for approximate solutions, such as asymptotic techniques [9, 10], analytical [1114], and numerical methods [1519]. Emden–Fowler-type equations have been solved using the Adomian polynomials [20] by Wazwaz in [8], using the homotopy perturbation method by Chowdhury and Hashim in [21], and using the variational iteration method by Wazwaz in [22]. In this paper, we use a new auxiliary homotopy parameter method using He’s polynomials, called the auxiliary Laplace homotopy parameter method (ALHPM) for solving Emden–Fowler equations.

2. Main Results

In this section, we present the new auxiliary Laplace homotopy parameter method (ALHPM). For this purpose, let us consider the following nonlinear nonhomogeneous PDE:which subject towhere , is a general nonlinear term, which may include nonlinear differential operators, and denotes the source term. Using the Laplace transform [23] and by its applying to the both sides of (4), we deduce

This quickly yields that

Now, we use the homotopy concept to deduce the following homotopy equation of two parameters:where and are nonzero auxiliary parameters, is the Laplace transform, and is an initial guess of the solution. The result of applying Laplace inverse on both sides of (8) is as follows:where

Now, we use the basis of the homotopy perturbation method to have a series expansion for as follows:where is an embedding parameter. To switch nonlinear operator , we use He’s polynomials [24] to obtainwhere

Substituting (11) and (12) in (9) and comparing the coefficient of like powers of , we obtain

3. Emden–Fowler Equations

Due to the singularity behavior at the origin, as well as other various linear and nonlinear singular IVPs, numerical solution of Emden–Fowler equations is a challenging issue. For numerical treatment, some authors have been forced to propose alternative approaches. These techniques often focus on the removal of the singularity of this equation. For example, by expanding the unknown function as different basis functions, the problem reduces to a set of algebraic equations to using operational matrices. However, it is clear that equivalent numerical techniques are more computationally expensive. So, often analytical methods have been considered. In this section, in order to show the efficiency of the auxiliary Laplace homotopy parameter method, presented in the previous section, here, we use this method to solve the initial value problems related to second order singular Emden–Fowler differential equations. The examples show that our method leads to the exact solution series of the problem, and in this case, it is possible to guess the closed form of the solution.

Example 3.1. Let us consider the following nonlinear singular Lane–Emden equation:withUsing the Laplace transform and by applying to the both sides of (15), we deducewhere . As the ALHPM technique, using the concept of the homotopy perturbation method and He’s polynomials for nonlinear term and constructing two parameters homotopy of equation (17), we get the following:and the first few components of He’s polynomials areUsing (12) and by comparing the same powers of , it is clear that the recursive relation isTherefore, we haveand so on. Consequently, using the Mathematica symbolic code, the series solution of (13) and (14) is given byBy putting , we get a series providing a closed form of the exact solution:

Example 3.2. Now, in this example, we consider the following Emden–Fowler equation:withwhich has been investigated in several studies from various points of view because of its interesting mathematical and physical properties.
Using the Laplace transform and by applying to the both sides of (24), we deducewhere . By constructing two parameters’ homotopy of equation (26) and applying the aforesaid method, we get the following:where are He’s polynomials. The first few components of He’s polynomials are as follows:Comparing the similar powers of , we deduceso that, using Laplace inverse, we obtainSubstituting in equation (30), one haswhich is the same as that obtained in [8]. We can easily show that if , one can obtain the exact solution of (24) from (31). As an example, substituting and into (31), we obtainwhich is the exact solution of the following equation:Similarly, for and , we obtain a series as a closed form ofwhich is the exact solution of the equation:More choices of and can be easily checked. However, for the values of and which the obtained series (31) is not defined, the solution of (21) must be investigated separately. If and , equation (24) is the well-known Euler equation and the solution can be obtained easily. For , , (and generally ), the exact solution of (21) can be represented by the Bessel functions of the first and the second kinds.

Example 3.3. Now, let us suppose and , so we considerwithSimilar to the previous examples, using the Laplace transform and by its applying to the both sides of (36), we deduceIn a similar way as above, one writesBy equating the powers of , one can obtainTherefore, we have the following series solution:Putting in equation (41) recovers the solution obtained in [8]. Then,In the direction of method efficiency, we consider the case and . The numerical results are compared with the solutions obtained via collocation scheme in Table 1. On the contrary, it is possible to improve the results of the series solution by using the Pade’ approximations. We have applied the Pade’ approximation to the obtained series solution (ALHPMP [11]).
The accuracy is obtained when the computational cost of our method is much less than the comparative method.

4. Conclusion

The new auxiliary Laplace homotopy parameter method was presented. Using the concept of homotopy and Laplace transformation, a two-parameter depending homotopy equation was constructed. Nonlinear terms were dealt with by He’s polynomials. Using of two auxiliary parameters in homotopy equation enables us to make the solutions more reliable. We used this method to solve the initial value problems related to second-order singular Emden–Fowler differential equations. In these equations, the singularity behavior at the origin prevents us from using the usual methods, and to overcome this problem, these methods must be changed. In our proposed method, this limitation is easily removed. Obtained solutions, constructed in the form of a convergent series, were in excellent agreement with the existing solutions.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

All authors contributed equally in writing this article. All authors read and approved the final manuscript.