The ellipses are now plotted using the fit_ellipse function. Every tenth position the variation results are plotted in black.

See license for source.