Abstract: A software package for automatic optical design is developed using the theory of canonical coordinates At present, its features include lay-tracing, aberration analysis, calculation of image criteria such as spot diagram, point spread function and optical transfer function, simulation and control of ghost hages and Narcissus effects (for amming infrated systems), damped-least-squares optimization, etc. The software is capable of handling all kinds of general optical systems, i. e., systems without rotational symmetry and/or with unconventional optical elements. Satisfactory results are achieved in practical applications.