Gamma random draws
[gb]=gamm_rnd(nrow,ncol,m,k)
* nrow,ncol = the size of the matrix drawn
* m = a parameter such that the mean of the gamma = m/k
* k = a parameter such that the variance of the gamma = m/(k^2)
* r = an nrow x ncol matrix of random numbers from the gamma distribution