global L var_list_
dynare RRendiss noclearall 

RRRss_vec = 0.01:0.01:0.04;
PHI_RR_vec = 0.1:0.1:0.3;
loss=cell(length(RRRss_vec),length(PHI_RR_vec));
first_time = 1;
  for i = 1: length(RRRss_vec)
      for j = 1: length(PHI_RR_vec)
     set_param_value('RRRss',RRRss_vec(i));
     set_param_value('PHI_RR',PHI_RR_vec(j));
        info = stoch_simul(var_list_);
      if info==0
      loss{i,j} = oo_;
      else
 fprintf('model cannot be solved for RRRss=%3.2f and PHI_RR=%3.2f\n!',RRRss_vec(i),PHI_RR_vec(j))
      end;
    end;
  end;
% How to plot?