Cubic spline and invariant imbedding for solving singular two-point boundary value problems (Q1077895)

A cubic spline approximation combined with an invariant imbedding method is applied to boundary value problems involving a regular singularity at one end of the interval of consideration. The procedure adopted by the paper is now more or less a standard approach for linear ordinary differential equations using a series expansion about the singular point and matched with a smooth solution (in this case spline), away from the singular point. However, the method has an advantage in that the coefficient matrix of the system is tridiagonal and the method has an order of convergence \(O(h^ 4)\), where h is the step size. The paper is very well written with a useful brief account of some numerical treatments of the class of problem considered.