ENH: Dirichlet distribution

Implements a dirichlet distribution in the stats module.
Is this logic correct?  Maybe np.allclose would be clearer here.