“…The main idea of this paper is to present an algorithm for computing the solutions of singular initial value problems including Lane-Emden-type equations of the form (1.1) where α, β, k 1 , k 2 are real constants, a(x), b(x) are continuous and maybe a(0) = 0, b(0) = 0, f (x, y) is a continuous real valued function, and g(x) ∈ C (0, 1], i.e., the case when the function g(x) may be undefined at the origin. Such problems have attracted much attention and have been studied by many authors [3][4][5][6][7][8][9][10][11][12]. There is considerable interest in numerical methods on singular initial value problems.…”