@Test
public void testRmultinom2() {
double[] b = {0.1, 0.2, 0.7};
Vector v = new DenseVector(b);
Vector t = v.like();
for (int i = 1; i <= 100; i++) {
Vector multinom = UncommonDistributions.rMultinom(100, v);
t = t.plus(multinom);
}
System.out.println("sum(rMultinom(" + 100 + ", [ 0.1, 0.2, 0.7 ]))/100="