The probability density function for a gamma(a) distributed random variable X, f(x), is:

Function List
 double Gamma_Random_Variate( double shape )
This function returns X where X is a random variate with the density given above. Note that shape must be positive.