Tau-leaping stochastic simulation