gamma_params function

Calculate shape and scale (or rate) parameters of a gamma distribution.