p(x) = k * x^(alpha-1) * e^(-x/beta) with k = 1/(g(alpha) * b^a)) and g(a) being the gamma function.
Valid parameter ranges: alpha > 0.
Note: For a Gamma distribution to have the mean mean and variance variance, set the parameters as follows:
alpha = mean*mean / variance; lambda = 1 / (variance / mean);
Instance methods operate on a user supplied uniform random number generator; they are unsynchronized.
Implementation:
J.H. Ahrens, U. Dieter (1974): Computer methods for sampling from gamma, beta, Poisson and binomial distributions, Computing 12, 223-246.
and
J.H. Ahrens, U. Dieter (1982): Generating gamma variates by a modified rejection technique, Communications of the ACM 25, 47-54. @author wolfgang.hoschek@cern.ch @version 1.0, 09/24/99
|
|
|
|