Interpolating numerical solutions of ordinary differential equations

Methods like the Runge-Kutta family for the solution of ordinary differential equations produce approximate solutions only at mesh points. The efficiency of such methods is greatly reduced if the user requests output too frequently. This paper justifies interpolating to resolve this difficulty. In addition the use of interpolation to approximate the derivatives of the solution is justified