まずは2次元の正規分布を描画する。簡便のために対角成分以外は0とする。 import numpy as np from scipy.stats import multivariate_normal import matplotlib.pyplot as plt np.random.seed(2) mu_x,mu_y = 5,10 sigma_x,sigma_y = 1,5 cov = np.array([[1, 0], [0, 1]]) # 共分散行列 n_samples = 100000 samples = np.random.multivariate_normal(np.array([0,…