Description of DASSL: a differential/algebraic system solver

Abstract

This paper describes a new code DASSL, for the numerical solution of implicit systems of differential/algebraic equations. These equations are written in the form F(t,y,y') = 0, and they can include systems which are substantially more complex than standard form ODE systems y' = f(t,y). Differential/algebraic equations occur in several diverse applications in the physical world. We outline the algorithms and strategies used in DASSL, and explain some of the features of the code. In addition, we outline briefly what needs to be done to solve a problem using DASSL.