master
simulate
def simulate(reaction_stochiometries, kinetic_laws, init_copy_number, time_max, time_step):