Relative tolerance: 0.001
See also:
Original radau5 with optimized evaluation of the Jacobian
Relative tolerance:
Modified radau5 to never skip evaluation of the Jacobian