Fig. 8
Final distributions of species S2 and S3 in the closed linear system (18) from the SSA, the HR hybrid method, and three remedies (Zero-Population, Zero-Reaction, and Zero-Time) based on 10,000 simulations. The parameters are f1=b1=kc=ks=1. The initial condition is m1=1 and the remaining species populations are zero. a Final population distributions of species S 2b Final population distributions of species S 3c Final population distribution of species S 2 under the Zero-Population remedy d Final population distribution of species S 3 under the Zero-Population remedy