Let's think why the simulation underestimate the expected population size at the steady state.
We intended to make a model:
Our implementation was:
popSizeIncrement = 0; if(drand48() < alpha) popSizeIncrement = 1; if(drand48() < beta) popSizeIncrement = popSizeIncrement - 1; n += popSizeIncrement;
Does this correspond to what we intended to do?
Exercises: Make the corrections, run the simulation, and compare with the analytical expectation.