ContourPlot3D[ (3 + x^2 + y^2 + z^2)^2 - 16 (x^2 + y^2) == 0, {x, -3, 3}, {y, -3, 3}, {z, -3, 3}]