Solves stiff differential algebraic equations (DAE) using variable stepsize backwards finite difference formula (BDF)