Spline Function Approximations for Solutions of Ordinary Differential Equations

A procedure for obtaining spline function approximations for solutions of the initial value problem in ordinary differential equations is presented. The proposed method with quadratic and cubic splines is shown to be related to the well-known trapezoidal rule and Milne-Simpson method, respectively. The method is shown to be divergent, however, when higher degree spline functions are used.