Stochastic simulation of coupled chemical reactions using recursive methods

In this paper, we present a new method for stochastic simulation of coupled chemical reactions. In this method we obtain recursive expressions for propagating the first two moments of the probability distributions over time. Its advantage over other simulation methods is that it does not require Monte Carlo simulations, and hence it performs several orders of magnitude faster than existing Monte Carlo methods. Simulation results are presented for some examples of coupled first-order reactions.