Update simulation.R

crtcoef1