normal random draws
[y]=norm_rnd(sig)
sig = a square-symmetric covariance matrix
n= # of vectors drawn
y = (nxm) matrix composed of n (1xm) random vector normal draws with mean 0, var-cov(sig)
NOTE: for mean b, var-cov sig use: b + norm_rnd(sig)