fix to discrete gamma distribution specification.  previous definition gave wrong assignment of each bin's rate, and caused E[r_k] != E[r]