How to simulate the quasistationary state

For a large class of processes with an absorbing state, statistical properties of the surviving sample attain time-independent values in the quasistationary (QS) regime. We propose a practical simulation method for studying quasistationary properties, based on the equation of motion governing the QS distribution. In applications to the contact process, the method is shown to reproduce exact results (for the process on a complete graph) and known scaling behavior to high precision. At the critical point, our method is about an order of magnitude more efficient than conventional simulation.