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