Abstract
An iterative method is discussed with respect to its effectiveness and capability of solving singular nonlinear Lane-Emden type equations using reproducing kernel Hilbert space method combined with the Picard iteration. Some new error estimates for application of the method are established. We prove the convergence of the combined method. The numerical examples demonstrates a good agreement between numerical results and analytical predictions.