The Stochastic Evolution of a Protocell: The Gillespie Algorithm in a Dynamically Varying Volume

Figure 9

Stochastic versus ODE SRM protocell (3). Case of two GMMs, (a) the amount of GMM at the beginning of each division cycle, (b) the division time as a function of the number of elapsed divisions. Parameters are 𝜂₁=𝜂₂=1, 𝛼₁=𝛼₂=2, 𝐿1=500, 𝐿2=600, 𝑃1=600, 𝑃₂=670, 𝑋₁(0)=𝑋₂(0)=100, 𝐶(0)=1000, 𝜌=200, and 𝛽=2/3.

(a)

(b)