Numerical Algorithm for Solution of the System of Emden–Fowler Type Equations

An efficient numerical approach for solving the system of the Lane–Emden–Fowler type equations with boundary conditions is proposed. To avoid singular behavior, we use the equivalent integral system of equations of the considered boundary value problems for constructing the recursive scheme for the approximate solution. In comparison to the Adomian decomposition method (Singh in Int J Appl Comput Math 6: 80, 2020), the present method involves adjustable parameters to control the convergence of the approximate series solution. Convergence of the new method is provided under quite general conditions. Several standard examples are considered to verify the accuracy of the new algorithm. Computational results reveal that the proposed approach computes the results more accurately than the existing iterative methods.