import matplotlib.pyplot as pl
import bluebell as bb
import bluebell.plot as bbplot

D = 3
C = np.eye(D)
mu = np.zeros(D)
std = np.sqrt(np.diag(C))

x = np.random.uniform(low=mu-std, high=mu+std, size=(100,D))
x[0] = mu
chi2 = np.sum(x.dot(np.linalg.inv(C))*x, axis=1)

fig, ax = bbplot.chi2_corner(x, chi2)