plot_complex(gamma(z), (z, -1-1j, 1+1j), coloring="b", n=500,
    grid=False, at_infinity=True, axis=False)
# Expected:
## Plot object containing:
## [0]: complex domain coloring: gamma(1/z) for re(z) over (-1.0, 1.0) and im(z) over (-1.0, 1.0)