dynare model1.mod
    load('model1_results.mat', 'oo_')
    irf1=oo_.irfs;
    save irf1
   


    dynare model2.mod 
    load('model2_results.mat', 'oo_')
    irf2=oo_.irfs;
    save irf2

     load irf1
     load irf2
    
    ending_cell={'_xi_M'};

    for ii=1:length(ending_cell)
HOR=1:options_.irf;
var={'c', 'y', 'i', 'q', 's', 'pi', 'pi_h', 'pi_f', 'psi_f'};
figure

for jj=1:length(var)
subplot(1,2,jj)
eval(['irf1.' var{1,jj},ending_cell{1,ii}]);
eval(['irf2.' var{1,jj},ending_cell{1,ii}]);
hold on
plot(HOR,[eval(['irf1.' var{1,jj},ending_cell{1,ii}])],'-k',HOR,[eval(['irf2.' var{1,jj},ending_cell{1,ii}])],'LineWidth',2.5,'--r') ;
title([var{1,jj}] )
end
end