Abstract
Traditional chemical kinetic models use ordinary differential equations (ODEs) to predict the concentrations of the involved molecule types. The evolution of the corresponding probability distribution is given by the chemical master equation (CME) which, in principle, can be solved by numerical integration.