simulation from Gillespie 1976 paper
redirect data to file, plot and display