is hybrid of tau-leaping method Runge-Kutta based method LSODA multistep method hybrid method hybrid tau-leaping method 2020-11-06 http://identifiers.org/biosimulators/gillespy2 https://github.com/StochSS/GillesPy2 A continuously coupled hybrid deterministic/stochastic simulation algorithm for biochemical networks. Biochemical species are classified as continuous, discrete, or switch. Tau-leaping is used to simulate stochastic species, and LSODA or another ODE integration method is used to simulate continuous species. Switch species are dynamically classified as either continuous or discrete at each timestep depending on a user defined error tolerance. BD JRK