Spline function approximations for solutions of ordinary differential equations