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

Figure 11

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=450, 𝑃₂=670, 𝑋₁(0)=𝑋₂(0)=5, 𝐶(0)=50, 𝜌=200, and 𝛽=2/3.

(a)

(b)