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