from sympy import symbols, gamma, I
from spb import *
z = symbols('z')
graphics(
    analytic_landscape(gamma(z), (z, -3 - 3*I, 3 + 3*I)),
    backend=PB, zlim=(-1, 6))