The optimization of an optical system benefits greatly from a study of its aberrations and an identification of each of its elements' contribution to the overall aberration figures. The matrix formalism developed by one of the authors was the object of a previous paper and allows the expression of image-space coordinates as high-order polynomials of object-space coordinates. In this paper we approach the question of aberrations, both through the evaluation of the wavefront evolution along the system and its departure from the ideal spherical shape and the use of ray density plots. Using seventh-order matrix modeling, we can calculate the optical path between any two points of a ray as it travels along the optical system and we define the wa...