function moffat

Create an elliptical moffat function.

Syntax
result = moffat(x, y, power, fwhm1 [, fwhm2 [, angle]] [, naxis1=value] [, naxis2 = value])

Arguments
X:   The center of the moffat in the x-axes
Y:   The center of the moffat in the y-axes
power:   The power index of the moffat function
fwhm1:   Full-width at half maximum of the major axes
fwhm2:   Optional, full-width at half maximum of the minor axes
angle:   Optional, position angle of the moffat
naxis1=value:   Optional, the size of the resulting array in the first dimension
naxis2=value:   Optional, the size of the resulting array in the second dimension