This model uses a standard elliptical exponential form for the galaxy where the ratio of the major to the minor axis is constant (before shear bias). The function centered around (xo,yo) is (note that it is not normalized):

f(x,y) = c*exp(-sqrt(((x-xo)/a)^2+((y-yo)/b)^2)) where c=mag*magnitude

"mag" is a random parameter that varies with each galaxy. "magnitude" is a user-provided constant overall intensity scaling factor. The value of "a" is the major axis length provided by the user. The value of "b" is given by b=a*e, where "e" is the axis ratio provided by the user.

The function is rasterized by sampling once per pixel.

A set of random transformations are applied to the function. These are:
mag = .2+getRand()*2.0
theta = pi*getRand()
First the model is rotated in the plane by theta, and finally every value is multiplied by mag. The random values are reselected every time drawGal() is called unless "SelectGal 0" has been specified as a parameter, in which case the calling function must call select() when it wants to choose new random values. Either of these random values can be made non-random by specifying the proper parameter as specified below.

Parameters

Name	Description
MajorAxis	The length of the major axis in pixels -- the value of "a" in the equation.
AxisRatio	The fixed ratio of the major to the minor axis. The value of "e" in the equation.
Magnitude	An overall intensity scaling factor in units of counts.
VaryRotation	If set to 0, never vary rotation of galaxy (it will be set to 0.0).
VaryMagnitude	If set to 0, never vary the mag parameter (it will be set to 1.0).
SelectGal	If set to 0, do not automatically change random parameters before drawing.