Plot3D[(x^2 + y^2)^3 + 0.001- 4x^2y^2, {x, -1.145, 1.145}, {y, -1.145, 1.145}]