add Marginal computation
create a class for sampling from Dirichlet distribution (with GSL)