ENH: Dirichlet distribution

Implements a dirichlet distribution in the stats module.
that would read:

    if np.allclose(np.sum(x, 0), np.ones(x.shape)):

Is that better?