Abstract
This paper uses APL to describe algorithms of classical Runge-Kutta functions up to the eighth-order with and without integration step size control for solving systems of ordinary differential equations of the first order. The results indicate that for high precision the best formula is the seventh-order Runge-Kutta. However, for low accuracy it is best to use the traditional Runge-Kutta formulas (fourth-order).
