polar plot r=1/(1 + 0.01671123 * cos theta), theta=0 to 2 pi