polar plot: r=(1-0.8^2)/(1-(0.8*sin(f))^2)^0.5