Plot r=\sqrt{a^2 b^2}*((1-(b^2*cos^2{\theta})/a^2-(a^2*sin^2{\theta})/b^2)/(1 (b^2*cos^2{\theta})/a^2 (a^2*sin^2{\theta})/b^2))(cos{\theta},sin{\theta})