PRGAMMA(x;k;t) probability density function for a Gamma distribution with shape parameter k (kappa) and scale parameter t (theta):
pr(x) = xk-1 × exp(-x/t) / { tk × Γ(k) }.
The mean of the distribution is k × t and the variance is k × t2 .