In recent papers [6]–[10] the technique for local and global errors estimation and the local-global step size control have been presented to solve both ordinary differential equations and semi-explicit index 1 differential-algebraic systems by multistep methods with any automatically obtained reasonable accuracy. Now we describe the object oriented library INTEGRATOR (ver. 1.0) built in C++ for portability and performance across a wide class of machine architectures. Also we give some computational tests and a comparison with the software implemented the well-known Gear’s method which has been included into MATLAB (ver. 6.2, release 12).