networks. We first derive the set of equations for the moments of the master probability
distribution. We then linearize these equations calculating the Jacobian matrix around the
stationary probability distribution. In order to demonstrate the method, we present examples
of stochastic reaction networks and compute their dynamic behavior. We find that the
calculations are accurate and significantly more efficient than stochastic simulation …