Star Field Simulation:

First a random position function was made for a variable canvas size and number of stars. At this point, the stars are just points.

Next, a 2D Gaussian function was used over the position of the star with mean values of x and y = 5. The array was then converted into a fits image.

The same was attempted for twenty stars on a larger canvas size, which worked, but took a long time. This is because the algorithm was looped over every pixel for every star, this can be fixed by limiting the number of pixels to iterate over around each star. The monte carlo part of the simulation hasn't yet been implemented.