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