In this paper, we establish higher-order numeric solutions for the IVP of the singular Lane-Emden type equation, including the Emden-Fowler equation.We use the multi-stage modified decomposition method to effectively treat these types of equations and develop numeric solutions that are effective in the large.The step-size and the order in our numeric solutions are two parameters that may be arbitrarily specified. Fast algorithms of the Adomian polynomials guarantee the efficiency of our approach, and a higher-order numeric solution can be readily generated at will. The proposed method overcomes the singular behavior at the origin x = 0 and exhibits approximations of high accuracy with a large effective region of convergence. Several numerical examples are examined to demonstrate the reliability of our new approach. In these examples, we have demonstrated A c c e p t e d M a n u s c r i p t that our numeric solutions are consistent by halving the step-size, i.e. the numeric solutions of different step-sizes nearly coincide.