Update simulation.R