combined rho_mult.R and sim_a.R to simulation.R