plot[(x^2+y^2+f^2-a^2)^2=4f^2x^2+a^4] where a=1, f=sqrt 2