rmvnorm function

Draw from multivariate normal distribution