Runge-Kutta-method implementation of implicit Runge-Kutta method with Richardson Extrapolation for automatic step size control