add simulation.R