simulation of Riccetti2012 roughly done