An iterative shooting method for a certain class of singular two-point boundary value problems

Many problems in physics and engineering give rise to singular differential equations. In this paper, we consider a simple shooting method coupled with an iterative method for the numerical solution of a certain class of singular two-point boundary value problems. Our approach enables us to take advantage of such highly accurate initial value solvers as multistep methods. The singularity of the problem is dealt with in a purely algebraic way so that no special care is required when applying an initial value solver. The accuracy of the method is the same as the accuracy of the IVP solver.