  We obtain exact polynomial solutions for two-dimensional coherent complex
scalar fields propagating through arbitrary aberrated shift-invariant linear
imaging systems. These are used to model nodal-line dynamics of coherent fields
output by such systems.
