global v var_list_
dynare housing noclearall 

rho_pi_vec = 1.1:0.1:2;
rho_y_vec  = 0.2:0.1:1;
welfare=cell(length(rho_pi_vec),length(rho_y_vec));
first_time = 1;
  for i = 1: length(rho_pi_vec)
      for j = 1: length(rho_y_vec)
     set_param_value('rho_pi',rho_pi_vec(i));
     set_param_value('rho_y',rho_y_vec(j));
        info = stoch_simul(var_list_);
      if info==0
      welfare{i,j} = oo_;
      else
 fprintf('model cannot be solved for rho_pi=%3.2f and rho_y=%3.2f\n',rho_pi_vec(i),rho_y_vec(j))
      end;
    end;
  end;
