Ray transfer matrices in NumPy: a beam expander and a laser cavity

Afterwards you can chain ray transfer matrices in NumPy, find focal and image planes, propagate a Gaussian beam, and test a laser cavity for stability.