Hi,

I am having a problem with the binomial distribution generator routine "bnldev", whenever the integer argument n exceeds a large value, typically 16700000. Often it gets stuck into a seemingly never-ending loop with the variable "t" taking either 0 or a very small value (typically 1e-13 or less), and the condition "while (ran1 (idum) > t)" remains more or less always true (I am using C-version). Did anyone on this forum have this problem, and can someone please point to a solution?

Thanks,
Manojit