tau function

Gillespie tau-leap algorithm