AM_find_gamma_NegBin function

Given that the prior on M is a Negative Binomial, find the γ\gamma hyperparameter of the weights prior to match E(K)=KE(K)=K*, where KK* is user-specified