This paper is devoted to one of the methods of symbolic computation, which based on perturbation theory and constructed in the framework of matrix formalism. This method can be used for solving optimization problems and preliminary modeling in accelerator physics. There are presented the main theoretical positions and demonstration the result of this method on example of one of classical problems.