Numerical analysis of the master equation

Abstract

Applied to the master equation, the usual numerical integration methods, such as Runge-Kutta, become inefficient when the rates associated with various transitions differ by several orders of magnitude. We introduce an integration scheme that remains stable with much larger time increments than can be used in standard methods. When only the stationary distribution is required, a direct iteration method is even more rapid; this method may be extended to construct the quasi-stationary distribution of a process with an absorbing state. Applications to birth-and-death processes reveal gains in efficiency of two or more orders of magnitude.

7 pages 3 figures

References in corpus (1)

• Quasi-stationary distributions for stochastic processes with an absorbing state

Cited by in corpus (10)

• Social consensus through the influence of committed minorities
• Quasi-Stationary Distributions for Models of Heterogeneous Catalysis
• An infinite-period phase transition versus nucleation in a stochastic model of collective oscillations
• Quasi-stationary distributions for the Domany-Kinzel stochastic cellular automaton
• Particle redistribution and slow decay of correlations in hard-core fluids on a half-driven ladder
• On singular probability densities generated by extremal dynamics
• The advantage of being slow: the quasi-neutral contact process
• Absorbing-state phase transitions: exact solutions of small systems
• Topological aspects of the multi-language phases of the Naming Game on community-based networks
• Sampling rare trajectories using stochastic bridges