phys = [0 2 4];
first_time = 1;
for i=1:length(phys)
    if first_time
        set_param_value('phy',phys(i));
        dynare rbc noclearall;
        first_time = 0;
    else
        set_param_value('phy',phys(i));
        info = stoch_simul(var_list_);
        if info;
          disp(['Computation fails for phy = ' num2str(phy)]);
        end;
    end
end