RegionPlot[(x^2+y^2)^2<=2*2^2(x^2-y^2),{x,-3,3},{y,-3,3}]