from sympy import I, symbols, exp, sqrt, cos, sin, pi, gamma
from spb import *
z = symbols("z")
expr = (z - 1) / (z**2 + z + 2)
graphics(riemann_sphere_2d(expr, coloring="b", n=800), grid=False)
# Expected:
## Plot object containing:
## [0]: complex domain coloring: (z - 1)/(z**2 + z + 2) for re(z) over (-1.25, 1.25) and im(z) over (-1.25, 1.25)
## [1]: parametric cartesian line: (cos(t), sin(t)) for t over (0.0, 6.283185307179586)