import matplotlib print(matplotlib.__version__) # 1.5.1 import numpy as np from scipy.stats import multivariate_normal # for plotting import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D m = 2 #dimension mean = np.zeros(m) sigma = np.eye(m) N = 1000 x1 = np.linspace(-5, 5, N) x2 = np.linspace(-5, 5, N) X1, X2 = np.meshgrid(x1, x2) X = np.c_[np.ravel(X1), np.ravel(X2)] Y_plot = m