Discrete invariant imbedding for singular boundary value problems with a regular singularity

We describe a discrete invariant imbedding method for solving a two point boundary value problem in the interval 0,b] for a linear second order ordinary differential equation with a singularity of the first kind at x = 0. By employing the series expansion on (0, δ], where δ is near the singularity, we first replace it by a regular problem on some interval δ, b]. The discrete invariant imbedding method is then described to solve the problem over the reduced interval. The stability analysis of the method is discussed. Some model problems are solved, and the numerical results compared with those of other methods.