ENH: Dirichlet distribution

Implements a dirichlet distribution in the stats module.
I tried to be consistent with the majority of the documentation. Then it should be changed from `Multivariate normal distribution` to `Multivariate normal`, but I have no strong opinions on this minor thing.