Many biochemical systems require stochastic descriptions. Unfortunately these can only be solved for the simplest cases and their direct simulation can become prohibitively expensive, precluding thorough analysis. As an alternative, moment closure approximation methods generate equations for the time-evolution of the systems moments and apply a closure ansatz to obtain a closed set of differential equations, that can become the basis for the deterministic analysis of the moments of the outputs of stochastic systems.