plot(1:100, tabulate(rnbinom(10000,3,0.1)+3,100))