Resolution:
## Figure 2.
Comparative simulation results for a reversible reaction with feedback controls. In all panels, the x-axis denotes time in seconds and the y-axis represents the number of molecules of species S_{1}. The upper and lower panels use two different sets of initial numbers of molecules,
namely: (x_{1}(0), x_{2}(0), x_{3}(0), U) = (5, 5, 6, 1μm^{3}) and (x_{1}(0), x_{2}(0), x_{3}(0), U) = (100, 100, 120, 20μm^{3}), respectively. Other simulation parameters are (f_{1}, f_{2}, f_{3}, g_{1}, g_{3}, k, _{f}k) = (1.3, 1.8, -1, 1, 1, 0.5, 0.5). In both the upper and lower panels, the first
column compares the time evolution of _{g}S_{1 }molecules by different methods: the black line shows the ODE solution of Equation
(52) for x_{1 }; the blue lines are the solutions of Equation (53) for μ_{1 }and for μ_{1 }± σ_{1}, respectively. The red dotted lines framing the mean indicate the 95% confidence
interval. The second column shows the propensity adjustment functions for the forward
reaction (solid line) and the backward reaction (dashed line). The third column shows
100 independent stochastic simulations with propensity adjustment (blue means and
error bars), in comparison with the ODE (Equation (52)) prediction (black line). The
fourth column shows a second set of 100 independent stochastic simulations without
propensity adjustment (blue means and error bars), in comparison with the ODE (Equation
(52)) prediction (black line). The red dotted lines framing the mean in columns 3
and 4 again indicate the 95% confidence intervals.
Wu |