A reliable efficient general-purpose method for automatic digital computer integration of systems of ordinary differential equations is described. The method operates with the current values of the higher derivatives of a polynomial approximating the solution. It is thoroughly stable under all circumstances, incorporates automatic starting and automatic choice and revision of elementary interval size, approximately minimizes the amount of computation for a specified accuracy of solution, and applies to any system of differential equations with derivatives continuous or piecewise continuous with finite jumps. ILLIAC library subroutine #F7, University of Illinois Digital Computer Laboratory, is a digital computer program applying this method.
Arnold Nordsieck (1962) studied this question.