Table 4 The Hybrid Simulation Algorithm (Algorithm 2)

1: initialize the simulation as for Algorithm 1;

2: while t < t _stop do

3:    pick a value for τ as in Algorithm 1;

4:    get θ _next = min{θ _i ∈ Θ | θ _i >t}

5:    if τ < θ _next then

6:       fire a reaction as in Algorithm 1;

7:    else

8:       update (T, E, I(t)) to (T, E + w _next , I(t + θ _next )) and change clock to t + θ _next ;

9:    end if

10: end while