Integration with the radau5-implementation in DotNumerics, over a discontinuity at 20 min.

See also:

• Original radau5 with optimized evaluation of the Jacobian

  Relative tolerance:

  • 0.01
  • 0.001

• Modified radau5 to never skip evaluation of the Jacobian

  Relative tolerance:

  • 0.01
  • 0.001