Cubic spline polynomial for non-linear singular two-point boundary value problems

We present a cubic spline polynomial for the solution of non-linear singular two-point boundary value problems. The quesilinearization technique is used to reduce the non-linear problem to a sequence of linear problems. The resulting sets of differential equations are modified at the singular point and are treated by using cubic spline for finding the numerical solution. The method is tested on three physical model problems from the literature.