def mv_normalrand(mu, sigma, size): lamb = cholesky(sigma) return mu + np.dot(lamb, np.random.randn(size))