dynare dsge_humancapital nostrict
phi = 0.1:0.1:1;
first_time = 1;
for i=1:length(phi)
    if first_time
       set_param_value('x',phi(i));
       dynare dsge_humancapital nostrict;
       first_time = 0;
    else
       set_param_value('x',phi(i));
       perfect_foresight_setup;
       perfect_foresight_solver;
    if oo_.deterministic_simulation.status ~= 1;
    disp(['Computation fails for x =' num2str(x)]);
    end;
    end
end
