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