lognormal random draws
[rnd]=logn_rnd(mu,sigma,m,n)
* mu = the mean (may be a matrix)
* sig = the standard deviation (may be a matrix)
* m,n = the size of r in the case of mu or sig a matrix
* rnd = a matrix of random numbers from the lognormal distribution