a = 1;

[r,theta] = ellipse_rt(a,0.99);
polar(theta,r)
axis manual
hold on
for e = 0.9 : -0.1 : 0 
    e
    [r,theta] = ellipse_rt(a,e);
    polar(theta,r)
end
hold off