This routine solves a system of differential/algebraic equations of the form
g(t,y,y'). Backward differentiation formulas are used to solve for y and
y' as t varies from an initial T to TOUT. This code is based on an
original code by Dr. Linda Petzold, with stepsize selection based on work
by Gustaf Söderlind, and a few other changes.
This software is part of the MATH77/mathc90 collection here.