import numpy as np

means = np.array([1, 2])
cov = np.array([
	[1, 1],
	[1, 1]
])

X1 = np.random.multivariate_normal(means, cov, )
import numpy as np

means = np.array([1, 2])
cov = np.array([
	[1, 1],
	[1, 1]
])

X1 = np.random.multivariate_normal(means, cov, )