bmixgamma function

Sampling algorithm for mixture of gamma distributions