%PLOT IRF (IN PERCENTAGE DEVIATION FROM STEADY STATE)

% ./dynare++ --per 320 --sim 1000 version1_logs.mod
% addpath /Applications/Dynare/4.4.2/matlab 
% /Applications/Dynare/4.2.1/dynare++/dynare++ --per 320 --sim 1000 --seed 123 rbc_cost.mod 
% /Applications/Dynare/4.2.1/dynare++/dynare++ --steps 5 rbc_cost_tax_debt.mod 

figure; % new figure
%title(figure,'Productivity shock (EPS>0)');

n = 6;

for i = 1:n
    s(i) = subplot(n,1,i);
end

Zss=dyn_ss(dyn_i_z);
plot(s(1), dyn_irfp_EPS_mean(dyn_i_z,:)/abs(Zss), '-ro');
ylabel(s(1),'delta_z');
set(s(1),'xlim',[0 10]);
set(s(1),'ylim',[-5 5]);

Css=dyn_ss(dyn_i_C);
plot(s(2), 100*dyn_irfp_EPS_mean(dyn_i_C,:)/abs(Css), '-ro');
ylabel(s(2),'delta_C');
set(s(2),'xlim',[0 10]);
%set(s(2),'ylim',[0 2]);

Iss=dyn_ss(dyn_i_I);
plot(s(3), 100*dyn_irfp_EPS_mean(dyn_i_I,:)/abs(Iss), '-ro');
ylabel(s(3),'delta_I');
set(s(3),'xlim',[0 10]);
set(s(3),'ylim',[-10 10]);

Mss=dyn_ss(dyn_i_M);
plot(s(4), dyn_irfp_EPS_mean(dyn_i_M,:)/abs(Mss), '-ro');  %  /abs(Mss)
ylabel(s(4),'M');
set(s(4),'xlim',[0 10]);
set(s(4),'ylim',[-10 0]);

erpss=dyn_ss(dyn_i_erp);
plot(s(5), dyn_irfp_EPS_mean(dyn_i_erp,:)/abs(erpss), '-ro');
ylabel(s(5),'Rex');
set(s(5),'xlim',[0 10]);

qss=dyn_ss(dyn_i_q);
plot(s(6), dyn_irfp_EPS_mean(dyn_i_q,:)/abs(qss), '-ro');
ylabel(s(6),'q');
set(s(6),'xlim',[0 10]);
%set(s(6),'ylim',[0 40]);