r,theta,z=var('r'),var('theta'),var('z');
S1=parametric_plot3d([sqrt(4)+4-r^2,r*sin(theta), r*cos(theta)],(r,0,sqrt(4)),(theta,0,2*pi), frame=False, color='blue',opacity=0.7);
S2=parametric_plot3d([sqrt(4)-sqrt(4-r^2),r*sin(theta), r*cos(theta)],(r,0,sqrt(4)),(theta,0,2*pi), frame=False, color='blue',opacity=0.7);
show(S1+S2)