A new efficient recursive numerical scheme is presented for solving a class of singular two-point boundary value problems that arise in various physical models. The approach is based on the homotopy perturbation method in which we establish a recursive scheme without any undetermined coefficients to approximate the singular boundary value problems. The convergence analysis of the present method is discussed. Several numerical examples are provided to show the efficiency of our method for obtaining approximate solutions and to analyze its accuracy. The numerical results reveal that the present method yields a very rapid convergence of the solution without requiring much computational effort. The approximate solution obtained by the present method shows its superiority over existing methods. The Mathematica codes for numerical computation of singular boundary value problems are provided in the paper.