added multinomial to simulate