% Steady State values
L_dec2_sec2=50;
L_dec2_sec3=50;
L_dec3_sec2=20;
L_dec3_sec3=48;

K_dec2_sec2=20;
K_dec2_sec3=20;
K_dec3_sec2=60;
K_dec3_sec3=20;

LS_sec2=70;
LS_sec3=98;
KS_sec2=80;
KS_sec3=40;

C_dec2_sec2=40;
C_dec2_sec3=5.775;
C_dec3_sec2=45;
C_dec3_sec3=5;

C_dec2=C_dec2_sec2+C_dec2_sec3;
C_dec3=C_dec3_sec2+C_dec3_sec3;

XS_sec2=187;
XS_sec3=176.8;

M_sec3=30;

D_sec2=178.875;
D_sec3=100.875;

I_sec2=16.775;
I_sec3=40;
I=56.775;

II_sec2_secB2=7;
II_sec2_secB3=15;
II_sec3_secB2=15;
II_sec3_secB3=10;

N_dec2_sec2=55;
N_dec2_sec3=55;
N_dec3_sec2=25;
N_dec3_sec3=53;

G_sec2=55.1;
G_sec3=55.1;

sigma_sec2=1.2;
sigma_sec3=1.2;

tauL_dec2=0.25;
tauL_dec3=0.3;
tauK_dec2=0.3;
tauK_dec3=0.3;
tauIVA_sec2=0.1;
tauIVA_sec3=0.1;


%Government calculations
ZG=tauL_dec2* (L_dec2_sec2+L_dec2_sec3 ) + tauL_dec3* (L_dec3_sec2+L_dec3_sec3 ) + tauK_dec2*(K_dec2_sec2 + K_dec2_sec3)+ tauK_dec3*(K_dec3_sec2 + K_dec3_sec3)+ tauIVA_sec2*(LS_sec2 + KS_sec2)+ tauIVA_sec3*(LS_sec3 + KS_sec3);
TR_dec1=ZG-(G_sec2 + G_sec3);


%Consumers calculations
MR_dec2=(L_dec2_sec2+L_dec2_sec3)*(1-tauL_dec2)+(K_dec2_sec2+K_dec2_sec3)*(1-tauK_dec2)+TR_dec1;
MR_dec3=(L_dec3_sec2+L_dec3_sec3)*(1-tauL_dec3)+(K_dec3_sec2+K_dec3_sec3)*(1-tauK_dec3);


%Investment calculations
phi_ME=I_sec3/I;
phi_C=I_sec2/I;
PK=((1/phi_ME)^(phi_ME))*((1/phi_C)^(phi_C));

I_dec2=(MR_dec2-C_dec2)/PK;
I_dec3=(MR_dec3-C_dec3)/PK;

delta_dec2=I_dec2/(K_dec2_sec2+K_dec2_sec3);
delta_dec3=I_dec3/(K_dec3_sec2+K_dec3_sec3);


% Parameters
alpha_sec2=KS_sec2/(KS_sec2+LS_sec2);
alpha_sec3=KS_sec3/(KS_sec3+LS_sec3);

c_sec2=(KS_sec2+LS_sec2)/((KS_sec2^alpha_sec2)*(LS_sec2^(1-alpha_sec2)));
c_sec3=(KS_sec3+LS_sec3)/((KS_sec3^alpha_sec3)*(LS_sec3^(1-alpha_sec3)));

gamma_dec2_sec2=(C_dec2_sec2)/((1-tauL_dec2)*(N_dec2_sec2-L_dec2_sec2)+C_dec2_sec2);
gamma_dec2_sec3=(C_dec2_sec3)/((1-tauL_dec2)*(N_dec2_sec3-L_dec2_sec3)+C_dec2_sec3);
gamma_dec3_sec2=(C_dec3_sec2)/((1-tauL_dec3)*(N_dec3_sec2-L_dec3_sec2)+C_dec3_sec2);
gamma_dec3_sec3=(C_dec3_sec3)/((1-tauL_dec3)*(N_dec3_sec3-L_dec3_sec3)+C_dec3_sec3);

betatilde_dec2=(PK)/((1-tauK_dec2)+(1-delta_dec2)*PK);
betatilde_dec3=(PK)/((1-tauK_dec3)+(1-delta_dec3)*PK);

beta_dec2_sec2=C_dec2_sec2/C_dec2;
beta_dec2_sec3=C_dec2_sec3/C_dec2;
beta_dec3_sec2=C_dec3_sec2/C_dec3;
beta_dec3_sec3=C_dec3_sec3/C_dec3;

va_sec2=(KS_sec2+LS_sec2)/XS_sec2;
va_sec3=(KS_sec3+LS_sec3)/XS_sec3;

a_sec2_secB2=II_sec2_secB2/XS_sec2;
a_sec2_secB3=II_sec3_secB2/XS_sec2;
a_sec3_secB2=II_sec2_secB3/XS_sec3;
a_sec3_secB3=II_sec3_secB3/XS_sec3;




mu_sec2=0;
mu_sec3=(M_sec3^(1/sigma_sec3))/(M_sec3^(1/sigma_sec3)+D_sec3^(1/sigma_sec3));

b_sec2=1;
b_sec3=(D_sec3+M_sec3)/((mu_sec3*M_sec3^((sigma_sec3-1)/sigma_sec3)+(1-mu_sec3)*D_sec3^((sigma_sec3-1)/sigma_sec3))^(sigma_sec3/(sigma_sec3-1)));


save pruebaparams a_sec2_secB2 a_sec2_secB3 a_sec3_secB2 a_sec3_secB3 alpha_sec2 alpha_sec3 b_sec2 b_sec3 beta_dec2_sec2 beta_dec2_sec3 beta_dec3_sec2 beta_dec3_sec3 betatilde_dec2 betatilde_dec3 c_sec2 c_sec3 delta_dec2 delta_dec3 gamma_dec2_sec2 gamma_dec2_sec3 gamma_dec3_sec2 gamma_dec3_sec3 mu_sec2 mu_sec3 phi_C phi_ME va_sec2 va_sec3





