from sympy import *
from spb import *
var("x:z")
plot3d_implicit(
    (x**2 + y**2 + z**2 >= 1) & (x**2 + y**2 + z**2 <= 3) & (y >= 0),
    (x, -2, 2), (y, -2, 2), (z, -2, 2),
    n=150, backend=KB
)