Orthonormal Bernstein coupled with Block-Pulse Functions for integro-differential form of the singular Emden-
Fowler initial value problems, is considered. First, we introduce the proposed method then the singular differential
equations are transformed to Volterra integro-differential equations. We transform the obtained equations to
algebraic system of equations. Solution to the problem is identified by solving this system and the constructed
solutions are of approximation form. We compare the results with exact solutions and other methods to show the
accuracy of proposed method. The obtained results guarantee the method provides an approximate solution to the
generalized singular Lane-Emden type equations.