Numerically solving generated Jacobian equations in freeform optical design

We present an efficient numerical algorithm that can be used to solve the generalized Monge-Ampère equations for a single freeform reflector and lens surface. These equations are instances of so-called ‘generated Jacobian equations’ which are characterized by associated generating functions. The algorithm has a wide applicability to any optical system that can be described by a smooth generating function