prior_gamma function

Prior gamma distribution