Implementation of two-step Runge-Kutta methods for ordinary differential equations