! with fixed-leading-coefficient implementation. ! this version is in double precision. ! dvode solves the initial value problem for stiff or nonstiff ! systems of first order odes, ! dy/dt = f(t,y) , ...