r==Sqrt[Cos[theta]^2+1]