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