In recent years sufficient success was achieved in the modeling and optimization of beamlines on the base of highorder maps. In this paper a new approach based on symbolic representation of high–order aberrations in matrix
forms is described. We discuss all the pros and cons of
such approach.