fanyue322 Update simulation.R