tpronk Update simulation.R