Implicit one-step Runge-Kutta solver implementation of the trapezoidal rule using a “free” interpolant.

Applicable for moderately stiff problems of low accuracy where you need a solution without numerical damping.

Applicable solver settings:

• Absolute tolerance
• Relative tolerance
• Initial step size
• Max step size
• Min step size
• Refine
• Limit data points to last
• Norm Control
• Allowed step size violations
• Enable saturation
• LU decomposition matrix format

• Shampine, L. F., M. W. Reichelt, and J.A. Kierzenka, “Solving Index-1 DAEs in MATLAB and Simulink,” SIAM Review, Vol. 41, 1999, pp 538-552.