Perturbation and Truncation of Probability Generating Function Methods for Stiff Chemical Reactions

Figure 2

Comparison of mean and standard deviation for the Brusselator model. The initial condition and are assumed. (a) shows a typical run of in 500 seconds from a Monte Carlo simulation of the Brusselator model. In (b), denotes the probability that the number of the substrate is . (c) and (d) compare mean and mean ± standard deviations () of , , and obtained from the SSA and the PPGF. The results from the SSA are based on 50,000 realizations.

(a)

(b)

(c)

(d)