Update simulate.R