The paper presents mathematical and software developments that are the basis for a suite of programs for the solution of initial value problems
, with initial conditions
. The solvers for stiff problems allow the more general form
with a nonsingular and sparse matrix
. The programs are developed for MATLAB, which influences the choice of methods and their implementation.