Nonlinear matrix integration as a mapping approach for ordinary differential equations solving provides powerful technique for particle dynamics simulation. The numerical implementation of this method provides by MODE software package. In this program one can build a matrix map for particles spin-orbit motion in arbitrary electromagnetic fields.