`Dirichlet::sample()` is incorrect

I noticed that #253 omitted the normalization step when refactoring the `sample()` method of the Dirichlet distribution. The `sum` local variable is not used after it is computed. For large dimensions of alpha, the magnitude of sampled vectors is noticeably shorter than those generated by statrs@0.17.1 or scipy.