Code to extract x- and y-extension and rotation (if the beam is elliptical) of a grayscale image taken from a laser beam, by fitting a 2D-Gaussian intensity profile.