Create simulate.R