%% Technology Shock in Germany

%First Scenario (Country-Variable-Scenario) 
load Buba_Fiskal_Erweiterung_FU_baseline_results;
G_Y_1   = oo_.irfs.dgdp2_a_t_nua_a;
G_C_1   = oo_.irfs.c2_a_obs_nua_a;
G_Pi_1  = oo_.irfs.pi_a_t_nua_a;
G_K_1   = oo_.irfs.k_a_t_nua_a;
G_B_1   = oo_.irfs.BG_a_t_nua_a;
G_R_1   = oo_.irfs.i_policy_t_nua_a;
G_In_1  = oo_.irfs.in_a_obs_nua_a;
G_U_1   = oo_.irfs.ur_a_obs_nua_a;
G_W_1   = oo_.irfs.wr_a_obs_nua_a;
G_RER_1 = oo_.irfs.rer_ba_obs_nua_a;

Eu_Y_1  = oo_.irfs.dgdp2_b_t_nua_a;
Eu_C_1  = oo_.irfs.c2_b_obs_nua_a;
Eu_Pi_1 = oo_.irfs.pi_b_t_nua_a;
Eu_K_1  = oo_.irfs.k_b_t_nua_a;
Eu_B_1  = oo_.irfs.BG_b_t_nua_a;
Eu_In_1 = oo_.irfs.in_b_obs_nua_a;
Eu_U_1  = oo_.irfs.ur_b_obs_nua_a;
Eu_W_1  = oo_.irfs.wr_b_obs_nua_a;
save G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1

%Second Scenario (Country-Variable-Scenario) 
load Buba_Fiskal_Erweiterung_FU_Eq_results;
G_Y_2   = oo_.irfs.dgdp2_a_t_nua_a;
G_C_2   = oo_.irfs.c2_a_obs_nua_a;
G_Pi_2  = oo_.irfs.pi_a_t_nua_a;
G_K_2   = oo_.irfs.k_a_t_nua_a;
G_B_2   = oo_.irfs.BG_a_t_nua_a;
G_R_2   = oo_.irfs.i_policy_t_nua_a;
G_In_2  = oo_.irfs.in_a_obs_nua_a;
G_U_2   = oo_.irfs.ur_a_obs_nua_a;
G_W_2   = oo_.irfs.wr_a_obs_nua_a;
G_RER_2 = oo_.irfs.rer_ba_obs_nua_a;


Eu_Y_2  = oo_.irfs.dgdp2_b_t_nua_a;
Eu_C_2  = oo_.irfs.c2_b_obs_nua_a;
Eu_Pi_2 = oo_.irfs.pi_b_t_nua_a;
Eu_K_2  = oo_.irfs.k_b_t_nua_a;
Eu_B_2  = oo_.irfs.BG_b_t_nua_a;
Eu_In_2 = oo_.irfs.in_b_obs_nua_a;
Eu_U_2  = oo_.irfs.ur_b_obs_nua_a;
Eu_W_2  = oo_.irfs.wr_b_obs_nua_a;

save G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2

%Third Scenario (Country-Variable-Scenario) 
load Buba_Fiskal_Erweiterung_FU_TH_results;
G_Y_3   = oo_.irfs.dgdp2_a_t_nua_a;
G_C_3   = oo_.irfs.c2_a_obs_nua_a;
G_Pi_3  = oo_.irfs.pi_a_t_nua_a;
G_K_3   = oo_.irfs.k_a_t_nua_a;
G_B_3   = oo_.irfs.BG_a_t_nua_a;
G_R_3   = oo_.irfs.i_policy_t_nua_a;
G_In_3  = oo_.irfs.in_a_obs_nua_a;
G_U_3   = oo_.irfs.ur_a_obs_nua_a;
G_W_3   = oo_.irfs.wr_a_obs_nua_a;
G_RER_3 = oo_.irfs.rer_ba_obs_nua_a;

Eu_Y_3  = oo_.irfs.dgdp2_b_t_nua_a;
Eu_C_3  = oo_.irfs.c2_b_obs_nua_a;
Eu_Pi_3 = oo_.irfs.pi_b_t_nua_a;
Eu_K_3  = oo_.irfs.k_b_t_nua_a;
Eu_B_3  = oo_.irfs.BG_b_t_nua_a;
Eu_In_3 = oo_.irfs.in_b_obs_nua_a;
Eu_U_3  = oo_.irfs.ur_b_obs_nua_a;
Eu_W_3  = oo_.irfs.wr_b_obs_nua_a;
save G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3

load Buba_Fiskal_Erweiterung_FU_CA_results;
G_Y_4   = oo_.irfs.dgdp2_a_t_nua_a;
G_C_4   = oo_.irfs.c2_a_obs_nua_a;
G_Pi_4  = oo_.irfs.pi_a_t_nua_a;
G_K_4   = oo_.irfs.k_a_t_nua_a;
G_B_4   = oo_.irfs.BG_a_t_nua_a;
G_R_4   = oo_.irfs.i_policy_t_nua_a;
G_In_4  = oo_.irfs.in_a_obs_nua_a;
G_U_4   = oo_.irfs.ur_a_obs_nua_a;
G_W_4   = oo_.irfs.wr_a_obs_nua_a;
G_RER_4 = oo_.irfs.rer_ba_obs_nua_a;

Eu_Y_4  = oo_.irfs.dgdp2_b_t_nua_a;
Eu_C_4  = oo_.irfs.c2_b_obs_nua_a;
Eu_Pi_4 = oo_.irfs.pi_b_t_nua_a;
Eu_K_4  = oo_.irfs.k_b_t_nua_a;
Eu_B_4  = oo_.irfs.BG_b_t_nua_a;
Eu_In_4 = oo_.irfs.in_b_obs_nua_a;
Eu_U_4  = oo_.irfs.ur_b_obs_nua_a;
Eu_W_4  = oo_.irfs.wr_b_obs_nua_a;
save G_Y_4 G_C_4 G_Pi_4 G_K_4 G_B_4 G_R_4 G_In_4 G_U_4 G_W_4 G_RER_4 Eu_Y_4 Eu_C_4 Eu_Pi_4 Eu_K_4 Eu_B_4 Eu_In_4 Eu_U_4 Eu_W_4


