The Algorithm

The Monte Carlo method implemented by the algorithm starts by normalizing rates so to have only one parameter: dead_rate, corresponding to a number between 0 and 1, representing a probability.

The process has two parameters: (i) birth and (ii) dead rates. Notice however, that birth has to be bigger than dead or the process goes extinct. So, particles must give birth before they die. This means we can take the time it takes to give birth (1/birth) to be the clock of our system (nothing happens faster than this). Now we measure rates in units of this clock. This way, particles give birth with probability birth/birth = 1 in one dt unit. At the same time, in this same period, particles can die with probability dead_rate = dead/birth.