Improved leap-size selection for accelerated stochastic simulation

In numerically simulating the time evolution of a well-stirred chemically reacting system, the recently introduced “tau-leaping” procedure attempts to accelerate the exact stochastic simulation algorithm by using a special Poisson approximation to leap over sequences of noncritical reaction events. Presented here is an improved procedure for determining the maximum leap size for a specified degree of accuracy.