rnb function

Simulate from a negative binomial distribution