expr = gamma(x)
plot(expr, (x, -5, 5))