The code MEBDFI is an extension of MEBDFDAE for the solution of implicit differential equations and
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.