Skip to main content
Figure 2 | BMC Systems Biology

Figure 2

From: Simulation methods with extended stability for stiff biochemical Kinetics

Figure 2

Histogram of X 1 in the Reversible isomerisation reaction. Histogram of X1 in the Reversible isomerisation reaction (106 samples used) solved by the SSA (grey background), Poisson τ-leap (dashed line), and optimal RK τ-leap methods with bounded relative variance. (a) τ = 0.05 (z = -1), Optimal RK τ-leap s = 3, ϵ = 0.1 (solid line) and s = 5, ϵ = 0.1 ("+" marks), (b) τ = 0.4 (z = -8), Optimal RK τ-leap s = 3, ϵ = 0.5 (solid line) and s = 5, ϵ = 0.1 ("+" marks), Poisson τ-leap is unstable for this time step. (c) τ = 0.6 (z = -12), Optimal RK τ-leap s = 5, ϵ = 0.5 ("+" marks), Poisson τ-leap and Optimal RK τ-leap s = 3 are unstable for this time step.

Back to article page