DiscreteGamma function

Discrete gamma distribution