This paper proposes an envelope rejection method for computer sampling from the negative binomial distribution, with index parameter α and probability parameter p, based on a simple probability distribution inequality. Comparative timings with existing methods of generating NB variates show that the proposed method is comparable or faster when α is not too large (α < 11), and it is fastest for ...