Sample random values from the multivariate normal distribution
sample_mvnorm(mu, sigma)
mean vector
covariance matrix
Samples multivariate normal variables by multiplying univariate random normal variables by the cholesky decomposition of the covariance matrix.
If mu is length 1 then just uses rnorm instead.