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

x = np.random.rand(10,2)
mu, C = bb.MVBE(x)

pl.plot(*x.T, 'o')
bbplot.cov_ellipse(mu, C, fc='none', ec='k')