%excute model_zanetti
clear all;
clc;


deltan_ts= 0.105;
delta_k= 0.025;
beta= 0.990;
sigma=1;
ns=0.65;%estimated
qs=0.92;%estimated
eta=0.95;   %estimated
mu=0.75;  % estimated
theta=0.28;   %estimated
h= 0.93;   %estimated
phi=1.99;  % estimated
kappa=6.01;   % estimated
chir_share=3.27;% estimated
vartheta=0.28; %estimated
zss=1;
Vs=1;
ass=1;
dss=1;
nuss=1;

%shocks estimated value
rho_z=0.91;
rho_a=0.95;
rho_nu=0.95;
rho_d=0.93;
rho_am=0.72;
rho_deltaN=0.83;

sigma_z=0.88;
sigma_a=0.47;
sigma_nu=0.48;
sigma_d=0.48;
sigma_am=0.22;
sigma_deltaN=0.20;
sdnewsam4=0.001;%
sdnewsam8=0.001;%
 sdnewsd4=4;
 sdnewsd8=0.24;
 sdnewsz4=0.28;
 sdnewsz8 =0.45;
sdnewsnu4=4;
 sdnewsnu8=4;
 sdnewsa4=0.72;
 sdnewsa8=0.25;
 sdnewsdeltaN4=0.001;%
 sdnewsdeltaN8= 0.001;%

alpha=eta/(1-eta);
qks=(1/beta)-(1-delta_k);
gamma1=qks;
gamma2=(vartheta/(1-vartheta))*qks;
us=1-ns;
ks=((theta/qks)^(1/(1-theta)))*ns;
ys=(ns^(1-theta))*(ks^theta);
kps=ks;
is=(1-(1-delta_k))*kps;
ms=deltan_ts*ns;
amss=(ms/(us^mu))*((qs/ms)^(1-mu));
vs=(ms/(amss*us^mu))^(1/(1-mu));
xs=vs/us;
B=(chir_share/100)*ys/((xs^alpha)*vs);
gs=B*(xs^alpha);
xis=gs/qs;
cs=ys-(gs*vs)-is;
lams=(1-h*beta)/((1-h)*cs);
ws=((1-theta)*ys/ns)-xis*(1-beta*(1-deltan_ts));
khi=(lams/((1-eta)*ns^phi))*(ws-eta*(((1-theta)*(ys/ns))+(beta*gs)));
Ss=lams;
ps=ms/us;
save parameter_main6  h beta sigma gamma1 gamma2 kappa zss delta_k theta eta phi khi mu  B alpha lams dss cs qks Vs Ss nuss is ys ass ks ns kps xis gs ws deltan_ts us  ps ms qs vs amss xs rho_a rho_nu rho_z rho_d rho_am rho_deltaN sigma_a sigma_nu sigma_z sigma_d sigma_am sigma_deltaN sdnewsam4 sdnewsam8  sdnewsd4 sdnewsd8  sdnewsz4 sdnewsz8 sdnewsnu4 sdnewsnu8 sdnewsa4 sdnewsa8  sdnewsdeltaN4 sdnewsdeltaN8;
dynare zanetti6