The code MEBDFDAE uses the Modified Extended Backward Differentiation Formulas of Cash, that
increase the absolute stability regions of the classical BDFs. These
methods are A-stable up to the order 4 and stiffly stable for orders up to 9. The orders of the implemented
formulae range from 1 to 8.