====== Trapezoidal ======

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|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

===== Reference =====

  * 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.