clear all;
phi_i = [0,0.3,0.6];
results = struct;
for i=1:length(phi_i)
    phi=phi_i(i);
    save phi;
    dynare home_tut33 noclearall;
    results.(sprintf('phi%d',i)) = oo_;
end

%Plotting every variable in one graph
figure;
for j=1:length(var_list_)
    subplot(3,2,j);
    for i=1:length(phi_i)
        plot(results.(sprintf('phi%d',i)).irfs.(sprintf('%s_z',char(var_list_(j)))));
        hold on;
        title(sprintf('%s_z',char(var_list_(j))));
    end
    legend(strcat('phi = ',string(phi_i)));
    hold on;
    legend;
end