clear; 
clc;
% rmpath('C:/dynare/rbc'); 
% addpath('C:/dynare/rbc');

alphha = 1/3;
deltta = 0.025;
sigma = 2;
rhoz = 0.95;
varphi = 2;

betta = 0.99;

dynare rbclinear.mod

figure(1)

subplot(3,3,1)
plot(Y_eps_z,'-ko','LineWidth',1.5)
xlabel('Production Quantity','FontSize',15,'FontWeight','bold')
hold on

subplot(3,3,2)
plot(C_eps_z,'-ko','LineWidth',1.5)
xlabel('Consumption','FontSize',15,'FontWeight','bold')
hold on

subplot(3,3,3)
plot(II_eps_z,'-ko','LineWidth',1.5)
xlabel('Investment','FontSize',15,'FontWeight','bold')
hold on

subplot(3,3,4)
plot(N_eps_z,'-ko','LineWidth',1.5)
xlabel('Labor','FontSize',15,'FontWeight','bold')
hold on

subplot(3,3,5)
plot(K_eps_z,'-ko','LineWidth',1.5)
xlabel('Capital','FontSize',15,'FontWeight','bold')
hold on

subplot(3,3,6)
plot(W_eps_z,'-ko','LineWidth',1.5)
xlabel('Wage','FontSize',15,'FontWeight','bold')
hold on

subplot(3,3,7)
plot(R_eps_z,'-ko','LineWidth',1.5)
xlabel('Interest Rate','FontSize',15,'FontWeight','bold')
hold on

subplot(3,3,8)
plot(Z_eps_z,'-ko','LineWidth',1.5)
xlabel('Technology Shock','FontSize',15,'FontWeight','bold')
hold on
sgtitle('RBC Calibration Result')

betta = 0.997;

dynare rbclinear.mod

figure(1)

subplot(3,3,1)
plot(Y_eps_z,'-bx','LineWidth',1.5)

subplot(3,3,2)
plot(C_eps_z,'-bx','LineWidth',1.5)

subplot(3,3,3)
plot(II_eps_z,'-bx','LineWidth',1.5)

subplot(3,3,4)
plot(N_eps_z,'-bx','LineWidth',1.5)

subplot(3,3,5)
plot(K_eps_z,'-bx','LineWidth',1.5)

subplot(3,3,6)
plot(W_eps_z,'-bx','LineWidth',1.5)

subplot(3,3,7)
plot(R_eps_z,'-bx','LineWidth',1.5)

subplot(3,3,8)
plot(Z_eps_z,'-bx','LineWidth',1.5)
legend('$\beta$=0.99','$\beta$=0.997','Interpreter','latex')

