Kj. Berry et al., AN ALGORITHM TO GENERATE DISCRETE PROBABILITY-DISTRIBUTIONS - BINOMIAL, HYPERGEOMETRIC, NEGATIVE BINOMIAL, INVERSE HYPERGEOMETRIC, AND POISSON, Behavior research methods, instruments, & computers, 26(3), 1994, pp. 366-367
A recursion procedure for generating point probability values, cumulat
ive probability values, and interval probability values of discrete di
stributions is described. Algorithms and associated FORTRAN functions
for five representative discrete probability distributions are present
ed: binomial, hypergeometric, negative binomial, inverse hypergeometri
c, and Poisson. The joint use of recursion and a small arbitrary initi
al value ensures computational efficiency.