demo_population_simulation()
run several population simulations
population_simulation()
simulate genetic drift in a finite population