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

x = bb.uniform_on_unit_sphere(10, 2)
pl.plot(*x.T, 'o')
bbplot.cov_ellipse(np.zeros(2), np.eye(2), fc='none', ec='k')