add sim_output and simulation.R code