kaizu

kaizu

kaizu

kaizu

kaizu

kaizu

# A simple implementation of the Gillespie's direct method
import numpy
from functools import partial

def step(t, s, events):
    a = numpy.array([ev[0](s) for ev in events])
    atot = a.sum()
    if atot == 0.0:
        return numpy.inf, s

kaizu

kaizu

kaizu

kaizu