The nonlinear singular two points boundary value problems arise in many modeling problems in biology, physiology and physics. In this work we propose an iterative method and show its effectiveness and capability of solving two point nonlinear singular boundary value problems arising in physiology. The proposed method is a combination of reproducing kernel Hilbert space method and the shooting technique which takes advantage of two powerful methods for solving nonlinear boundary value problems. An analytical approximated solution is represented in the form of series with easily computable components. Also, the convergence of the proposed method is proved. The method is computationally very simple and attractive. To demonstrate the computational efficiency, the mentioned method is implemented for two problems and the results have been compared with reported results in the literatures.
