A reported analysis of the stability of some digital simulation methods is investigated by numerical experiments and the results are consistent with the analysis. Traditional stability conditions need to be modified slightly in the presence of homogeneous reactions, though not to a degree that has practical significance. The oscillatory response of the Crank-Nicolson method compared with the Laasonen method can be reduced or eliminated by a preliminary expansion of the time steps, as suggested by Feldberg.