sample_mvnorm function

Sample random values from the multivariate normal distribution