powered by
A modified dgamma function
dgamma_l(x, p1, p2, lower, upper, lg = FALSE)
quantile
shape parameter
scale parameter
lower bound
upper bound
log density?