Spline approximations for systems of ordinary differential equations