A general method for numerically simulating the stochastic time evolution of coupled chemical reactions