are real constants. In recent years, singular boundary value problems for ordinary differential equations have been studied extensively[3, 4, 5, 6, 7, 8, 9, 10, 11, 1, 2, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21] and references therein. However, we find only the following results on coupled Lane-Emden equations.
Recently, in [22, 23, 24] authors studied (1.1) with boundary conditions and that relates the concentration of the carbon substrate and the concentration of oxygen. In [25, 26], authors considered the coupled Lane–Emden equations (1.1) with boundary conditions , and occurs in catalytic diffusion reactions. In [26, 24], the Adomian decomposition method was applied to obtain a convergent analytic approximate solution of (1.1) with . Later, in , the variational iteration method was applied to obtain approximations to solutions of (1.1) for shape factors . In , the Sinc-collocation method was used to obtain the solution of (1.1). In  authors used the reproducing kernel Hilbert space method for solving to obtain the solution of (1.1).
2 Adomian decomposition method
Recently, many researchers [30, 31, 32, 33, 34, 35, 36] have applied the Adomian decomposition method to deal with many different scientific models. According to the Adomian decomposition method we rewrite (1.1) in a operator form as
where and are differential operators and their inverse integral operators are defined as
Operating on (2.1) and using , we get
According to the ADM, we decompose and as
where are Adomian’s polynomials  are given
Upon comparing both sides of (2.6), we have
where are unknown constants to be determined. The -term series solutions are given as
The unknown constants may be obtained by imposing boundary condition at on , that leads to
Solving above transcendental equations for require additional computational work, and may not be uniquely determined.
To avoid solving the above sequence of difficult transcendental equations, the Adomian decomposition method with Green’s function was introduced in [1, 2]. This technique relies on constructing Green’s function before establishing the recursive scheme for the solution components. Unlike the standard Adomian decomposition method, this avoids solving a sequence of transcendental equations for the undetermined coefficients.
3 Green’s function and decomposition technique
In this section, we extend the application of the Adomian decomposition method with Green’s function [2, 12, 37, 38, 39], where we transformed the singular boundary value problem into the integral equation before establishing the recursive scheme for the approximate solution. To apply this technique to coupled Lane-Emden boundary value problems (1.1), we first consider the equivalent integral form of coupled Lane-Emden equation (1.1) as
where are given by
Comparing components from both sides of (3.4) we have the following recursive scheme
Then, we obtain the approximate series solutions as
Unlike ADM or MADM, the proposed recursive schemes (3.6) do not require any computation of unknown constants.
4 Convergence and error analysis
Let be a Banach space with norm
where, and .
The function satisfy Lipschitz condition as
where , and , are Lipschitz constants.
If the approximate solution converges to , then the maximum absolute truncated error is estimated
5 Numerical Results
In this section, we consider three coupled Lane-Emden type boundary value problems to examine the accuracy of the present method. Since the exact solution of the problems is not known, we examine the accuracy and applicability of the present method by the absolute residual error
where is the absolute residual error and is the present approximate solution. The maximum residual errors are defined as
In Tables 1 and 3, we list the numerical results of the approximate solution and the absolute error obtained by the proposed method of Example 5.1 for () and (), respectively. We also compare the numerical results of the maximum residual error obtained by the present method and the results obtained by the modified ADM  in Table 2. In Table 4, we list the numerical results of the maximum residual error.
By applying the proposed scheme (3.5) with the initial guesses , we obtain the 5-terms series solutions as
|Modified ADM |
On applying the proposed scheme (3.5) with the initial guesses , we obtain the 5-terms series solutions as
5.2 The concentration of the carbon substrate and the concentration of oxygen problem 
In Table 5, 6 and 7, we list the numerical results of the approximate solution and the absolute error obtained by the proposed method of Example 5.2 for ), and , respectively. We also compare the numerical results of the maximum residual error obtained by the present method and the results obtained by the modified ADM  in Table 8 for .
5.2.1 When the shape factors
By applying the proposed scheme (3.5) with the initial guesses , we obtain the 4-terms series solutions as