The program SIMODE for solution of ODE systems with singular matrix multiplying the derivative

The paper presents an algorithm for the numerical solution of the initial value problems for systems of ordinary differential equations with singular matrix multiplying the derivative. The algorithm uses the (m,k) scheme of the Rosenbrock type with time-lagging derivative matrices, and the adaptive step size control for the global error. Some examples of solution of test problems are presented.