Solves stiff ordinary differential equations (ODE) using variable stepsize backwards finite difference formula (BDF)