from sympy import *
from spb import *
z = symbols("z")
expr = (z - 1) / (z**2 + z + 1)
graphics(
    riemann_sphere_3d(expr, n=150,
        coloring="b"),
    backend=KB, legend=False, grid=False)