rdirichlet_mat function

Random generation for multiple Dirichlet distributions