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