Loop over shock values in deterministic simulation

The easiest way is to change the entries of

oo_.exo_simul

and then loop over

perfect_foresight_solver