powered by
Random multivariate normal (standard N(0,1) margins)
rmvn(n,rmat)
nxd matrix, where d=nrow(rmat)
simulation sample size
correlation matrix