2×2 matrices compose ray transfer through lenses, mirrors, media; Gaussian beam q(z) transforms as q' = (Aq+B)/(Cq+D).