next up previous
Next: Overview Up: Results (don't read this Previous: Results (don't read this

Revisiting the implementation

Let's think why the simulation underestimate the expected population size at the steady state.

We intended to make a model:

Image immigEmigProcess

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.



Naoki Takebayashi 2008-02-11