Added a sample_ellipsoid function

In problems with strong covariances it can speed burn-in quite a bit to initialize walker positions more carefully than emcee.utils.sample_ball allows.  This new function can be used to initalize a collection of walkers with a specified mean and covariance.