Dormand-Prince 8(5,3) adaptive Runge-Kutta integrator.
Butcher tableau coefficients follow the Hairer and Wanner DOP853 scheme.