Open buggsley opened 7 years ago
This can be done by having the waiting time computed using the following equation:
\int_t^{t+\tau} \alpha_0(t)dt + log(p1) = 0
where p1 is a random number and you solve for \tau.
When you reach \tau, find reaction as follows:
\Sum_{i=1}^{k-1} \alpha_i(t+\tau) < p2 * \alpha0(t+\tau) <= \sum{i=1}^k \alpha_i(t + \tau)
This is being done by Leandro's Java simulator.
Monte Carlo methods currently use an Euler style update for continuous dynamics. Should add support for other ODE update methods to create a variety of hybrid simulators.