From: Simulation methods with extended stability for stiff biochemical Kinetics
Reactions
Parameters
Propensities
A + 2 X → k 1 3 X
k1 = 3·10-7
k1x(x - 1)A/2
3 X → k 2 A + 2 X
k2 = 10-4
k2x(x - 1)(x - 2)/6
X → k 3 B
k3 = 3.5
k 3 x
B → k 4 X
k4 = 10-3
k 4 B