Simulation of multivariate normal distribution