dist_gamma function

Gamma distribution