gillespie function

Gillespie exact algorithm