bigthetas = 0.95:0.05:0.50;
first_time = 1;
for i=1:length(bigthetas)
    if first_time
        set_param_value('bigtheta',bigthetas(i));
        dynare testmodel noclearall;
        first_time = 0;
    else
        set_param_value('bigtheta',bigthetas(i));
        info = stoch_simul(var_list_);
        if info;
          disp(['Computation fails for bigtheta = ' num2str(bigtheta)]);
        end;
    end
end