The program NODE for solution of ODE stiff systems

The paper introduces an algorithm for the numerical solution of initial value problems for systems of ordinary differential equations (ODE). The algorithm uses the Rosenbrock-type and the Runge-Kutta-type schemes with Jacobian freezing and automatic step size control policy based on the global error estimation. Some examples of solution of test problems are given.