 Simulating a Poisson process with a constant rate formula_18 for the number of events formula_19 that occur in interval formula_20 can be carried out with the following algorithm.Published by Dan Gillespie in 1977, and is a linear search on the cumulative array.See Gillespie algorithm.