import pykat from pykat.utilities.optics.gaussian_beams import * gx = beam_param(w0=2e-3, z=0) gy = beam_param(w0=1e-3, z=0) beam = HG_beam(gx,gy,0,0) beam.n = 5 beam.m = 6 beam.plot()