PolarPlot[r == 2/(Sqrt[1 - Sin[2 \[Theta]]] + Sqrt[1 + Sin[2 \[Theta]]])]