Implicit multistep-solver of variable order (1-5) based on the Numerical Differentiation Formulas (NDFs).

Applicable for **stiff** problems of **low** to **medium** accuracy.

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
  * Max order
  * LU decomposition matrix format

====Reference====

* Shampine, L. F. and M. W. Reichelt, “The MATLAB ODE Suite,” _SIAM Journal on Scientific Computing_, Vol. 18, 1997, pp 1-22.

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