Optimized simulation and working on Fig 3