Star Field Simulation
Stot distribution changed to a skewed Gaussian so there are more dim stars than bright stars
Shape function fixed to be more realistic. Adds the effect of an unfocused image.
Also, a globular option was implemented.
Tested the ability of the star finding algorithm with dense star fields
Profiles of Stars
Simulated Data

The sigma for these two stars when using two separate sigma parameters in fit are 8.987 and 9.072, which corresponds to a size of 5.47'' and 5.52''. In comparison, the FWHM corresponds to 6.51''.
This was repeated for real data
The sigma for these stars are 1.8413 and 1.8473 pixels, which correspond to a size of 1.12'' and 1.13''.