The program NODEI for solution of differential-algebraic ODE systems

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