Finite difference methods for certain singular two-point boundary value problems

The aim of this paper is to study the outcome of finite difference methods of mainly second order for the numerical solution of certain singular two-point boundary value problems. The method is mainly an extension of the method of Gustafsson (1973) from the linear to nonlinear differential equation and is particularly suitable when the differential equation in question cannot be reduced easily to an elegant system of first-order equations.