decided to use stochpy's NextReactionMethod for simulation