powered by
This routine samples multivarate normal distribution of mean 0 from precision matrix
rmvnorm(n, Omega)
sample size
**precision** matrix of dimension p by p
a matrix with dimension n by p, each row is a sample