rmult function

Generate multinomial samples from a probability matrix