pcout2: changed from eigen to svd
simlulation.R: used closed form solution to calculate rho