Create simulate_population.Rscript