from sympy import *
from spb import *
var("x:z")
graphics(
    implicit_3d(
        (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
)