I have been working on a ray tracing python script as a way to teach myself some vector math and learn something about optics. The
My program calculates how light bounces through a 200mm f/5 Newtonian tube (but currently ignores that the secondary mirror would create a shadow on the primary mirror). My program specifies that light rays start at a random point at the opening of the OTA and is given one of four discrete directions (pointing at the corners of an imaging chip, upper left, upper right, lower left, lower right).
I was expecting the coma to look like the "RealComa' picture (From wikipedia.
Instead I got the plot shown in the 'CalculatedComa' picture. Instead of one cone shape, there appears to be two. Is that the real shape, or is it indicative of an error in my program?