clear
clc
close all
dynare fondamental.mod
irf_1=oo_.irfs;
dynare fondamentaljojo.mod
irf_2=oo_.irfs;
ending_cell={'_epsA','_epsG','_epsU','_epsR','_epsHC','_epsm','_epsi','_epspi','_epsT','_epsCNR','_epsCR'};
for ii=1:length(ending_cell)
    HOR=1:options_.irf;
    var={'R','Y','R_SPREAD','R_IB','R_B','R_D' ,'PRIM'};
        figure
        for  jj = 1:length(var)
        subplot(4,4,jj)
        eval(['irf_1.' var{1,jj},ending_cell{1,ii}]);
        eval(['irf_2.' var{1,jj},ending_cell{1,ii}]);
        plot(HOR,irf_1.([var{1,jj},ending_cell{1,ii}]),'LineWidth',5);
        plot(HOR,irf_2.([var{1,jj},ending_cell{1,ii}]),'LineWidth',5);
        grid;
        title(['Reponse de  ',var{1,jj},'  Au choc  ',ending_cell{1,ii}] )
        legend('Sans frictions financières','Avec frictions financières')
        xlabel('Trimestres');
        ylabel('IRFs');
    end
end
