from sympy import *
from spb import *
z, u, a, b = symbols("z u a b")
graphics(
    complex_vector_field(
        log(gamma(u * z)), prange(z, -5*a - b*5j, 5*a + b*5j),
        params={
            u: (1, 0, 2),
            a: (1, 0, 2),
            b: (1, 0, 2)
        }, quiver_kw=dict(color="orange", headwidth=4)),
    n=20, grid=False)