covar_normal function

Simulate from multivariate normal distribution