expr = gamma(x)
plot(expr, (x, -5, 5), adaptive=True)
# Expected:
## Plot object containing:
## [0]: cartesian line: gamma(x) for x over (-5.0, 5.0)