mvgauss function

Simulate from a multivariate normal distribution with mean zero