Loop in perfect foresight context

That depends on how you compute the steady state. In principle, you can just change the first and last entry of

oo_.endo_simul
and then loop over

perfect_foresight_solver