Abstract
The use of automatically computed aberration functions for an exact determination of the geometric illuminance and caustic surfaces is demonstrated. Numerical results for two optical systems studied using aberration coefficients to the thirteenth order are presented.