Update simulation.R

replace Lmatrix with Amatrix