load G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2
load G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1
load G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3
load G_Y_4 G_C_4 G_Pi_4 G_K_4 G_B_4 G_R_4 G_In_4 G_U_4 G_W_4 G_RER_4 Eu_Y_4 Eu_C_4 Eu_Pi_4 Eu_K_4 Eu_B_4 Eu_In_4 Eu_U_4 Eu_W_4
 
% Figure with Impulse Responses after a technology shock
T=size(G_Y_1,2);
t=1:T;
Zeroline=zeros(size(t,2),1);
f = figure;
set(f,'InvertHardcopy','off','Color',[1 1 1],'Position',[1 1 800 1000]);
subplot (4,4,1);
   plot(t,G_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
 subplot (4,4,2); 
   plot(t,G_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_C_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_C_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,3);
   plot(t,G_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_In_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_In_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,4);
   plot(t,G_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('GE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,5);
   plot(t,G_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_U_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_U_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('GE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,6);
   plot(t,G_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_W_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_W_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on    
   plot(t,Zeroline,'k','LineWidth',1);
   title('GW Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,7); 
   plot(t,G_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_B_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_B_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,8);
   plot(t,Eu_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_B_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_B_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,9); 
  plot(t,Eu_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,10); 
   plot(t,Eu_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_C_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_C_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,11); 
   plot(t,Eu_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_In_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_In_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,12); 
   plot(t,Eu_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,13); 
   plot(t,Eu_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_U_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_U_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,14); 
   plot(t,Eu_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_W_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_W_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,15);
   plot(t,G_RER_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_RER_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_RER_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE ToT to DE','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,16);
   plot(t,G_R_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_R_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_R_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('Interest Rate','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
   h_legend=legend('SQ','RE','TH','CA','Orientation','horizontal','Location','SouthOutside');
%    legend('boxoff');
    pos = get(h_legend,'Position');
    pos(1) = 0.253;
    pos(2) = 0.038;
    pos(3) = 0.53;
    pos(4) = 0.033;
    set(h_legend,'Position',pos);
    set(h_legend,'FontSize',13,'FontWeight','Bold');


%%
% 

load Buba_Fiskal_Erweiterung_FU_baseline_results;
G_Y_1   = oo_.irfs.dgdp2_a_t_nub_a;
G_C_1   = oo_.irfs.c2_a_obs_nub_a;
G_Pi_1  = oo_.irfs.pi_a_t_nub_a;
G_K_1   = oo_.irfs.k_a_t_nub_a;
G_B_1   = oo_.irfs.BG_a_t_nub_a;
G_R_1   = oo_.irfs.i_policy_t_nub_a;
G_In_1  = oo_.irfs.in_a_obs_nub_a;
G_U_1   = oo_.irfs.ur_a_obs_nub_a;
G_W_1   = oo_.irfs.wr_a_obs_nub_a;
G_RER_1 = oo_.irfs.rer_ba_obs_nub_a;

Eu_Y_1  = oo_.irfs.dgdp2_b_t_nub_a;
Eu_C_1  = oo_.irfs.c2_b_obs_nub_a;
Eu_Pi_1 = oo_.irfs.pi_b_t_nub_a;
Eu_K_1  = oo_.irfs.k_b_t_nub_a;
Eu_B_1  = oo_.irfs.BG_b_t_nub_a;
Eu_In_1 = oo_.irfs.in_b_obs_nub_a;
Eu_U_1  = oo_.irfs.ur_b_obs_nub_a;
Eu_W_1  = oo_.irfs.wr_b_obs_nub_a;
save G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1


load Buba_Fiskal_Erweiterung_FU_Eq_results;
G_Y_2   = oo_.irfs.dgdp2_a_t_nub_a;
G_C_2   = oo_.irfs.c2_a_obs_nub_a;
G_Pi_2  = oo_.irfs.pi_a_t_nub_a;
G_K_2   = oo_.irfs.k_a_t_nub_a;
G_B_2   = oo_.irfs.BG_a_t_nub_a;
G_R_2   = oo_.irfs.i_policy_t_nub_a;
G_In_2  = oo_.irfs.in_a_obs_nub_a;
G_U_2   = oo_.irfs.ur_a_obs_nub_a;
G_W_2   = oo_.irfs.wr_a_obs_nub_a;
G_RER_2 = oo_.irfs.rer_ba_obs_nub_a;


Eu_Y_2  = oo_.irfs.dgdp2_b_t_nub_a;
Eu_C_2  = oo_.irfs.c2_b_obs_nub_a;
Eu_Pi_2 = oo_.irfs.pi_b_t_nub_a;
Eu_K_2  = oo_.irfs.k_b_t_nub_a;
Eu_B_2  = oo_.irfs.BG_b_t_nub_a;
Eu_In_2 = oo_.irfs.in_b_obs_nub_a;
Eu_U_2  = oo_.irfs.ur_b_obs_nub_a;
Eu_W_2  = oo_.irfs.wr_b_obs_nub_a;

save G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2

 
load Buba_Fiskal_Erweiterung_FU_TH_results;
G_Y_3   = oo_.irfs.dgdp2_a_t_nub_a;
G_C_3   = oo_.irfs.c2_a_obs_nub_a;
G_Pi_3  = oo_.irfs.pi_a_t_nub_a;
G_K_3   = oo_.irfs.k_a_t_nub_a;
G_B_3   = oo_.irfs.BG_a_t_nub_a;
G_R_3   = oo_.irfs.i_policy_t_nub_a;
G_In_3  = oo_.irfs.in_a_obs_nub_a;
G_U_3   = oo_.irfs.ur_a_obs_nub_a;
G_W_3   = oo_.irfs.wr_a_obs_nub_a;
G_RER_3 = oo_.irfs.rer_ba_obs_nub_a;

Eu_Y_3  = oo_.irfs.dgdp2_b_t_nub_a;
Eu_C_3  = oo_.irfs.c2_b_obs_nub_a;
Eu_Pi_3 = oo_.irfs.pi_b_t_nub_a;
Eu_K_3  = oo_.irfs.k_b_t_nub_a;
Eu_B_3  = oo_.irfs.BG_b_t_nub_a;
Eu_In_3 = oo_.irfs.in_b_obs_nub_a;
Eu_U_3  = oo_.irfs.ur_b_obs_nub_a;
Eu_W_3  = oo_.irfs.wr_b_obs_nub_a;

save G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3


load Buba_Fiskal_Erweiterung_FU_CA_results;
G_Y_4   = oo_.irfs.dgdp2_a_t_nub_a;
G_C_4   = oo_.irfs.c2_a_obs_nub_a;
G_Pi_4  = oo_.irfs.pi_a_t_nub_a;
G_K_4   = oo_.irfs.k_a_t_nub_a;
G_B_4   = oo_.irfs.BG_a_t_nub_a;
G_R_4   = oo_.irfs.i_policy_t_nub_a;
G_In_4  = oo_.irfs.in_a_obs_nub_a;
G_U_4   = oo_.irfs.ur_a_obs_nub_a;
G_W_4   = oo_.irfs.wr_a_obs_nub_a;
G_RER_4 = oo_.irfs.rer_ba_obs_nub_a;

Eu_Y_4  = oo_.irfs.dgdp2_b_t_nub_a;
Eu_C_4  = oo_.irfs.c2_b_obs_nub_a;
Eu_Pi_4 = oo_.irfs.pi_b_t_nub_a;
Eu_K_4  = oo_.irfs.k_b_t_nub_a;
Eu_B_4  = oo_.irfs.BG_b_t_nub_a;
Eu_In_4 = oo_.irfs.in_b_obs_nub_a;
Eu_U_4  = oo_.irfs.ur_b_obs_nub_a;
Eu_W_4  = oo_.irfs.wr_b_obs_nub_a;

save G_Y_4 G_C_4 G_Pi_4 G_K_4 G_B_4 G_R_4 G_In_4 G_U_4 G_W_4 G_RER_4 Eu_Y_4 Eu_C_4 Eu_Pi_4 Eu_K_4 Eu_B_4 Eu_In_4 Eu_U_4 Eu_W_4

load G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3
load G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2
load G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1
 

T=size(G_Y_1,2);
t=1:T;
Zeroline=zeros(size(t,2),1);
f = figure;
set(f,'InvertHardcopy','off','Color',[1 1 1],'Position',[1 1 800 1000]);
subplot (4,4,1);
   plot(t,G_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
 subplot (4,4,2); 
   plot(t,G_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,3);
   plot(t,G_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,4);
   plot(t,G_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,5);
   plot(t,G_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,6);
   plot(t,G_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,7); 
   plot(t,G_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,8);
   plot(t,Eu_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,9); 
  plot(t,Eu_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,10); 
   plot(t,Eu_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,11); 
   plot(t,Eu_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,12); 
   plot(t,Eu_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,13); 
   plot(t,Eu_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,14); 
   plot(t,Eu_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,15);
   plot(t,G_RER_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE ToT to Germany','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,16);
   plot(t,G_R_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('Interest Rate','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
   h_legend=legend('SQ','RE','TH','CA','Orientation','horizontal','Location','SouthOutside');
%    legend('boxoff');
    pos = get(h_legend,'Position');
    pos(1) = 0.253;
    pos(2) = 0.038;
    pos(3) = 0.53;
    pos(4) = 0.033;
    set(h_legend,'Position',pos);
    set(h_legend,'FontSize',13,'FontWeight','Bold');
   
%% 
% 

load Buba_Fiskal_Erweiterung_FU_baseline_results;
G_Y_1   = oo_.irfs.dgdp2_a_t_nua_eM;
G_C_1   = oo_.irfs.c2_a_obs_nua_eM;
G_Pi_1  = oo_.irfs.pi_a_t_nua_eM;
G_K_1   = oo_.irfs.k_a_t_nua_eM;
G_B_1   = oo_.irfs.BG_a_t_nua_eM;
G_R_1   = oo_.irfs.i_policy_t_nua_eM;
G_In_1  = oo_.irfs.in_a_obs_nua_eM;
G_U_1   = oo_.irfs.ur_a_obs_nua_eM;
G_W_1   = oo_.irfs.wr_a_obs_nua_eM;
G_RER_1 = oo_.irfs.rer_ba_obs_nua_eM;

Eu_Y_1  = oo_.irfs.dgdp2_b_t_nua_eM;
Eu_C_1  = oo_.irfs.c2_b_obs_nua_eM;
Eu_Pi_1 = oo_.irfs.pi_b_t_nua_eM;
Eu_K_1  = oo_.irfs.k_b_t_nua_eM;
Eu_B_1  = oo_.irfs.BG_b_t_nua_eM;
Eu_In_1 = oo_.irfs.in_b_obs_nua_eM;
Eu_U_1  = oo_.irfs.ur_b_obs_nua_eM;
Eu_W_1  = oo_.irfs.wr_b_obs_nua_eM;
save G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1


load Buba_Fiskal_Erweiterung_FU_Eq_results;
G_Y_2   = oo_.irfs.dgdp2_a_t_nua_eM;
G_C_2   = oo_.irfs.c2_a_obs_nua_eM;
G_Pi_2  = oo_.irfs.pi_a_t_nua_eM;
G_K_2   = oo_.irfs.k_a_t_nua_eM;
G_B_2   = oo_.irfs.BG_a_t_nua_eM;
G_R_2   = oo_.irfs.i_policy_t_nua_eM;
G_In_2  = oo_.irfs.in_a_obs_nua_eM;
G_U_2   = oo_.irfs.ur_a_obs_nua_eM;
G_W_2   = oo_.irfs.wr_a_obs_nua_eM;
G_RER_2 = oo_.irfs.rer_ba_obs_nua_eM;


Eu_Y_2  = oo_.irfs.dgdp2_b_t_nua_eM;
Eu_C_2  = oo_.irfs.c2_b_obs_nua_eM;
Eu_Pi_2 = oo_.irfs.pi_b_t_nua_eM;
Eu_K_2  = oo_.irfs.k_b_t_nua_eM;
Eu_B_2  = oo_.irfs.BG_b_t_nua_eM;
Eu_In_2 = oo_.irfs.in_b_obs_nua_eM;
Eu_U_2  = oo_.irfs.ur_b_obs_nua_eM;
Eu_W_2  = oo_.irfs.wr_b_obs_nua_eM;

save G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2

 
load Buba_Fiskal_Erweiterung_FU_TH_results;
G_Y_3   = oo_.irfs.dgdp2_a_t_nua_eM;
G_C_3   = oo_.irfs.c2_a_obs_nua_eM;
G_Pi_3  = oo_.irfs.pi_a_t_nua_eM;
G_K_3   = oo_.irfs.k_a_t_nua_eM;
G_B_3   = oo_.irfs.BG_a_t_nua_eM;
G_R_3   = oo_.irfs.i_policy_t_nua_eM;
G_In_3  = oo_.irfs.in_a_obs_nua_eM;
G_U_3   = oo_.irfs.ur_a_obs_nua_eM;
G_W_3   = oo_.irfs.wr_a_obs_nua_eM;
G_RER_3 = oo_.irfs.rer_ba_obs_nua_eM;

Eu_Y_3  = oo_.irfs.dgdp2_b_t_nua_eM;
Eu_C_3  = oo_.irfs.c2_b_obs_nua_eM;
Eu_Pi_3 = oo_.irfs.pi_b_t_nua_eM;
Eu_K_3  = oo_.irfs.k_b_t_nua_eM;
Eu_B_3  = oo_.irfs.BG_b_t_nua_eM;
Eu_In_3 = oo_.irfs.in_b_obs_nua_eM;
Eu_U_3  = oo_.irfs.ur_b_obs_nua_eM;
Eu_W_3  = oo_.irfs.wr_b_obs_nua_eM;

save G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3


load Buba_Fiskal_Erweiterung_FU_CA_results;
G_Y_4   = oo_.irfs.dgdp2_a_t_nua_eM;
G_C_4   = oo_.irfs.c2_a_obs_nua_eM;
G_Pi_4  = oo_.irfs.pi_a_t_nua_eM;
G_K_4   = oo_.irfs.k_a_t_nua_eM;
G_B_4   = oo_.irfs.BG_a_t_nua_eM;
G_R_4   = oo_.irfs.i_policy_t_nua_eM;
G_In_4  = oo_.irfs.in_a_obs_nua_eM;
G_U_4   = oo_.irfs.ur_a_obs_nua_eM;
G_W_4   = oo_.irfs.wr_a_obs_nua_eM;
G_RER_4 = oo_.irfs.rer_ba_obs_nua_eM;

Eu_Y_4  = oo_.irfs.dgdp2_b_t_nua_eM;
Eu_C_4  = oo_.irfs.c2_b_obs_nua_eM;
Eu_Pi_4 = oo_.irfs.pi_b_t_nua_eM;
Eu_K_4  = oo_.irfs.k_b_t_nua_eM;
Eu_B_4  = oo_.irfs.BG_b_t_nua_eM;
Eu_In_4 = oo_.irfs.in_b_obs_nua_eM;
Eu_U_4  = oo_.irfs.ur_b_obs_nua_eM;
Eu_W_4  = oo_.irfs.wr_b_obs_nua_eM;

save G_Y_4 G_C_4 G_Pi_4 G_K_4 G_B_4 G_R_4 G_In_4 G_U_4 G_W_4 G_RER_4 Eu_Y_4 Eu_C_4 Eu_Pi_4 Eu_K_4 Eu_B_4 Eu_In_4 Eu_U_4 Eu_W_4

load G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3
load G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2
load G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1
 
% 
T=size(G_Y_1,2);
t=1:T;
Zeroline=zeros(size(t,2),1);
f = figure;
set(f,'InvertHardcopy','off','Color',[1 1 1],'Position',[1 1 800 1000]);
subplot (4,4,1);
   plot(t,G_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
 subplot (4,4,2); 
   plot(t,G_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,3);
   plot(t,G_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,4);
   plot(t,G_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,5);
   plot(t,G_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_4,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,6);
   plot(t,G_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,7); 
   plot(t,G_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,8);
   plot(t,Eu_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,9); 
  plot(t,Eu_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,10); 
   plot(t,Eu_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,11); 
   plot(t,Eu_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,12); 
   plot(t,Eu_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,13); 
   plot(t,Eu_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,14); 
   plot(t,Eu_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,15);
   plot(t,G_RER_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE ToT to Germany','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,16);
   plot(t,G_R_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('Interest Rate','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
   h_legend=legend('SQ','RE','TH','CA','Orientation','horizontal','Location','SouthOutside');
%    legend('boxoff');
    pos = get(h_legend,'Position');
    pos(1) = 0.253;
    pos(2) = 0.038;
    pos(3) = 0.53;
    pos(4) = 0.033;
    set(h_legend,'Position',pos);
    set(h_legend,'FontSize',13,'FontWeight','Bold');
    
    
%% 
% 

load Buba_Fiskal_Erweiterung_FU_baseline_results;
G_Y_1   = oo_.irfs.dgdp2_a_t_eps_y_c;
G_C_1   = oo_.irfs.c2_a_obs_eps_y_c;
G_Pi_1  = oo_.irfs.pi_a_t_eps_y_c;
G_K_1   = oo_.irfs.k_a_t_eps_y_c;
G_B_1   = oo_.irfs.BG_a_t_eps_y_c;
G_R_1   = oo_.irfs.i_policy_t_eps_y_c;
G_In_1  = oo_.irfs.in_a_obs_eps_y_c;
G_U_1   = oo_.irfs.ur_a_obs_eps_y_c;
G_W_1   = oo_.irfs.wr_a_obs_eps_y_c;
G_RER_1 = oo_.irfs.rer_ba_obs_eps_y_c;

Eu_Y_1  = oo_.irfs.dgdp2_b_t_eps_y_c;
Eu_C_1  = oo_.irfs.c2_b_obs_eps_y_c;
Eu_Pi_1 = oo_.irfs.pi_b_t_eps_y_c;
Eu_K_1  = oo_.irfs.k_b_t_eps_y_c;
Eu_B_1  = oo_.irfs.BG_b_t_eps_y_c;
Eu_In_1 = oo_.irfs.in_b_obs_eps_y_c;
Eu_U_1  = oo_.irfs.ur_b_obs_eps_y_c;
Eu_W_1  = oo_.irfs.wr_b_obs_eps_y_c;
save G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1


load Buba_Fiskal_Erweiterung_FU_Eq_results;
G_Y_2   = oo_.irfs.dgdp2_a_t_eps_y_c;
G_C_2   = oo_.irfs.c2_a_obs_eps_y_c;
G_Pi_2  = oo_.irfs.pi_a_t_eps_y_c;
G_K_2   = oo_.irfs.k_a_t_eps_y_c;
G_B_2   = oo_.irfs.BG_a_t_eps_y_c;
G_R_2   = oo_.irfs.i_policy_t_eps_y_c;
G_In_2  = oo_.irfs.in_a_obs_eps_y_c;
G_U_2   = oo_.irfs.ur_a_obs_eps_y_c;
G_W_2   = oo_.irfs.wr_a_obs_eps_y_c;
G_RER_2 = oo_.irfs.rer_ba_obs_eps_y_c;


Eu_Y_2  = oo_.irfs.dgdp2_b_t_eps_y_c;
Eu_C_2  = oo_.irfs.c2_b_obs_eps_y_c;
Eu_Pi_2 = oo_.irfs.pi_b_t_eps_y_c;
Eu_K_2  = oo_.irfs.k_b_t_eps_y_c;
Eu_B_2  = oo_.irfs.BG_b_t_eps_y_c;
Eu_In_2 = oo_.irfs.in_b_obs_eps_y_c;
Eu_U_2  = oo_.irfs.ur_b_obs_eps_y_c;
Eu_W_2  = oo_.irfs.wr_b_obs_eps_y_c;

save G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2


load Buba_Fiskal_Erweiterung_FU_TH_results;
G_Y_3   = oo_.irfs.dgdp2_a_t_eps_y_c;
G_C_3   = oo_.irfs.c2_a_obs_eps_y_c;
G_Pi_3  = oo_.irfs.pi_a_t_eps_y_c;
G_K_3   = oo_.irfs.k_a_t_eps_y_c;
G_B_3   = oo_.irfs.BG_a_t_eps_y_c;
G_R_3   = oo_.irfs.i_policy_t_eps_y_c;
G_In_3  = oo_.irfs.in_a_obs_eps_y_c;
G_U_3   = oo_.irfs.ur_a_obs_eps_y_c;
G_W_3   = oo_.irfs.wr_a_obs_eps_y_c;
G_RER_3 = oo_.irfs.rer_ba_obs_eps_y_c;

Eu_Y_3  = oo_.irfs.dgdp2_b_t_eps_y_c;
Eu_C_3  = oo_.irfs.c2_b_obs_eps_y_c;
Eu_Pi_3 = oo_.irfs.pi_b_t_eps_y_c;
Eu_K_3  = oo_.irfs.k_b_t_eps_y_c;
Eu_B_3  = oo_.irfs.BG_b_t_eps_y_c;
Eu_In_3 = oo_.irfs.in_b_obs_eps_y_c;
Eu_U_3  = oo_.irfs.ur_b_obs_eps_y_c;
Eu_W_3  = oo_.irfs.wr_b_obs_eps_y_c;

save G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3


load Buba_Fiskal_Erweiterung_FU_CA_results;
G_Y_4   = oo_.irfs.dgdp2_a_t_eps_y_c;
G_C_4   = oo_.irfs.c2_a_obs_eps_y_c;
G_Pi_4  = oo_.irfs.pi_a_t_eps_y_c;
G_K_4   = oo_.irfs.k_a_t_eps_y_c;
G_B_4   = oo_.irfs.BG_a_t_eps_y_c;
G_R_4   = oo_.irfs.i_policy_t_eps_y_c;
G_In_4  = oo_.irfs.in_a_obs_eps_y_c;
G_U_4   = oo_.irfs.ur_a_obs_eps_y_c;
G_W_4   = oo_.irfs.wr_a_obs_eps_y_c;
G_RER_4 = oo_.irfs.rer_ba_obs_eps_y_c;

Eu_Y_4  = oo_.irfs.dgdp2_b_t_eps_y_c;
Eu_C_4  = oo_.irfs.c2_b_obs_eps_y_c;
Eu_Pi_4 = oo_.irfs.pi_b_t_eps_y_c;
Eu_K_4  = oo_.irfs.k_b_t_eps_y_c;
Eu_B_4  = oo_.irfs.BG_b_t_eps_y_c;
Eu_In_4 = oo_.irfs.in_b_obs_eps_y_c;
Eu_U_4  = oo_.irfs.ur_b_obs_eps_y_c;
Eu_W_4  = oo_.irfs.wr_b_obs_eps_y_c;

save G_Y_4 G_C_4 G_Pi_4 G_K_4 G_B_4 G_R_4 G_In_4 G_U_4 G_W_4 G_RER_4 Eu_Y_4 Eu_C_4 Eu_Pi_4 Eu_K_4 Eu_B_4 Eu_In_4 Eu_U_4 Eu_W_4

load G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3
load G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2
load G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1
 
% 
T=size(G_Y_1,2);
t=1:T;
Zeroline=zeros(size(t,2),1);
f = figure;
set(f,'InvertHardcopy','off','Color',[1 1 1],'Position',[1 1 800 1000]);
subplot (4,4,1);
   plot(t,G_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
 subplot (4,4,2); 
   plot(t,G_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,3);
   plot(t,G_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,4);
   plot(t,G_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,5);
   plot(t,G_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_4,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,6);
   plot(t,G_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,7); 
   plot(t,G_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,8);
   plot(t,Eu_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,9); 
  plot(t,Eu_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,10); 
   plot(t,Eu_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,11); 
   plot(t,Eu_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,12); 
   plot(t,Eu_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,13); 
   plot(t,Eu_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,14); 
   plot(t,Eu_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,15);
   plot(t,G_RER_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE ToT to Germany','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,16);
   plot(t,G_R_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_3,'g-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_4,'m-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('Interest Rate','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
   h_legend=legend('SQ','RE','TH','CA','Orientation','horizontal','Location','SouthOutside');
%    legend('boxoff');
    pos = get(h_legend,'Position');
    pos(1) = 0.253;
    pos(2) = 0.038;
    pos(3) = 0.53;
    pos(4) = 0.033;
    set(h_legend,'Position',pos);
    set(h_legend,'FontSize',13,'FontWeight','Bold');
    
    
%% 


load Buba_Fiskal_Erweiterung_FU_baseline_results;
G_Y_1   = oo_.irfs.dgdp2_a_t_nua_eb;
G_C_1   = oo_.irfs.c2_a_obs_nua_eb;
G_Pi_1  = oo_.irfs.pi_a_t_nua_eb;
G_K_1   = oo_.irfs.k_a_t_nua_eb;
G_B_1   = oo_.irfs.BG_a_t_nua_eb;
G_R_1   = oo_.irfs.i_policy_t_nua_eb;
G_In_1  = oo_.irfs.in_a_obs_nua_eb;
G_U_1   = oo_.irfs.ur_a_obs_nua_eb;
G_W_1   = oo_.irfs.wr_a_obs_nua_eb;
G_RER_1 = oo_.irfs.rer_ba_obs_nua_eb;

Eu_Y_1  = oo_.irfs.dgdp2_b_t_nua_eb;
Eu_C_1  = oo_.irfs.c2_b_obs_nua_eb;
Eu_Pi_1 = oo_.irfs.pi_b_t_nua_eb;
Eu_K_1  = oo_.irfs.k_b_t_nua_eb;
Eu_B_1  = oo_.irfs.BG_b_t_nua_eb;
Eu_In_1 = oo_.irfs.in_b_obs_nua_eb;
Eu_U_1  = oo_.irfs.ur_b_obs_nua_eb;
Eu_W_1  = oo_.irfs.wr_b_obs_nua_eb;
save G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1


load Buba_Fiskal_Erweiterung_FU_Eq_results;
G_Y_2   = oo_.irfs.dgdp2_a_t_nua_eb;
G_C_2   = oo_.irfs.c2_a_obs_nua_eb;
G_Pi_2  = oo_.irfs.pi_a_t_nua_eb;
G_K_2   = oo_.irfs.k_a_t_nua_eb;
G_B_2   = oo_.irfs.BG_a_t_nua_eb;
G_R_2   = oo_.irfs.i_policy_t_nua_eb;
G_In_2  = oo_.irfs.in_a_obs_nua_eb;
G_U_2   = oo_.irfs.ur_a_obs_nua_eb;
G_W_2   = oo_.irfs.wr_a_obs_nua_eb;
G_RER_2 = oo_.irfs.rer_ba_obs_nua_eb;


Eu_Y_2  = oo_.irfs.dgdp2_b_t_nua_eb;
Eu_C_2  = oo_.irfs.c2_b_obs_nua_eb;
Eu_Pi_2 = oo_.irfs.pi_b_t_nua_eb;
Eu_K_2  = oo_.irfs.k_b_t_nua_eb;
Eu_B_2  = oo_.irfs.BG_b_t_nua_eb;
Eu_In_2 = oo_.irfs.in_b_obs_nua_eb;
Eu_U_2  = oo_.irfs.ur_b_obs_nua_eb;
Eu_W_2  = oo_.irfs.wr_b_obs_nua_eb;

save G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2


load Buba_Fiskal_Erweiterung_FU_TH_results;
G_Y_3   = oo_.irfs.dgdp2_a_t_nua_eb;
G_C_3   = oo_.irfs.c2_a_obs_nua_eb;
G_Pi_3  = oo_.irfs.pi_a_t_nua_eb;
G_K_3   = oo_.irfs.k_a_t_nua_eb;
G_B_3   = oo_.irfs.BG_a_t_nua_eb;
G_R_3   = oo_.irfs.i_policy_t_nua_eb;
G_In_3  = oo_.irfs.in_a_obs_nua_eb;
G_U_3   = oo_.irfs.ur_a_obs_nua_eb;
G_W_3   = oo_.irfs.wr_a_obs_nua_eb;
G_RER_3 = oo_.irfs.rer_ba_obs_nua_eb;

Eu_Y_3  = oo_.irfs.dgdp2_b_t_nua_eb;
Eu_C_3  = oo_.irfs.c2_b_obs_nua_eb;
Eu_Pi_3 = oo_.irfs.pi_b_t_nua_eb;
Eu_K_3  = oo_.irfs.k_b_t_nua_eb;
Eu_B_3  = oo_.irfs.BG_b_t_nua_eb;
Eu_In_3 = oo_.irfs.in_b_obs_nua_eb;
Eu_U_3  = oo_.irfs.ur_b_obs_nua_eb;
Eu_W_3  = oo_.irfs.wr_b_obs_nua_eb;
save G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3

load Buba_Fiskal_Erweiterung_FU_CA_results;
G_Y_4   = oo_.irfs.dgdp2_a_t_nua_eb;
G_C_4   = oo_.irfs.c2_a_obs_nua_eb;
G_Pi_4  = oo_.irfs.pi_a_t_nua_eb;
G_K_4   = oo_.irfs.k_a_t_nua_eb;
G_B_4   = oo_.irfs.BG_a_t_nua_eb;
G_R_4   = oo_.irfs.i_policy_t_nua_eb;
G_In_4  = oo_.irfs.in_a_obs_nua_eb;
G_U_4   = oo_.irfs.ur_a_obs_nua_eb;
G_W_4   = oo_.irfs.wr_a_obs_nua_eb;
G_RER_4 = oo_.irfs.rer_ba_obs_nua_eb;

Eu_Y_4  = oo_.irfs.dgdp2_b_t_nua_eb;
Eu_C_4  = oo_.irfs.c2_b_obs_nua_eb;
Eu_Pi_4 = oo_.irfs.pi_b_t_nua_eb;
Eu_K_4  = oo_.irfs.k_b_t_nua_eb;
Eu_B_4  = oo_.irfs.BG_b_t_nua_eb;
Eu_In_4 = oo_.irfs.in_b_obs_nua_eb;
Eu_U_4  = oo_.irfs.ur_b_obs_nua_eb;
Eu_W_4  = oo_.irfs.wr_b_obs_nua_eb;
save G_Y_4 G_C_4 G_Pi_4 G_K_4 G_B_4 G_R_4 G_In_4 G_U_4 G_W_4 G_RER_4 Eu_Y_4 Eu_C_4 Eu_Pi_4 Eu_K_4 Eu_B_4 Eu_In_4 Eu_U_4 Eu_W_4


load G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2
load G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1
load G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3
load G_Y_4 G_C_4 G_Pi_4 G_K_4 G_B_4 G_R_4 G_In_4 G_U_4 G_W_4 G_RER_4 Eu_Y_4 Eu_C_4 Eu_Pi_4 Eu_K_4 Eu_B_4 Eu_In_4 Eu_U_4 Eu_W_4
 
% 
T=size(G_Y_1,2);
t=1:T;
Zeroline=zeros(size(t,2),1);
f = figure;
set(f,'InvertHardcopy','off','Color',[1 1 1],'Position',[1 1 800 1000]);
subplot (4,4,1);
   plot(t,G_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
 subplot (4,4,2); 
   plot(t,G_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_C_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_C_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,3);
   plot(t,G_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_In_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_In_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,4);
   plot(t,G_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('GE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,5);
   plot(t,G_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_U_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_U_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('GE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,6);
   plot(t,G_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_W_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_W_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on    
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,7); 
   plot(t,G_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_B_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_B_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,8);
   plot(t,Eu_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_B_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_B_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,9); 
  plot(t,Eu_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,10); 
   plot(t,Eu_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_C_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_C_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,11); 
   plot(t,Eu_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_In_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_In_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,12); 
   plot(t,Eu_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,13); 
   plot(t,Eu_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_U_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_U_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,14); 
   plot(t,Eu_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_W_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_W_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,15);
   plot(t,G_RER_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_RER_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_RER_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE ToT to Germany','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,16);
   plot(t,G_R_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_R_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_R_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('Interest Rate','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
   h_legend=legend('SQ','RE','TH','CA','Orientation','horizontal','Location','SouthOutside');
%    legend('boxoff');
    pos = get(h_legend,'Position');
    pos(1) = 0.253;
    pos(2) = 0.038;
    pos(3) = 0.53;
    pos(4) = 0.033;
    set(h_legend,'Position',pos);
    set(h_legend,'FontSize',13,'FontWeight','Bold');



load Buba_Fiskal_Erweiterung_FU_baseline_results;
G_Y_1   = oo_.irfs.dgdp2_a_t_nub_eb;
G_C_1   = oo_.irfs.c2_a_obs_nub_eb;
G_Pi_1  = oo_.irfs.pi_a_t_nub_eb;
G_K_1   = oo_.irfs.k_a_t_nub_eb;
G_B_1   = oo_.irfs.BG_a_t_nub_eb;
G_R_1   = oo_.irfs.i_policy_t_nub_eb;
G_In_1  = oo_.irfs.in_a_obs_nub_eb;
G_U_1   = oo_.irfs.ur_a_obs_nub_eb;
G_W_1   = oo_.irfs.wr_a_obs_nub_eb;
G_RER_1 = oo_.irfs.rer_ba_obs_nub_eb;

Eu_Y_1  = oo_.irfs.dgdp2_b_t_nub_eb;
Eu_C_1  = oo_.irfs.c2_b_obs_nub_eb;
Eu_Pi_1 = oo_.irfs.pi_b_t_nub_eb;
Eu_K_1  = oo_.irfs.k_b_t_nub_eb;
Eu_B_1  = oo_.irfs.BG_b_t_nub_eb;
Eu_In_1 = oo_.irfs.in_b_obs_nub_eb;
Eu_U_1  = oo_.irfs.ur_b_obs_nub_eb;
Eu_W_1  = oo_.irfs.wr_b_obs_nub_eb;
save G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1


load Buba_Fiskal_Erweiterung_FU_Eq_results;
G_Y_2   = oo_.irfs.dgdp2_a_t_nub_eb;
G_C_2   = oo_.irfs.c2_a_obs_nub_eb;
G_Pi_2  = oo_.irfs.pi_a_t_nub_eb;
G_K_2   = oo_.irfs.k_a_t_nub_eb;
G_B_2   = oo_.irfs.BG_a_t_nub_eb;
G_R_2   = oo_.irfs.i_policy_t_nub_eb;
G_In_2  = oo_.irfs.in_a_obs_nub_eb;
G_U_2   = oo_.irfs.ur_a_obs_nub_eb;
G_W_2   = oo_.irfs.wr_a_obs_nub_eb;
G_RER_2 = oo_.irfs.rer_ba_obs_nub_eb;


Eu_Y_2  = oo_.irfs.dgdp2_b_t_nub_eb;
Eu_C_2  = oo_.irfs.c2_b_obs_nub_eb;
Eu_Pi_2 = oo_.irfs.pi_b_t_nub_eb;
Eu_K_2  = oo_.irfs.k_b_t_nub_eb;
Eu_B_2  = oo_.irfs.BG_b_t_nub_eb;
Eu_In_2 = oo_.irfs.in_b_obs_nub_eb;
Eu_U_2  = oo_.irfs.ur_b_obs_nub_eb;
Eu_W_2  = oo_.irfs.wr_b_obs_nub_eb;

save G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2


load Buba_Fiskal_Erweiterung_FU_TH_results;
G_Y_3   = oo_.irfs.dgdp2_a_t_nub_eb;
G_C_3   = oo_.irfs.c2_a_obs_nub_eb;
G_Pi_3  = oo_.irfs.pi_a_t_nub_eb;
G_K_3   = oo_.irfs.k_a_t_nub_eb;
G_B_3   = oo_.irfs.BG_a_t_nub_eb;
G_R_3   = oo_.irfs.i_policy_t_nub_eb;
G_In_3  = oo_.irfs.in_a_obs_nub_eb;
G_U_3   = oo_.irfs.ur_a_obs_nub_eb;
G_W_3   = oo_.irfs.wr_a_obs_nub_eb;
G_RER_3 = oo_.irfs.rer_ba_obs_nub_eb;

Eu_Y_3  = oo_.irfs.dgdp2_b_t_nub_eb;
Eu_C_3  = oo_.irfs.c2_b_obs_nub_eb;
Eu_Pi_3 = oo_.irfs.pi_b_t_nub_eb;
Eu_K_3  = oo_.irfs.k_b_t_nub_eb;
Eu_B_3  = oo_.irfs.BG_b_t_nub_eb;
Eu_In_3 = oo_.irfs.in_b_obs_nub_eb;
Eu_U_3  = oo_.irfs.ur_b_obs_nub_eb;
Eu_W_3  = oo_.irfs.wr_b_obs_nub_eb;
save G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3

load Buba_Fiskal_Erweiterung_FU_CA_results;
G_Y_4   = oo_.irfs.dgdp2_a_t_nub_eb;
G_C_4   = oo_.irfs.c2_a_obs_nub_eb;
G_Pi_4  = oo_.irfs.pi_a_t_nub_eb;
G_K_4   = oo_.irfs.k_a_t_nub_eb;
G_B_4   = oo_.irfs.BG_a_t_nub_eb;
G_R_4   = oo_.irfs.i_policy_t_nub_eb;
G_In_4  = oo_.irfs.in_a_obs_nub_eb;
G_U_4   = oo_.irfs.ur_a_obs_nub_eb;
G_W_4   = oo_.irfs.wr_a_obs_nub_eb;
G_RER_4 = oo_.irfs.rer_ba_obs_nub_eb;

Eu_Y_4  = oo_.irfs.dgdp2_b_t_nub_eb;
Eu_C_4  = oo_.irfs.c2_b_obs_nub_eb;
Eu_Pi_4 = oo_.irfs.pi_b_t_nub_eb;
Eu_K_4  = oo_.irfs.k_b_t_nub_eb;
Eu_B_4  = oo_.irfs.BG_b_t_nub_eb;
Eu_In_4 = oo_.irfs.in_b_obs_nub_eb;
Eu_U_4  = oo_.irfs.ur_b_obs_nub_eb;
Eu_W_4  = oo_.irfs.wr_b_obs_nub_eb;
save G_Y_4 G_C_4 G_Pi_4 G_K_4 G_B_4 G_R_4 G_In_4 G_U_4 G_W_4 G_RER_4 Eu_Y_4 Eu_C_4 Eu_Pi_4 Eu_K_4 Eu_B_4 Eu_In_4 Eu_U_4 Eu_W_4


load G_Y_2 G_C_2 G_Pi_2 G_K_2 G_B_2 G_R_2 G_In_2 G_U_2 G_W_2 G_RER_2 Eu_Y_2 Eu_C_2 Eu_Pi_2 Eu_K_2 Eu_B_2 Eu_In_2 Eu_U_2 Eu_W_2
load G_Y_1 G_C_1 G_Pi_1 G_K_1 G_B_1 G_R_1 G_In_1 G_U_1 G_W_1 G_RER_1 Eu_Y_1 Eu_C_1 Eu_Pi_1 Eu_K_1 Eu_B_1 Eu_In_1 Eu_U_1 Eu_W_1
load G_Y_3 G_C_3 G_Pi_3 G_K_3 G_B_3 G_R_3 G_In_3 G_U_3 G_W_3 G_RER_3 Eu_Y_3 Eu_C_3 Eu_Pi_3 Eu_K_3 Eu_B_3 Eu_In_3 Eu_U_3 Eu_W_3
load G_Y_4 G_C_4 G_Pi_4 G_K_4 G_B_4 G_R_4 G_In_4 G_U_4 G_W_4 G_RER_4 Eu_Y_4 Eu_C_4 Eu_Pi_4 Eu_K_4 Eu_B_4 Eu_In_4 Eu_U_4 Eu_W_4
 
% 
T=size(G_Y_1,2);
t=1:T;
Zeroline=zeros(size(t,2),1);
f = figure;
set(f,'InvertHardcopy','off','Color',[1 1 1],'Position',[1 1 800 1000]);
subplot (4,4,1);
   plot(t,G_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
 subplot (4,4,2); 
   plot(t,G_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_C_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_C_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,3);
   plot(t,G_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_In_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_In_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,4);
   plot(t,G_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('GE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,5);
   plot(t,G_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_U_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_U_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('GE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,6);
   plot(t,G_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_W_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_W_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on    
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,7); 
   plot(t,G_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_B_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_B_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('GE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,8);
   plot(t,Eu_B_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_B_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_B_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_B_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Debt','FontSize',11,'FontWeight','Bold');
   xlabel ('time');   
subplot (4,4,9); 
  plot(t,Eu_Y_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Y_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_Y_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_Y_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Output','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,10); 
   plot(t,Eu_C_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_C_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_C_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_C_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Consumption','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,11); 
   plot(t,Eu_In_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_In_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_In_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_In_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Investment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,12); 
   plot(t,Eu_Pi_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_Pi_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_Pi_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_Pi_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Inflation','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,13); 
   plot(t,Eu_U_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_U_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_U_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_U_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Unemployment','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,14); 
   plot(t,Eu_W_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,Eu_W_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,Eu_W_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,Eu_W_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE Real Wage','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,15);
   plot(t,G_RER_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_RER_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_RER_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_RER_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('RoE ToT to Germany','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
subplot (4,4,16);
   plot(t,G_R_1,'b-+','LineWidth',1.5,'Markersize',2);
   hold on
   plot(t,G_R_2,'r-*','LineWidth',1.5,'Markersize',2);
   hold on
    plot(t,G_R_3,'g-+','LineWidth',1.5,'Markersize',2);
    hold on
    plot(t,G_R_4,'m-+','LineWidth',1.5,'Markersize',2);
    hold on
   plot(t,Zeroline,'k','LineWidth',1);
   title('Interest Rate','FontSize',11,'FontWeight','Bold');
   xlabel ('time');
   h_legend=legend('SQ','RE','TH','CA','Orientation','horizontal','Location','SouthOutside');
%    legend('boxoff');
    pos = get(h_legend,'Position');
    pos(1) = 0.253;
    pos(2) = 0.038;
    pos(3) = 0.53;
    pos(4) = 0.033;
    set(h_legend,'Position',pos);
    set(h_legend,'FontSize',13,'FontWeight','Bold');
