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