x, y, z = var('x, y, z')
sfera = implicit_plot3d(x**2 + y**2 + z**2 - 2*x, (x, -3, 3), (y, -3, 3), (z, 0, 3), rgbcolor = 'red')
x1, y1, z1 = var('x1, y1, z1')
cilindr = implicit_plot3d(x1**2 + y1**2 - x1, (x1, -3, 3), (y1, -3, 3), (z1, -3, 3))
xz = sfera + cilindr
xz.show(aspect_ratio=[1,1,1], axes=true)