Question: numeric stiff DAE solving

I have read a few papers on accuracy and numerical methods to solve problems.  My personal interest is in very stiff DAE equations of state (eos) for gaseous expansions/flow problems in very high differential starting conditions for the convergence to occur in the initial time steps.

A long time goal has been to convert the FORTRAN real*16 floating point DAE solver (initially ran on CRAY) into Maple proc() using arbitrary precision that actually solves at the hardware float ieee precision.

Any reference material on stiff solvers and maple implementation would be great.


Please Wait...