import numpy as np mean = [3.5, 3.5] cov = [[0.95, 0.6], [0.6, 0.95]] x, y = np.random.multivariate_normal(mean, cov, 10000).T print(np.mean(x), np.mean(y)) print(np.cov(x, y))