clear all; close all; clc;

Phi = [0.25 0.5 0.75 1];
multiplier_zlb_mat = zeros(size(Phi,1),1);
multiplier_nozlb_mat = zeros(size(Phi,1),1);

addpath c:\dynare\4.5.3\matlab
%addpath /Applications/Dynare/4.5.4/matlab

for kk = 1:size(Phi,1)
   
    dynare reis_ni_zlb noclearall;
    
    y_ni = oo_.endo_simul(1,:);
    pi_ni = oo_.endo_simul(2,:);
    i_ni = oo_.endo_simul(3,:);
    rn_ni = oo_.endo_simul(4,:);
    p_ni = oo_.endo_simul(5,:);
    g_ni = oo_.endo_simul(6,:);
    tau_ni = oo_.endo_simul(8,:);
    save g_ni_vals.mat y_ni pi_ni i_ni rn_ni p_ni g_ni;
 
    dynare reis_i_zlb noclearall;
   
    y_i = oo_.endo_simul(1,:);
    pi_i = oo_.endo_simul(2,:);
    i_i = oo_.endo_simul(3,:);
    rn_i = oo_.endo_simul(4,:);
    p_i = oo_.endo_simul(5,:);
    g_i = oo_.endo_simul(6,:);
    tau_i = oo_.endo_simul(8,:);
    save g_i_vals.mat y_i pi_i i_i rn_i p_i g_i;
    
    load g_ni_vals.mat;
    load g_i_vals.mat;
    
    multiplier_zlb = (y_i-y_ni)./(g_i-g_ni);
    multiplier_zlb_mat(kk) = multiplier_zlb(2);

end

%% G multipliers in one plot

f9 = figure(1);
plot( Phi, multiplier_zlb_mat, 'b');
legend('ZLB binding','location','southeast');
%title('Government spending multiplier (SI)');
%saveTightFigure(f9, 'g_multiplier_combined_reis.png');

