rmvn function

Fast simulation of multivariate normal random variables