I have equations and the 2D curve drawing in polar coordinates:
Equations and 2D curve drawing in polar coordinates.
var('beta, phi, R, d')
a(phi,beta)=(cos(phi)+beta)/(1+betacos(phi)) # Transformation function
r1(phi,beta,R,d) = da + sqrt(d^2(a^2-1)+R^2) # First equation in polar coordinates
r2(phi,beta,R,d) = da - sqrt(d^2*(a^2-1)+R^2) # second equation
submission of both charts for parameters: beta=0.8, R=1, d=2
polar_plot(r1(phi,0.8,1,2),(-pi/2+0.16,-0.16+pi/2),color='red') + polar_plot(r2(phi,0.8,1,2),(-pi/2+0.18, pi/2-0.18))