graphics(
    domain_coloring(gamma(z), (z, -3-3j, 3+3j), coloring="b", n=500),
    grid=False)
# Expected:
## Plot object containing:
## [0]: complex domain coloring: gamma(z) for re(z) over (-3.0, 3.0) and im(z) over (-3.0, 3.0)
