powered by
rmultin(n, p)
for(i in 1:10) print(rmultin(10, c(1/3, 1/3, 1/3) ))
Run the code above in your browser using DataLab