Internal function to sample from multivariate normal distribution.
mvrnormArma(n, mu, sigma)
Useful links