ENH: Dirichlet distribution

Implements a dirichlet distribution in the stats module.
since alpha is ensured to be an array here, can we use np.min instead of min ?