Skip to content
Snippets Groups Projects

Compare revisions

Changes are shown as if the source revision was being merged into the target revision. Learn more about comparing revisions.

Source

Select target project
No results found
Select Git revision
Loading items

Target

Select target project
  • giovanma/dynare
  • giorgiomas/dynare
  • Vermandel/dynare
  • Dynare/dynare
  • normann/dynare
  • MichelJuillard/dynare
  • wmutschl/dynare
  • FerhatMihoubi/dynare
  • sebastien/dynare
  • lnsongxf/dynare
  • rattoma/dynare
  • CIMERS/dynare
  • FredericKarame/dynare
  • SumuduK/dynare
  • MinjeJeon/dynare
  • camilomrch/dynare
  • DoraK/dynare
  • avtishin/dynare
  • selma/dynare
  • claudio_olguin/dynare
  • jeffjiang07/dynare
  • EthanSystem/dynare
  • stepan-a/dynare
  • wjgatt/dynare
  • JohannesPfeifer/dynare
  • gboehl/dynare
  • chskcau/dynare-doc-fixes
27 results
Select Git revision
Loading items
Show changes
Show whitespace changes
Inline Side-by-side
Showing
with 0 additions and 1992 deletions
dynare++/tests/example1.mod deleted 100644 → 0
View file @ 219d2bb3
// this is an example from the tutorial
var Y, C, K, A, H, B;
varexo EPS, NU;
parameters rho, beta, alpha, delta, theta, psi, tau;
alpha = 0.36;
rho = 0.95;
tau = 0.025;
beta = 1/(1.03^0.25);
delta = 0.025;
psi = 0;
theta = 2.95;
model;
C*theta*H^(1+psi) = (1-alpha)*Y;
beta*exp(B)*C/exp(B(1))/C(1)*
(exp(B(1))*alpha*Y(1)/K(1)+1-delta) = 1;
Y = exp(A)*K^alpha*H^(1-alpha);
K = exp(B(-1))*(Y(-1)-C(-1)) + (1-delta)*K(-1);
A = rho*A(-1) + tau*B(-1) + EPS;
B = tau*A(-1) + rho*B(-1) + NU;
end;
initval;
A = 0;
B = 0;
H = ((1-alpha)/(theta*(1-(delta*alpha)/(1/beta-1+delta))))^(1/(1+psi));
Y = (alpha/(1/beta-1+delta))^(alpha/(1-alpha))*H;
K = alpha/(1/beta-1+delta)*Y;
C = Y - delta*K;
end;
vcov = [
0.0002 0.00005;
0.00005 0.0001
];
order = 7;
dynare++/tests/example1_optim.mod deleted 100644 → 0
View file @ 219d2bb3
// this is a file trying to replicate example1.mod as optimization of a social planner
// it serves also as an example combining +2 lead and optimal policy
var Y, C, K, A, H, B;
varexo EPS, NU;
parameters rho, beta, alpha, delta, theta, psi, tau;
alpha = 0.36;
rho = 0.95;
tau = 0.025;
beta = 1/(1.03^0.25);
delta = 0.025;
psi = 0;
theta = 2.95;
planner_objective log(C)-theta*H^(1+psi)/(1+psi);
planner_discount beta;
model;
//Y = exp(A)*K^alpha*H^(1-alpha);
Y = exp(A)*exp(A(+1))*exp(A(+2))*K^alpha*H^(1-alpha);
K = exp(B(-1))*(Y(-1)-C(-1)) + (1-delta)*K(-1);
A = rho*A(-1) + tau*B(-1) + EPS;
B = tau*A(-1) + rho*B(-1) + NU;
end;
initval;
A = 0;
B = 0;
H = ((1-alpha)/(theta*(1-(delta*alpha)/(1/beta-1+delta))))^(1/(1+psi));
Y = (alpha/(1/beta-1+delta))^(alpha/(1-alpha))*H;
K = alpha/(1/beta-1+delta)*Y;
C = Y - delta*K;
end;
vcov = [
0.0002 0.00005;
0.00005 0.0001
];
order = 2;
dynare++/tests/gentay1a.dyn deleted 100644 → 0
View file @ 219d2bb3
var AF AH BIGGAMF BIGGAMH BIGGAMIMPF BIGGAMIMPH BIGGAMMF BIGGAMMH BIGGAMNF BIGGAMNH BIGGAMQF BIGGAMQH BIGGAM_MONF BIGGAM_MONH BIGGAM_MOTF BIGGAM_MOTH BIGGAM_O_NF BIGGAM_O_NH BIGGAM_O_TF BIGGAM_O_TH CAPAF CAPAH CF CF_NAT CH CH_NAT CURBALF_RAT CURBALH_RAT DEEF DEEH DEPEX EXPORTSF EXPORTSF_NAT EXPORTSH EXPORTSH_NAT EYEF EYEH GAF GAH GAMMAF GAMMAH GDPF GDPF_NAT GDPGAPF GDPGAPH GDPH GDPH_NAT GF_NAT GH_NAT GNF GNH HF HH IMPORTSF IMPORTSF_NAT IMPORTSH IMPORTSH_NAT KF KH KNF KNF_RAT KNH KNH_RAT KTF KTF_RAT KTH KTH_RAT K_OF K_OF_RAT K_OH K_OH_RAT LANDF LANDH LF LH LNF LNH LTF LTH L_OF L_OH MARGUTF MARGUTH MF MF_NAT MH MH_NAT M_ONF M_ONH M_OTF M_OTH NF NH NNF NNH O_NF O_NH O_TF O_TH PIE4F PIE4H PIE4TARF PIE4TARH PIEBARMF PIEBARMH PIEBARQF PIEBARQH PIEF PIEH PIENF PIENH PIEWF PIEWH PSIF PSIH PSIPRIMEF PSIPRIMEH QF QH Q_ONF Q_ONH Q_OTF Q_OTH REALBF REALBH REALEX REALEXF REALEXH REALFINH REALMCNF REALMCNH REALMCTF REALMCTH REALMONEYF REALMONEYH REALPBARMF REALPBARMH REALPBARQF REALPBARQH REALPLANDF REALPLANDH REALPMF REALPMH REALPNF REALPNH REALPQF REALPQH REALPXF REALPXH REALP_MOF REALP_MOH REALP_ONF REALP_ONH REALP_OTF REALP_OTH REALP_QOF REALP_QOH REALRF REALRH REALTBALF REALTBALH REALWF REALWH RNOMF RNOMH SHOPF SHOPH SHOPPRIMEF SHOPPRIMEH TF TH T_OF T_OH VELOF VELOH VPRIMEF VPRIMEH XF_NAT XH_NAT XI ZBF ZBH ZEYEF ZEYEH ZNF ZNH ZTF ZTH ZUF ZUH ZZ_CF ZZ_CH ZZ_CURBALF_RAT ZZ_CURBALH_RAT ZZ_C_SHARF ZZ_C_SHARH ZZ_DPIE4F ZZ_DPIE4H ZZ_DRNOMF ZZ_DRNOMH ZZ_EXPORTSF ZZ_EXPORTSH ZZ_EYEF ZZ_EYEH ZZ_EYE_SHARF ZZ_EYE_SHARH ZZ_GDPF ZZ_GDPH ZZ_GDP_PPP_RATH ZZ_GF ZZ_GH ZZ_IMPORTSF ZZ_IMPORTSH ZZ_M_ON_SHARF ZZ_M_ON_SHARH ZZ_M_OT_SHARF ZZ_M_OT_SHARH ZZ_M_O_SHARF ZZ_M_O_SHARH ZZ_M_SHARF ZZ_M_SHARH ZZ_N_SHARF ZZ_N_SHARH ZZ_PIE4F ZZ_PIE4H ZZ_REALEX ZZ_RNOMF ZZ_RNOMH ZZ_UTILITYF ZZ_UTILITYH ZZ_XBALF_TOT_RAT ZZ_XBALH_TOT_RAT Z_OF Z_OH PIEF1 PIEF2 PIEH1 PIEH2 REALBH1;
varexo E_ZBH E_ZUH E_ZUF E_ZEYEH E_ZEYEF E_GAMMAH E_GAMMAF E_LANDH E_LANDF E_GAH E_GAF E_CAPAH E_CAPAF;
parameters COSTLF COSTLH EPSF EPSH EPSQMF EPSQMH GLAMBDAF GLAMBDAH SIGMAF SIGMAH SSF SSH XR1F XR1H XR2F XR2H XR3F XR3H XR4F XR4H ALPHANF ALPHANH ALPHATF ALPHATH ALPHA_OF ALPHA_OH A_ONEF A_ONEH A_ZEROF A_ZEROH B0F B0H B1F B1H B2F B2H BET CAPAF_SS CAPAH_SS CHI0 CHI1 CHI2 CHI3 COSTF COSTH COST_MONF COST_MONH COST_MOTF COST_MOTH COST_O_NF COST_O_NH COST_O_TF COST_O_TH DELTAF DELTAH EPS_ONF EPS_ONH EPS_OTF EPS_OTH ETAF ETAH E_PIE4TARF E_PIE4TARH E_ZBF GAMA_NF GAMA_NH GAMA_TF GAMA_TH GAMMAF_SS GAMMAH_SS GAMMA_LANDF GAMMA_LANDH GA_RATF GA_RATH GDPF_EXOG GDPH_EXOG GN_RATF GN_RATH LANDF_SS LANDH_SS LILCAPPA1F LILCAPPA1H LILCAPPA2F LILCAPPA2H LILCAPPA3F LILCAPPA3H LILCAPPA4F LILCAPPA4H LILCAPPA5F LILCAPPA5H LILCAPPA6F LILCAPPA6H LILCAPPA7F LILCAPPA7H LILCAPPA8F LILCAPPA8H NYF NYH NY_NF NY_NH NY_TF NY_TH OMEGA0F OMEGA0H OMEGAF OMEGAH PHIF PHIH PIE4TARF_SS PIE4TARH_SS REALEX_EXOG REALPBARMF_EXOG REALPBARMH_EXOG REALPNF_EXOG REALPNH_EXOG REALP_MOF_EXOG REALP_MOH_EXOG RNOM_EXOGF RNOM_EXOGH THETAF THETAH XDUMF XDUMH XIXI_NF XIXI_NH XIXI_OF XIXI_OH XIXI_TF XIXI_TH ZBF_SS ZBH_SS ZEDF ZEDH ZEYEF_SS ZEYEH_SS ZNF_SS ZNH_SS ZTF_SS ZTH_SS ZUF_SS ZUH_SS Z_OF_SS Z_OH_SS;
COSTLF=5;
COSTLH=5;
EPSF=1.1;
EPSH=1.1;
EPSQMF=4;
EPSQMH=1.1;
GLAMBDAF=0;
GLAMBDAH=0;
SIGMAF=0.333333333333333;
SIGMAH=0.333333333333333;
SSF=0.95;
SSH=0.05;
XR1F=1;
XR1H=0.40;
XR2F=0;
XR2H=0.1;
XR3F=0.5;
XR3H=0.84;
XR4F=0;
XR4H=0;
ALPHANF=0.33;
ALPHANH=0.4;
ALPHATF=0.33;
ALPHATH=0.4;
ALPHA_OF=0.2;
ALPHA_OH=0.2;
A_ONEF=0.075;
A_ONEH=0.075;
A_ZEROF=0.011;
A_ZEROH=0.011;
B0F=1;
B0H=1;
B1F=0.95;
B1H=0.95;
B2F=1;
B2H=1;
BET=0.99263753614514;
CAPAF_SS=11;
CAPAH_SS=11;
CHI0=0.05;
CHI1=.1;
CHI2=0;
CHI3=0;
COSTF=5;
COSTH=5;
COST_MONF=5;
COST_MONH=5;
COST_MOTF=5;
COST_MOTH=5;
COST_O_NF=5;
COST_O_NH=5;
COST_O_TF=5;
COST_O_TH=5;
DELTAF=0.025;
DELTAH=0.025;
EPS_ONF=3;
EPS_ONH=3;
EPS_OTF=3;
EPS_OTH=3;
ETAF=0.35;
ETAH=0.2;
E_PIE4TARF=0;
E_PIE4TARH=0;
E_ZBF=0;
GAMA_NF=0.3;
GAMA_NH=0.1;
GAMA_TF=0.3;
GAMA_TH=0.5;
GAMMAF_SS=0.5;
GAMMAH_SS=0.25;
GAMMA_LANDF=0.1;
GAMMA_LANDH=0.1;
GA_RATF=0.05;
GA_RATH=0.05;
GDPF_EXOG=2.41500497257461;
GDPH_EXOG=2.90512477822209;
GN_RATF=0.1;
GN_RATH=0.1;
LANDF_SS=0.1;
LANDH_SS=0.1;
LILCAPPA1F=0;
LILCAPPA1H=0;
LILCAPPA2F=800;
LILCAPPA2H=400;
LILCAPPA3F=0;
LILCAPPA3H=0;
LILCAPPA4F=800;
LILCAPPA4H=400;
LILCAPPA5F=0;
LILCAPPA5H=0;
LILCAPPA6F=800;
LILCAPPA6H=400;
LILCAPPA7F=0;
LILCAPPA7H=0;
LILCAPPA8F=0;
LILCAPPA8H=0;
NYF=0.98;
NYH=0.0499999999999999;
NY_NF=0.98;
NY_NH=0.0499999999999999;
NY_TF=0.98;
NY_TH=0.02;
OMEGA0F=60;
OMEGA0H=60;
OMEGAF=5;
OMEGAH=5;
PHIF=3;
PHIH=3;
PIE4TARF_SS=1.125;
PIE4TARH_SS=1.125;
REALEX_EXOG=1.3734519289908;
REALPBARMF_EXOG=0.87146958398196;
REALPBARMH_EXOG=1.19072687148694;
REALPNF_EXOG=0.840675522925242;
REALPNH_EXOG=0.902486321747893;
REALP_MOF_EXOG=0.966533486000563;
REALP_MOH_EXOG=1.63690883121281;
RNOM_EXOGF=1.00741707177773;
RNOM_EXOGH=1.00741707177773;
THETAF=6;
THETAH=6;
XDUMF=1;
XDUMH=1;
XIXI_NF=0.75;
XIXI_NH=0.75;
XIXI_OF=0.75;
XIXI_OH=0.75;
XIXI_TF=0.75;
XIXI_TH=0.75;
ZBF_SS=0;
ZBH_SS=0;
ZEDF=2.5;
ZEDH=2.5;
ZEYEF_SS=0;
ZEYEH_SS=0;
ZNF_SS=1;
ZNH_SS=1;
ZTF_SS=1;
ZTH_SS=0.6;
ZUF_SS=1;
ZUH_SS=1;
Z_OF_SS=1;
Z_OH_SS=1;
model;
ZZ_UTILITYH = (ZUH*(CH-HH)^(1-SIGMAH)-1)/(1-SIGMAH)-CAPAH*LH^(1+ZEDH)/(1+ZEDH) ;
ZZ_GH = 100*log(GH_NAT) ;
ZZ_CURBALH_RAT = CURBALH_RAT*100 ;
ZZ_M_SHARH = REALPBARMH*MH/AH ;
ZZ_M_O_SHARH = (REALP_MOH*M_ONH+REALP_MOH*M_OTH)/AH ;
ZZ_M_ON_SHARH = REALP_MOH*M_ONH/AH ;
ZZ_M_OT_SHARH = REALP_MOH*M_OTH/AH ;
ZZ_N_SHARH = NH*REALPNH/AH ;
ZZ_EYE_SHARH = EYEH/GDPH ;
ZZ_C_SHARH = CH/GDPH ;
ZZ_GDPH = 100*log(GDPH_NAT) ;
ZZ_CH = 100*log(CH_NAT) ;
ZZ_EYEH = 100*log(EYEH) ;
ZZ_EXPORTSH = 100*log(EXPORTSH_NAT) ;
ZZ_IMPORTSH = 100*log(IMPORTSH_NAT) ;
ZZ_XBALH_TOT_RAT = 100*(EXPORTSH_NAT-IMPORTSH_NAT)/GDPH_NAT ;
ZZ_PIE4H = 100*(PIE4H-1) ;
ZZ_DPIE4H = ZZ_PIE4H-ZZ_PIE4H(-1) ;
ZZ_RNOMH = 100*(RNOMH^4-1) ;
ZZ_DRNOMH = ZZ_RNOMH-ZZ_RNOMH(-1) ;
100*(PIE4TARH-1) = 1*100*(PIE4TARH_SS-1)+(1-1)*100*(PIE4TARH(-1)-1)+E_PIE4TARH ;
log(ZUH) = 0.3*log(ZUH_SS)+0.7*log(ZUH(-1))+E_ZUH ;
ZBH = 0.3*ZBH_SS+0.7*ZBH(-1)+E_ZBH ;
log(LANDH) = 0.05*log(LANDH_SS)+0.95*log(LANDH(-1))+E_LANDH ;
log(ZTH) = 0.05*log(ZTH_SS)+0.95*log(ZTH(-1))+E_LANDH ;
log(ZNH) = 0.05*log(ZNH_SS)+0.95*log(ZNH(-1))+E_LANDH ;
log(Z_OH) = 0.05*log(Z_OH_SS)+0.95*log(Z_OH(-1))+E_LANDH ;
ZEYEH = 0.05*ZEYEH_SS+0.95*ZEYEH(-1)+E_ZEYEH ;
CAPAH = 0.05*CAPAH_SS+0.95*CAPAH(-1)+E_CAPAH ;
log(GAMMAH) = 0.05*log(GAMMAH_SS)+0.95*log(GAMMAH(-1))+E_GAMMAH ;
BIGGAM_O_NH = 1-COST_O_NH/2*(O_NH/NH/(O_NH(-1)/NH(-1))-1)^2 ;
BIGGAM_O_TH = 1-COST_O_TH/2*(O_TH/TH/(O_TH(-1)/TH(-1))-1)^2 ;
O_NH = GAMA_NH*NH/ZNH*(REALP_ONH/(REALMCNH*ZNH))^(-XIXI_NH)*(BIGGAM_O_NH-COST_O_NH*(O_NH/NH/(O_NH(-1)/NH(-1))-1)*O_NH/NH/(O_NH(-1)/NH(-1)))^XIXI_NH/BIGGAM_O_NH ;
O_TH = GAMA_TH*TH/ZTH*(REALP_OTH/(REALMCTH*ZTH))^(-XIXI_TH)*(BIGGAM_O_TH-COST_O_TH*(O_TH/TH/(O_TH(-1)/TH(-1))-1)*O_TH/TH/(O_TH(-1)/TH(-1)))^XIXI_NH/BIGGAM_O_TH ;
NH = ZNH*((1-ALPHANH-GAMA_NH)^(1/XIXI_NH)*LNH^(1-1/XIXI_NH)+ALPHANH^(1/XIXI_NH)*KNH^(1-1/XIXI_NH)+GAMA_NH^(1/XIXI_NH)*(BIGGAM_O_NH*O_NH)^(1-1/XIXI_NH))^(XIXI_NH/(XIXI_NH-1)) ;
TH = ZTH*((1-ALPHATH-GAMA_TH)^(1/XIXI_TH)*LTH^(1-1/XIXI_TH)+ALPHATH^(1/XIXI_TH)*KTH^(1-1/XIXI_TH)+GAMA_TH^(1/XIXI_TH)*(BIGGAM_O_TH*O_TH)^(1-1/XIXI_TH))^(XIXI_TH/(XIXI_TH-1)) ;
REALMCNH = 1/ZNH*((1-ALPHANH-GAMA_NH)*REALWH^(1-XIXI_NH)+ALPHANH*REALRH^(1-XIXI_NH)+GAMA_NH*REALP_ONH^(1-XIXI_NH)*(BIGGAM_O_NH-COST_O_NH*(O_NH/NH/(O_NH(-1)/NH(-1))-1)*O_NH/NH/(O_NH(-1)/NH(-1)))^(XIXI_NH-1))^(1/(1-XIXI_NH)) ;
REALMCTH = 1/ZTH*((1-ALPHATH-GAMA_TH)*REALWH^(1-XIXI_TH)+ALPHATH*REALRH^(1-XIXI_TH)+GAMA_TH*REALP_OTH^(1-XIXI_TH)*(BIGGAM_O_TH-COST_O_TH*(O_TH/TH/(O_TH(-1)/TH(-1))-1)*O_TH/TH/(O_TH(-1)/TH(-1)))^(XIXI_TH-1))^(1/(1-XIXI_TH)) ;
MARGUTH = (CH-B2H*HH)^(-SIGMAH)*ZUH ;
HH = (1-B0H)*HH(-1)+B0H*B1H*CH(-1) ;
VPRIMEH = CAPAH*LH^ZEDH ;
AH^(1-1/EPSH) = GAMMAH^(1/EPSH)*(NYH^(1/EPSQMH)*QH^(1-1/EPSQMH)+(1-NYH)^(1/EPSQMH)*(MH*BIGGAMIMPH)^(1-1/EPSQMH))^(EPSQMH/(EPSQMH-1)*(1-1/EPSH))+(1-GAMMAH)^(1/EPSH)*NNH^(1-1/EPSH) ;
QH = GAMMAH*NYH*REALPQH^(-EPSQMH)*AH*REALPXH^(EPSQMH-EPSH) ;
MH = GAMMAH*(1-NYH)*REALPMH^(-EPSQMH)*AH*REALPXH^(EPSQMH-EPSH)*1/BIGGAMIMPH*(BIGGAMIMPH-COSTH*(MH/AH/(MH(-1)/AH(-1))-1)*MH/AH/(MH(-1)/AH(-1)))^EPSQMH ;
REALPXH = (NYH*REALPQH^(1-EPSQMH)+(1-NYH)*REALPMH^(1-EPSQMH)*(BIGGAMIMPH-COSTH*(MH/AH/(MH(-1)/AH(-1))-1)*MH/AH/(MH(-1)/AH(-1)))^(EPSQMH-1))^(1/(1-EPSQMH)) ;
BIGGAMIMPH = 1-COSTH/2*(MH/AH/(MH(-1)/AH(-1))-1)^2 ;
NNH = (1-GAMMAH)*REALPNH^(-EPSH)*AH ;
NH = NNH+ETAH*MH+ETAH*QH+GNH ;
PIENH = REALPNH/REALPNH(-1)*PIEH ;
BIGGAMNH = LILCAPPA3H/2*(PIENH/PIE4TARH^0.25-1)^2+LILCAPPA4H/2*(PIENH/PIENH(-1)-1)^2 ;
-(1-BIGGAMNH)*(REALPNH*(1-THETAH)+THETAH*REALMCNH) = -(REALPNH-REALMCNH)*(LILCAPPA3H*PIENH/PIE4TARH^0.25*(PIENH/PIE4TARH^0.25-1)+LILCAPPA4H*PIENH/PIENH(-1)*(PIENH/PIENH(-1)-1))+DEEH*PIEH(+1)*(REALPNH(+1)-REALMCNH(+1))*NH(+1)/NH*(LILCAPPA3H*PIENH(+1)/PIE4TARH^0.25*(PIENH(+1)/PIE4TARH^0.25-1)+LILCAPPA4H*PIENH(+1)/PIENH*(PIENH(+1)/PIENH-1)) ;
PIEBARQH = PIEH*REALPBARQH/REALPBARQH(-1) ;
BIGGAMQH = LILCAPPA5H/2*(PIEBARQH/PIE4TARH^0.25-1)^2+LILCAPPA6H/2*(PIEBARQH/PIEBARQH(-1)-1)^2 ;
-(1-BIGGAMQH)*REALPBARQH/REALPQH*(REALPBARQH*(1-THETAH)+ETAH*REALPNH+THETAH*REALMCTH) = -(REALPBARQH-REALMCTH)*(LILCAPPA5H*PIEBARQH/PIE4TARH^0.25*(PIEBARQH/PIE4TARH^0.25-1)+LILCAPPA6H*PIEBARQH/PIEBARQH(-1)*(PIEBARQH/PIEBARQH(-1)-1))+DEEH*PIEH(+1)*(REALPBARQH(+1)-REALMCTH(+1))*QH(+1)/QH*(LILCAPPA5H*PIEBARQH(+1)/PIE4TARH^0.25*(PIEBARQH(+1)/PIE4TARH^0.25-1)+LILCAPPA6H*PIEBARQH(+1)/PIEBARQH*(PIEBARQH(+1)/PIEBARQH-1)) ;
REALPQH = REALPBARQH+ETAH*REALPNH ;
KH = KH(-1)*(1-DELTAH)+PSIH(-1)*KH(-1) ;
PSIH = EYEH/KH-OMEGAH/2*(EYEH/KH-DELTAH*(1+ZEYEH))^2-OMEGA0H/2*(EYEH/KH-EYEH(-1)/KH(-1))^2 ;
PSIPRIMEH = 1-OMEGAH*(EYEH/KH-DELTAH*(1+ZEYEH))-OMEGA0H*(EYEH/KH-EYEH(-1)/KH(-1)) ;
1/PSIPRIMEH = DEEH*PIEH(+1)*(REALRH(+1)+1/PSIPRIMEH(+1)*(1-DELTAH+PSIH(+1)*(1-PSIPRIMEH(+1)*EYEH(+1)/(PSIH(+1)*KH(+1))))) ;
BIGGAMH = LILCAPPA1H/2*(PIEWH/PIE4TARH^0.25-1)^2+LILCAPPA2H/2*(PIEWH/PIEWH(-1)-1)^2 ;
PIEH*REALWH/REALWH(-1) = PIEWH ;
REALWH = PHIH*VPRIMEH/MARGUTH*((PHIH-1)*(1-BIGGAMH)+PIEWH*LILCAPPA1H/PIE4TARH^0.25*(PIEWH/PIE4TARH^0.25-1)+PIEWH/PIEWH(-1)*LILCAPPA2H*(PIEWH/PIEWH(-1)-1)-DEEH*PIEWH(+1)*LH(+1)/LH*LILCAPPA1H*PIEWH(+1)/PIE4TARH^0.25*(PIEWH(+1)/PIE4TARH^0.25-1)-DEEH*PIEWH(+1)*LH(+1)/LH*LILCAPPA2H*PIEWH(+1)/(REALWH/REALWH(-1))*(PIEWH(+1)/PIEWH-1))^(-1) ;
DEEH = BET*MARGUTH(+1)/PIEH(+1)/MARGUTH*(1+SHOPH+SHOPPRIMEH*VELOH)/(1+SHOPH(+1)+SHOPPRIMEH(+1)*VELOH(+1)) ;
SHOPH = A_ZEROH*VELOH+A_ONEH/VELOH-2*(A_ZEROH*A_ONEH)^0.5 ;
SHOPPRIMEH = A_ZEROH-A_ONEH*VELOH^(-2) ;
VELOH = CH/REALMONEYH ;
DEEH = 1-SHOPPRIMEH*VELOH^2 ;
1 = RNOMH*DEEH ;
100*(RNOMH^4-1) = (1-XDUMH)*100*(RNOM_EXOGH^4-1)+XDUMH*(XR3H*100*(RNOMH(-1)^4-1)+(1-XR3H)*(100*((1/BET*PIE4H^0.25)^4-1))+XR1H*(100*(PIE4H-1)-100*(PIE4TARH-1))+XR4H*100*(DEPEX^4-1)+XR2H*GDPGAPH) ;
GDPGAPH = 100*(GDPH_NAT-GDPH_EXOG)/GDPH_EXOG ;
PIE4H = PIEH*PIEH(-1)*PIEH1(-1)*PIEH2(-1) ;
AH = CH*(1+SHOPH)+EYEH+GAH ;
GAH = .05*(GA_RATH*(GLAMBDAH*GDPH_NAT+(1-GLAMBDAH)*GDPH_EXOG))+.95*GAH(-1)+E_GAH;
GNH = GN_RATH*(GLAMBDAH*GDPH_NAT+(1-GLAMBDAH)*GDPH_EXOG)/REALPNH_EXOG ;
PIEBARMH = PIEH*REALPBARMH/REALPBARMH(-1) ;
BIGGAMMH = LILCAPPA7H/2*(PIEBARMH/PIE4TARH^0.25-1)^2+LILCAPPA8H/2*(PIEBARMH/PIEBARMH(-1)-1)^2 ;
REALPMH = REALPBARMH+ETAH*REALPNH ;
KNH_RAT = ALPHANH/(1-ALPHANH-GAMA_NH)*(REALWH/REALRH)^XIXI_NH ;
KTH_RAT = ALPHATH/(1-ALPHATH-GAMA_TH)*(REALWH/REALRH)^XIXI_TH ;
KNH_RAT = KNH/LNH ;
KTH_RAT = KTH/LTH ;
KH = KTH+KNH+K_OH ;
LH = (LNH+LTH+L_OH)*(1-COSTLH/2*(LNH/(LTH+L_OH)/(LNH(-1)/(LTH(-1)+L_OH(-1)))-1)^2) ;
T_OH = Z_OH*((1-ALPHA_OH-GAMMA_LANDH)^(1/XIXI_OH)*L_OH^(1-1/XIXI_OH)+ALPHA_OH^(1/XIXI_OH)*K_OH^(1-1/XIXI_OH)+GAMMA_LANDH^(1/XIXI_OH)*LANDH^(1-1/XIXI_OH))^(XIXI_OH/(XIXI_OH-1)) ;
Q_ONH = NY_NH*(REALP_QOH/REALP_ONH)^(-EPS_ONH)*O_NH ;
Q_OTH = NY_TH*(REALP_QOH/REALP_OTH)^(-EPS_OTH)*O_TH ;
M_ONH = (1-NY_NH)*(REALP_MOH/REALP_ONH)^(-EPS_ONH)*O_NH*(BIGGAM_MONH-COST_MONH*(M_ONH/O_NH/(M_ONH(-1)/O_NH(-1))-1)*M_ONH/O_NH/(M_ONH(-1)/O_NH(-1)))^EPS_ONH/BIGGAM_MONH ;
M_OTH = (1-NY_TH)*(REALP_MOH/REALP_OTH)^(-EPS_OTH)*O_TH*(BIGGAM_MOTH-COST_MOTH*(M_OTH/O_TH/(M_OTH(-1)/O_TH(-1))-1)*M_OTH/O_TH/(M_OTH(-1)/O_TH(-1)))^EPS_OTH/BIGGAM_MOTH ;
BIGGAM_MONH = 1-COST_MONH/2*(M_ONH/O_NH/(M_ONH(-1)/O_NH(-1))-1)^2 ;
BIGGAM_MOTH = 1-COST_MOTH/2*(M_OTH/O_TH/(M_OTH(-1)/O_TH(-1))-1)^2 ;
K_OH_RAT = ALPHA_OH/(1-ALPHA_OH-GAMMA_LANDH)*(REALWH/REALRH)^XIXI_OH ;
K_OH_RAT = K_OH/L_OH ;
REALP_QOH = 1/Z_OH*((1-ALPHA_OH-GAMMA_LANDH)*REALWH^(1-XIXI_OH)+ALPHA_OH*REALRH^(1-XIXI_OH)+GAMMA_LANDH*REALPLANDH^(1-XIXI_OH))^(1/(1-XIXI_OH)) ;
LANDH = GAMMA_LANDH*(REALPLANDH/(REALP_QOH*Z_OH))^(-XIXI_OH)*T_OH/Z_OH ;
REALP_ONH = (NY_NH*REALP_QOH^(1-EPS_ONH)+(1-NY_NH)*REALP_MOH^(1-EPS_ONH)*(BIGGAM_MONH-COST_MONH*(M_ONH/O_NH/(M_ONH(-1)/O_NH(-1))-1)*M_ONH/O_NH/(M_ONH(-1)/O_NH(-1)))^(EPS_ONH-1))^(1/(1-EPS_ONH)) ;
REALP_OTH = (NY_TH*REALP_QOH^(1-EPS_OTH)+(1-NY_TH)*REALP_MOH^(1-EPS_OTH)*(BIGGAM_MOTH-COST_MOTH*(M_OTH/O_TH/(M_OTH(-1)/O_TH(-1))-1)*M_OTH/O_TH/(M_OTH(-1)/O_TH(-1)))^(EPS_OTH-1))^(1/(1-EPS_OTH)) ;
SSH*TH = SSH*QH+SSF*MF ;
SSH*T_OH = SSH*Q_ONH+SSH*Q_OTH+SSF*M_ONF+SSF*M_OTF ;
REALP_MOH = REALP_QOF*REALEXH ;
ZZ_GDP_PPP_RATH = GDPH/REALEX/GDPF ;
XI = CHI0*(exp(CHI1*REALEX*REALBH)+CHI2*(REALEX*(REALBH-REALBH(-1)/PIEF)/GDPH)^2+CHI3*(REALEX*(REALBH-REALBH(-1)/PIEF)/GDPH-REALEX(-1)*(REALBH(-1)-REALBH1(-1)/PIEF(-1))/GDPH(-1))^2-1)/(exp(CHI1*REALEX*REALBH)+CHI2*(REALEX*(REALBH-REALBH(-1)/PIEF)/GDPH)^2+CHI3*(REALEX*(REALBH-REALBH(-1)/PIEF)/GDPH-REALEX(-1)*(REALBH(-1)-REALBH1(-1)/PIEF(-1))/GDPH(-1))^2+1)+ZBH ;
1 = RNOMF*(1-XI)*DEEH*DEPEX(+1) ;
DEPEX = PIEH/PIEF*REALEX/REALEX(-1) ;
REALFINH = RNOMF(-1)*(1-XI(-1))*REALEX*REALBH(-1)/PIEF ;
SSH*DEEH*PIEH(+1)*REALFINH(+1) = SSH*REALFINH+SSH*RNOMF(-1)*XI(-1)*REALEX*REALBH(-1)/PIEF+REALTBALH ;
REALEXH = REALEX ;
REALEXF = 1/REALEXH ;
ZZ_REALEX = 100*log(REALEX) ;
-(1-BIGGAMMH)*REALPBARMH/REALPMH*(REALPBARMH/REALEX*(1-THETAF)+ETAH*REALPNH/REALEX+THETAF*REALMCTF) = -(REALPBARMH/REALEX-REALMCTF)*(LILCAPPA7H*PIEBARMH/PIE4TARH^0.25*(PIEBARMH/PIE4TARH^0.25-1)+LILCAPPA8H*PIEBARMH/PIEBARMH(-1)*(PIEBARMH/PIEBARMH(-1)-1))+DEEF*PIEF(+1)*(REALPBARMH(+1)/REALEX(+1)-REALMCTF(+1))*MH(+1)/MH*(LILCAPPA7H*PIEBARMH(+1)/PIE4TARH^0.25*(PIEBARMH(+1)/PIE4TARH^0.25-1)+LILCAPPA8H*PIEBARMH(+1)/PIEBARMH*(PIEBARMH(+1)/PIEBARMH-1)) ;
GDPH = AH+REALPNH*GNH+EXPORTSH-IMPORTSH+(RNOMF(-1)-1)*REALEX*REALBH(-1)/PIEF ;
GDPH_NAT = AH+REALPNH_EXOG*GNH+EXPORTSH_NAT-IMPORTSH_NAT ;
CH_NAT = CH*(1+SHOPH) ;
GH_NAT = GAH+REALPNH_EXOG*GNH ;
XH_NAT = SSF/SSH*REALEX_EXOG*REALPBARMF_EXOG*MF ;
MH_NAT = REALPBARMH_EXOG*MH ;
CURBALH_RAT = REALEX*(REALBH-REALBH(-1)/PIEF)/GDPH ;
REALTBALH = SSF*(REALPBARMF*MF+REALP_MOF*M_ONF+REALP_MOF*M_OTF)*REALEX-SSH*(REALPBARMH*MH+REALP_MOH*M_ONH+REALP_MOH*M_OTH) ;
EXPORTSH = SSF/SSH*(REALPBARMF*MF+REALP_MOF*M_ONF+REALP_MOF*M_OTF)*REALEX ;
IMPORTSH = REALPBARMH*MH+REALP_MOH*M_ONH+REALP_MOH*M_OTH ;
EXPORTSH_NAT = SSF/SSH*(REALPBARMF_EXOG*MF+REALP_MOF_EXOG*M_ONF+REALP_MOF_EXOG*M_OTF)*REALEX_EXOG ;
IMPORTSH_NAT = REALPBARMH_EXOG*MH+REALP_MOH_EXOG*M_ONH+REALP_MOH_EXOG*M_OTH ;
ZZ_UTILITYF = (ZUF*(CF-HF)^(1-SIGMAF)-1)/(1-SIGMAF)-CAPAF*LF^(1+ZEDF)/(1+ZEDF) ;
ZZ_GF = 100*log(GF_NAT) ;
ZZ_CURBALF_RAT = CURBALF_RAT*100 ;
ZZ_M_SHARF = REALPBARMF*MF/AF ;
ZZ_M_O_SHARF = (REALP_MOF*M_ONF+REALP_MOF*M_OTF)/AF ;
ZZ_M_ON_SHARF = REALP_MOF*M_ONF/AF ;
ZZ_M_OT_SHARF = REALP_MOF*M_OTF/AF ;
ZZ_N_SHARF = NF*REALPNF/AF ;
ZZ_EYE_SHARF = EYEF/GDPF ;
ZZ_C_SHARF = CF/GDPF ;
ZZ_GDPF = 100*log(GDPF_NAT) ;
ZZ_CF = 100*log(CF_NAT) ;
ZZ_EYEF = 100*log(EYEF) ;
ZZ_EXPORTSF = 100*log(EXPORTSF_NAT) ;
ZZ_IMPORTSF = 100*log(IMPORTSF_NAT) ;
ZZ_XBALF_TOT_RAT = 100*(EXPORTSF_NAT-IMPORTSF_NAT)/GDPF_NAT ;
ZZ_PIE4F = 100*(PIE4F-1) ;
ZZ_DPIE4F = ZZ_PIE4F-ZZ_PIE4F(-1) ;
ZZ_RNOMF = 100*(RNOMF^4-1) ;
ZZ_DRNOMF = ZZ_RNOMF-ZZ_RNOMF(-1) ;
100*(PIE4TARF-1) = 1*100*(PIE4TARF_SS-1)+(1-1)*100*(PIE4TARF(-1)-1)+E_PIE4TARF ;
log(ZUF) = 0.3*log(ZUF_SS)+0.7*log(ZUF(-1))+E_ZUF ;
ZBF = 0.3*ZBF_SS+0.7*ZBF(-1)+E_ZBF ;
log(LANDF) = 0.05*log(LANDF_SS)+0.95*log(LANDF(-1))+E_LANDF ;
log(ZTF) = 0.05*log(ZTF_SS)+0.95*log(ZTF(-1))+E_LANDF ;
log(ZNF) = 0.05*log(ZNF_SS)+0.95*log(ZNF(-1))+E_LANDF ;
log(Z_OF) = 0.05*log(Z_OF_SS)+0.95*log(Z_OF(-1))+E_LANDF ;
ZEYEF = 0.05*ZEYEF_SS+0.95*ZEYEF(-1)+E_ZEYEF ;
CAPAF = 0.05*CAPAF_SS+0.95*CAPAF(-1)+E_CAPAF ;
log(GAMMAF) = 0.05*log(GAMMAF_SS)+0.95*log(GAMMAF(-1))+E_GAMMAF ;
BIGGAM_O_NF = 1-COST_O_NF/2*(O_NF/NF/(O_NF(-1)/NF(-1))-1)^2 ;
BIGGAM_O_TF = 1-COST_O_TF/2*(O_TF/TF/(O_TF(-1)/TF(-1))-1)^2 ;
O_NF = GAMA_NF*NF/ZNF*(REALP_ONF/(REALMCNF*ZNF))^(-XIXI_NF)*(BIGGAM_O_NF-COST_O_NF*(O_NF/NF/(O_NF(-1)/NF(-1))-1)*O_NF/NF/(O_NF(-1)/NF(-1)))^XIXI_NF/BIGGAM_O_NF ;
O_TF = GAMA_TF*TF/ZTF*(REALP_OTF/(REALMCTF*ZTF))^(-XIXI_TF)*(BIGGAM_O_TF-COST_O_TF*(O_TF/TF/(O_TF(-1)/TF(-1))-1)*O_TF/TF/(O_TF(-1)/TF(-1)))^XIXI_NF/BIGGAM_O_TF ;
NF = ZNF*((1-ALPHANF-GAMA_NF)^(1/XIXI_NF)*LNF^(1-1/XIXI_NF)+ALPHANF^(1/XIXI_NF)*KNF^(1-1/XIXI_NF)+GAMA_NF^(1/XIXI_NF)*(BIGGAM_O_NF*O_NF)^(1-1/XIXI_NF))^(XIXI_NF/(XIXI_NF-1)) ;
TF = ZTF*((1-ALPHATF-GAMA_TF)^(1/XIXI_TF)*LTF^(1-1/XIXI_TF)+ALPHATF^(1/XIXI_TF)*KTF^(1-1/XIXI_TF)+GAMA_TF^(1/XIXI_TF)*(BIGGAM_O_TF*O_TF)^(1-1/XIXI_TF))^(XIXI_TF/(XIXI_TF-1)) ;
REALMCNF = 1/ZNF*((1-ALPHANF-GAMA_NF)*REALWF^(1-XIXI_NF)+ALPHANF*REALRF^(1-XIXI_NF)+GAMA_NF*REALP_ONF^(1-XIXI_NF)*(BIGGAM_O_NF-COST_O_NF*(O_NF/NF/(O_NF(-1)/NF(-1))-1)*O_NF/NF/(O_NF(-1)/NF(-1)))^(XIXI_NF-1))^(1/(1-XIXI_NF)) ;
REALMCTF = 1/ZTF*((1-ALPHATF-GAMA_TF)*REALWF^(1-XIXI_TF)+ALPHATF*REALRF^(1-XIXI_TF)+GAMA_TF*REALP_OTF^(1-XIXI_TF)*(BIGGAM_O_TF-COST_O_TF*(O_TF/TF/(O_TF(-1)/TF(-1))-1)*O_TF/TF/(O_TF(-1)/TF(-1)))^(XIXI_TF-1))^(1/(1-XIXI_TF)) ;
MARGUTF = (CF-B2F*HF)^(-SIGMAF)*ZUF ;
HF = (1-B0F)*HF(-1)+B0F*B1F*CF(-1) ;
VPRIMEF = CAPAF*LF^ZEDF ;
AF^(1-1/EPSF) = GAMMAF^(1/EPSF)*(NYF^(1/EPSQMF)*QF^(1-1/EPSQMF)+(1-NYF)^(1/EPSQMF)*(MF*BIGGAMIMPF)^(1-1/EPSQMF))^(EPSQMF/(EPSQMF-1)*(1-1/EPSF))+(1-GAMMAF)^(1/EPSF)*NNF^(1-1/EPSF) ;
QF = GAMMAF*NYF*REALPQF^(-EPSQMF)*AF*REALPXF^(EPSQMF-EPSF) ;
MF = GAMMAF*(1-NYF)*REALPMF^(-EPSQMF)*AF*REALPXF^(EPSQMF-EPSF)*1/BIGGAMIMPF*(BIGGAMIMPF-COSTF*(MF/AF/(MF(-1)/AF(-1))-1)*MF/AF/(MF(-1)/AF(-1)))^EPSQMF ;
REALPXF = (NYF*REALPQF^(1-EPSQMF)+(1-NYF)*REALPMF^(1-EPSQMF)*(BIGGAMIMPF-COSTF*(MF/AF/(MF(-1)/AF(-1))-1)*MF/AF/(MF(-1)/AF(-1)))^(EPSQMF-1))^(1/(1-EPSQMF)) ;
BIGGAMIMPF = 1-COSTF/2*(MF/AF/(MF(-1)/AF(-1))-1)^2 ;
NNF = (1-GAMMAF)*REALPNF^(-EPSF)*AF ;
NF = NNF+ETAF*MF+ETAF*QF+GNF ;
PIENF = REALPNF/REALPNF(-1)*PIEF ;
BIGGAMNF = LILCAPPA3F/2*(PIENF/PIE4TARF^0.25-1)^2+LILCAPPA4F/2*(PIENF/PIENF(-1)-1)^2 ;
-(1-BIGGAMNF)*(REALPNF*(1-THETAF)+THETAF*REALMCNF) = -(REALPNF-REALMCNF)*(LILCAPPA3F*PIENF/PIE4TARF^0.25*(PIENF/PIE4TARF^0.25-1)+LILCAPPA4F*PIENF/PIENF(-1)*(PIENF/PIENF(-1)-1))+DEEF*PIEF(+1)*(REALPNF(+1)-REALMCNF(+1))*NF(+1)/NF*(LILCAPPA3F*PIENF(+1)/PIE4TARF^0.25*(PIENF(+1)/PIE4TARF^0.25-1)+LILCAPPA4F*PIENF(+1)/PIENF*(PIENF(+1)/PIENF-1)) ;
PIEBARQF = PIEF*REALPBARQF/REALPBARQF(-1) ;
BIGGAMQF = LILCAPPA5F/2*(PIEBARQF/PIE4TARF^0.25-1)^2+LILCAPPA6F/2*(PIEBARQF/PIEBARQF(-1)-1)^2 ;
-(1-BIGGAMQF)*REALPBARQF/REALPQF*(REALPBARQF*(1-THETAF)+ETAF*REALPNF+THETAF*REALMCTF) = -(REALPBARQF-REALMCTF)*(LILCAPPA5F*PIEBARQF/PIE4TARF^0.25*(PIEBARQF/PIE4TARF^0.25-1)+LILCAPPA6F*PIEBARQF/PIEBARQF(-1)*(PIEBARQF/PIEBARQF(-1)-1))+DEEF*PIEF(+1)*(REALPBARQF(+1)-REALMCTF(+1))*QF(+1)/QF*(LILCAPPA5F*PIEBARQF(+1)/PIE4TARF^0.25*(PIEBARQF(+1)/PIE4TARF^0.25-1)+LILCAPPA6F*PIEBARQF(+1)/PIEBARQF*(PIEBARQF(+1)/PIEBARQF-1)) ;
REALPQF = REALPBARQF+ETAF*REALPNF ;
KF = KF(-1)*(1-DELTAF)+PSIF(-1)*KF(-1) ;
PSIF = EYEF/KF-OMEGAF/2*(EYEF/KF-DELTAF*(1+ZEYEF))^2-OMEGA0F/2*(EYEF/KF-EYEF(-1)/KF(-1))^2 ;
PSIPRIMEF = 1-OMEGAF*(EYEF/KF-DELTAF*(1+ZEYEF))-OMEGA0F*(EYEF/KF-EYEF(-1)/KF(-1)) ;
1/PSIPRIMEF = DEEF*PIEF(+1)*(REALRF(+1)+1/PSIPRIMEF(+1)*(1-DELTAF+PSIF(+1)*(1-PSIPRIMEF(+1)*EYEF(+1)/(PSIF(+1)*KF(+1))))) ;
BIGGAMF = LILCAPPA1F/2*(PIEWF/PIE4TARF^0.25-1)^2+LILCAPPA2F/2*(PIEWF/PIEWF(-1)-1)^2 ;
PIEF*REALWF/REALWF(-1) = PIEWF ;
REALWF = PHIF*VPRIMEF/MARGUTF*((PHIF-1)*(1-BIGGAMF)+PIEWF*LILCAPPA1F/PIE4TARF^0.25*(PIEWF/PIE4TARF^0.25-1)+PIEWF/PIEWF(-1)*LILCAPPA2F*(PIEWF/PIEWF(-1)-1)-DEEF*PIEWF(+1)*LF(+1)/LF*LILCAPPA1F*PIEWF(+1)/PIE4TARF^0.25*(PIEWF(+1)/PIE4TARF^0.25-1)-DEEF*PIEWF(+1)*LF(+1)/LF*LILCAPPA2F*PIEWF(+1)/(REALWF/REALWF(-1))*(PIEWF(+1)/PIEWF-1))^(-1) ;
DEEF = BET*MARGUTF(+1)/PIEF(+1)/MARGUTF*(1+SHOPF+SHOPPRIMEF*VELOF)/(1+SHOPF(+1)+SHOPPRIMEF(+1)*VELOF(+1)) ;
SHOPF = A_ZEROF*VELOF+A_ONEF/VELOF-2*(A_ZEROF*A_ONEF)^0.5 ;
SHOPPRIMEF = A_ZEROF-A_ONEF*VELOF^(-2) ;
VELOF = CF/REALMONEYF ;
DEEF = 1-SHOPPRIMEF*VELOF^2 ;
1 = RNOMF*DEEF ;
100*(RNOMF^4-1) = (1-XDUMF)*100*(RNOM_EXOGF^4-1)+XDUMF*(XR3F*100*(RNOMF(-1)^4-1)+(1-XR3F)*(100*((1/BET*PIE4F^0.25)^4-1)+XR1F*(100*(PIE4F-1)-100*(PIE4TARF-1))+XR4F*100*(DEPEX^4-1)+XR2F*GDPGAPF)) ;
GDPGAPF = 100*(GDPF_NAT-GDPF_EXOG)/GDPF_EXOG ;
PIE4F = PIEF*PIEF(-1)*PIEF1(-1)*PIEF2(-1) ;
AF = CF*(1+SHOPF)+EYEF+GAF ;
GAF = .05*(GA_RATF*(GLAMBDAF*GDPF_NAT+(1-GLAMBDAF)*GDPF_EXOG))+.95*GAF(-1)+E_GAF;
GNF = GN_RATF*(GLAMBDAF*GDPF_NAT+(1-GLAMBDAF)*GDPF_EXOG)/REALPNF_EXOG ;
PIEBARMF = PIEF*REALPBARMF/REALPBARMF(-1) ;
BIGGAMMF = LILCAPPA7F/2*(PIEBARMF/PIE4TARF^0.25-1)^2+LILCAPPA8F/2*(PIEBARMF/PIEBARMF(-1)-1)^2 ;
REALPMF = REALPBARMF+ETAF*REALPNF ;
KNF_RAT = ALPHANF/(1-ALPHANF-GAMA_NF)*(REALWF/REALRF)^XIXI_NF ;
KTF_RAT = ALPHATF/(1-ALPHATF-GAMA_TF)*(REALWF/REALRF)^XIXI_TF ;
KNF_RAT = KNF/LNF ;
KTF_RAT = KTF/LTF ;
KF = KTF+KNF+K_OF ;
LF = (LNF+LTF+L_OF)*(1-COSTLF/2*(LNF/(LTF+L_OF)/(LNF(-1)/(LTF(-1)+L_OF(-1)))-1)^2) ;
T_OF = Z_OF*((1-ALPHA_OF-GAMMA_LANDF)^(1/XIXI_OF)*L_OF^(1-1/XIXI_OF)+ALPHA_OF^(1/XIXI_OF)*K_OF^(1-1/XIXI_OF)+GAMMA_LANDF^(1/XIXI_OF)*LANDF^(1-1/XIXI_OF))^(XIXI_OF/(XIXI_OF-1)) ;
Q_ONF = NY_NF*(REALP_QOF/REALP_ONF)^(-EPS_ONF)*O_NF ;
Q_OTF = NY_TF*(REALP_QOF/REALP_OTF)^(-EPS_OTF)*O_TF ;
M_ONF = (1-NY_NF)*(REALP_MOF/REALP_ONF)^(-EPS_ONF)*O_NF*(BIGGAM_MONF-COST_MONF*(M_ONF/O_NF/(M_ONF(-1)/O_NF(-1))-1)*M_ONF/O_NF/(M_ONF(-1)/O_NF(-1)))^EPS_ONF/BIGGAM_MONF ;
M_OTF = (1-NY_TF)*(REALP_MOF/REALP_OTF)^(-EPS_OTF)*O_TF*(BIGGAM_MOTF-COST_MOTF*(M_OTF/O_TF/(M_OTF(-1)/O_TF(-1))-1)*M_OTF/O_TF/(M_OTF(-1)/O_TF(-1)))^EPS_OTF/BIGGAM_MOTF ;
BIGGAM_MONF = 1-COST_MONF/2*(M_ONF/O_NF/(M_ONF(-1)/O_NF(-1))-1)^2 ;
BIGGAM_MOTF = 1-COST_MOTF/2*(M_OTF/O_TF/(M_OTF(-1)/O_TF(-1))-1)^2 ;
K_OF_RAT = ALPHA_OF/(1-ALPHA_OF-GAMMA_LANDF)*(REALWF/REALRF)^XIXI_OF ;
K_OF_RAT = K_OF/L_OF ;
REALP_QOF = 1/Z_OF*((1-ALPHA_OF-GAMMA_LANDF)*REALWF^(1-XIXI_OF)+ALPHA_OF*REALRF^(1-XIXI_OF)+GAMMA_LANDF*REALPLANDF^(1-XIXI_OF))^(1/(1-XIXI_OF)) ;
LANDF = GAMMA_LANDF*(REALPLANDF/(REALP_QOF*Z_OF))^(-XIXI_OF)*T_OF/Z_OF ;
REALP_ONF = (NY_NF*REALP_QOF^(1-EPS_ONF)+(1-NY_NF)*REALP_MOF^(1-EPS_ONF)*(BIGGAM_MONF-COST_MONF*(M_ONF/O_NF/(M_ONF(-1)/O_NF(-1))-1)*M_ONF/O_NF/(M_ONF(-1)/O_NF(-1)))^(EPS_ONF-1))^(1/(1-EPS_ONF)) ;
REALP_OTF = (NY_TF*REALP_QOF^(1-EPS_OTF)+(1-NY_TF)*REALP_MOF^(1-EPS_OTF)*(BIGGAM_MOTF-COST_MOTF*(M_OTF/O_TF/(M_OTF(-1)/O_TF(-1))-1)*M_OTF/O_TF/(M_OTF(-1)/O_TF(-1)))^(EPS_OTF-1))^(1/(1-EPS_OTF)) ;
SSF*TF = SSF*QF+SSH*MH ;
SSF*T_OF = SSF*Q_ONF+SSF*Q_OTF+SSH*M_ONH+SSH*M_OTH ;
REALP_MOF = REALP_QOH*REALEXF ;
SSH*REALBH+SSF*REALBF = 0 ;
REALTBALF = SSF*(REALPBARMF*MF+REALP_MOF*M_ONF+REALP_MOF*M_OTF)-SSH*(REALPBARMH*MH+REALP_MOH*M_ONH+REALP_MOH*M_OTH)*1/REALEX ;
EXPORTSF = SSH/SSF*(REALPBARMH*MH+REALP_MOH*M_ONH+REALP_MOH*M_OTH)*1/REALEX ;
IMPORTSF = REALPBARMF*MF+REALP_MOF*M_ONF+REALP_MOF*M_OTF ;
EXPORTSF_NAT = SSH/SSF*(REALPBARMH_EXOG*MH+REALP_MOH_EXOG*M_ONH+REALP_MOH_EXOG*M_OTH)*1/REALEX_EXOG ;
IMPORTSF_NAT = REALPBARMF_EXOG*MF+REALP_MOF_EXOG*M_ONF+REALP_MOF_EXOG*M_OTF ;
-(1-BIGGAMMF)*REALPBARMF/REALPMF*(REALPBARMF*REALEX*(1-THETAH)+ETAF*REALPNF*REALEX+THETAH*REALMCTH) = -(REALPBARMF*REALEX-REALMCTH)*(LILCAPPA7F*PIEBARMF/PIE4TARF^0.25*(PIEBARMF/PIE4TARF^0.25-1)+LILCAPPA8F*PIEBARMF/PIEBARMF(-1)*(PIEBARMF/PIEBARMF(-1)-1))+DEEH*PIEH(+1)*(REALPBARMF(+1)*REALEX(+1)-REALMCTH(+1))*MF(+1)/MF*(LILCAPPA7F*PIEBARMF(+1)/PIE4TARF^0.25*(PIEBARMF(+1)/PIE4TARF^0.25-1)+LILCAPPA8F*PIEBARMF(+1)/PIEBARMF*(PIEBARMF(+1)/PIEBARMF-1)) ;
GDPF = AF+REALPNF*GNF+EXPORTSF-IMPORTSF+(RNOMF(-1)-1)*REALBF(-1)/PIEF ;
GDPF_NAT = AF+REALPNF_EXOG*GNF+EXPORTSF_NAT-IMPORTSF_NAT ;
CF_NAT = CF*(1+SHOPF) ;
GF_NAT = GAF+REALPNF_EXOG*GNF ;
XF_NAT = SSH/SSF*1/REALEX_EXOG*REALPBARMH_EXOG*MH ;
MF_NAT = REALPBARMF_EXOG*MF ;
CURBALF_RAT = -(REALTBALH/REALEX/SSF/GDPF)+(RNOMF(-1)-1)*REALBF(-1)/PIEF/GDPF ;
PIEF1 = PIEF(-1);
PIEF2 = PIEF1(-1);
PIEH1 = PIEH(-1);
PIEH2 = PIEH1(-1);
REALBH1 = REALBH(-1);
end;
initval;
AF=2.17350447531715;
AH=2.61461230039988;
BIGGAMF=0;
BIGGAMH=0;
BIGGAMIMPF=1;
BIGGAMIMPH=1;
BIGGAMMF=0;
BIGGAMMH=0;
BIGGAMNF=0;
BIGGAMNH=0;
BIGGAMQF=0;
BIGGAMQH=0;
BIGGAM_MONF=1;
BIGGAM_MONH=1;
BIGGAM_MOTF=1;
BIGGAM_MOTH=1;
BIGGAM_O_NF=1;
BIGGAM_O_NH=1;
BIGGAM_O_TF=1;
BIGGAM_O_TH=1;
CAPAF=11;
CAPAH=11;
CF=1.77599320017707;
CF_NAT=1.77797456682707;
CH=2.10139281352027;
CH_NAT=2.10373720855446;
CURBALF_RAT=2.20209042676066e-018;
CURBALH_RAT=0;
DEEF=0.963834712172592;
DEEH=0.963834712172592;
DEPEX=1;
EXPORTSF=0.0374229290542059;
EXPORTSF_NAT=0.0374229290542059;
EXPORTSH=0.976573287861717;
EXPORTSH_NAT=0.976573287861717;
EYEF=0.27477965986135;
EYEH=0.365618852934316;
GAF=0.12075024862873;
GAH=0.145256238911104;
GAMMAF=0.5;
GAMMAH=0.25;
GDPF=2.41500497257461;
GDPF_NAT=2.41500497257461;
GDPGAPF=0;
GDPGAPH=0;
GDPH=2.90512477822209;
GDPH_NAT=2.90512477822209;
GF_NAT=0.362250745886191;
GH_NAT=0.435768716733313;
GNF=0.287269571519256;
GNH=0.321902361090147;
HF=1.68719354016822;
HH=1.99632317284426;
IMPORTSF=0.0374229290542059;
IMPORTSF_NAT=0.0374229290542059;
IMPORTSH=0.976573287861718;
IMPORTSH_NAT=0.976573287861718;
KF=10.991186394454;
KH=14.6247541173726;
KNF=6.33686501417153;
KNF_RAT=22.6981730731029;
KNH=11.034700665508;
KNH_RAT=22.8755992006951;
KTF=2.97137434524903;
KTF_RAT=22.6981730731029;
KTH=2.23720856941572;
KTH_RAT=114.377996003476;
K_OF=1.68294703503345;
K_OF_RAT=7.27127622255245;
K_OH=1.35284488244891;
K_OH_RAT=8.16985685739111;
LANDF=0.1;
LANDH=0.1;
LF=0.64153899810027;
LH=0.667528221502678;
LNF=0.279179517830034;
LNH=0.482378650224502;
LTF=0.130908083909629;
LTH=0.019559781143112;
L_OF=0.231451396360608;
L_OH=0.165589790135064;
MARGUTF=2.24145263303312;
MARGUTH=2.11921125101343;
MF=0.0196445696804563;
MF_NAT=0.0171196449669319;
MH=0.438784845846124;
MH_NAT=0.522472906750236;
M_ONF=0.0143006671963624;
M_ONH=0.134410532365428;
M_OTF=0.00670562423725087;
M_OTH=0.143002828997546;
NF=1.91582345366461;
NH=2.609674642079;
NNF=1.31534385473198;
NNH=2.19524942542191;
O_NF=0.387338325509274;
O_NH=0.147043832240678;
O_TF=0.18162406186278;
O_TH=0.148205762233076;
PIE4F=1.125;
PIE4H=1.125;
PIE4TARF=1.125;
PIE4TARH=1.125;
PIEBARMF=1.02988357195356;
PIEBARMH=1.02988357195356;
PIEBARQF=1.02988357195356;
PIEBARQH=1.02988357195356;
PIEF=1.02988357195356;
PIEF1=1.02988357195356;
PIEF2=1.02988357195356;
PIEH=1.02988357195356;
PIEH1=1.02988357195356;
PIEH2=1.02988357195356;
PIENF=1.02988357195356;
PIENH=1.02988357195356;
PIEWF=1.02988357195356;
PIEWH=1.02988357195356;
PSIF=0.025;
PSIH=0.025;
PSIPRIMEF=1;
PSIPRIMEH=1;
QF=0.875241222929181;
QH=0.0238294319885835;
Q_ONF=0.373740369418894;
Q_ONH=0.0132636199615755;
Q_OTF=0.175247940896905;
Q_OTH=0.00547180886242481;
REALBF=0;
REALBH=0;
REALBH1=0;
REALEX=1.3734519289908;
REALEXF=0.728092464608345;
REALEXH=1.3734519289908;
REALFINH=0;
REALMCNF=0.700562935771035;
REALMCNH=0.752071934789911;
REALMCTF=0.700562935771035;
REALMCTH=0.930081384894704;
REALMONEYF=0.558667031035572;
REALMONEYH=0.661026677383566;
REALPBARMF=0.87146958398196;
REALPBARMH=1.19072687148694;
REALPBARQF=0.899522809530009;
REALPBARQH=1.15219711474356;
REALPLANDF=0.554831427212494;
REALPLANDH=0.414697221827051;
REALPMF=1.16570601700579;
REALPMH=1.37122413583652;
REALPNF=0.840675522925242;
REALPNH=0.902486321747893;
REALPQF=1.19375924255384;
REALPQH=1.33269437909314;
REALPXF=1.19317131724075;
REALPXH=1.36926881180313;
REALP_MOF=0.966533486000563;
REALP_MOH=1.63690883121281;
REALP_ONF=1.18566549908199;
REALP_ONH=1.61601524261254;
REALP_OTF=1.18566549908199;
REALP_OTH=1.62845456685201;
REALP_QOF=1.1918209852569;
REALP_QOH=1.32748728078168;
REALRF=0.0324170717777328;
REALRH=0.0324170717777329;
REALTBALF=-6.93889390390723e-018;
REALTBALH=-6.93889390390723e-018;
REALWF=2.42667732699502;
REALWH=2.83454771236558;
RNOMF=1.03752229232945;
RNOMH=1.03752229232945;
SHOPF=0.00111563864647424;
SHOPH=0.00111563864647424;
SHOPPRIMEF=0.00357861859467432;
SHOPPRIMEH=0.00357861859467432;
TF=0.89833516218424;
TH=0.397076255917254;
T_OF=0.563589013545429;
T_OH=0.417854966062653;
VELOF=3.17898336847443;
VELOH=3.17898336847443;
VPRIMEF=3.62618818940983;
VPRIMEH=4.00467026905301;
XF_NAT=0.0200215045456245;
XH_NAT=0.446747178665936;
XI=0;
ZBF=0;
ZBH=0;
ZEYEF=0;
ZEYEH=0;
ZNF=1;
ZNH=1;
ZTF=1;
ZTH=0.6;
ZUF=1;
ZUH=1;
ZZ_CF=57.5474832617676;
ZZ_CH=74.3715386197541;
ZZ_CURBALF_RAT=2.20209042676066e-016;
ZZ_CURBALH_RAT=0;
ZZ_C_SHARF=0.735399396831762;
ZZ_C_SHARH=0.723339950584259;
ZZ_DPIE4F=0;
ZZ_DPIE4H=0;
ZZ_DRNOMF=0;
ZZ_DRNOMH=0;
ZZ_EXPORTSF=-328.547168610049;
ZZ_EXPORTSH=-2.37054799079326;
ZZ_EYEF=-129.17857393452;
ZZ_EYEH=-100.616387362469;
ZZ_EYE_SHARF=0.113780163180538;
ZZ_EYE_SHARH=0.12585306341233;
ZZ_GDPF=88.1701346139521;
ZZ_GDPH=106.647634229781;
ZZ_GDP_PPP_RATH=0.875857186130553;
ZZ_GF=-101.541863874636;
ZZ_GH=-83.0643642588075;
ZZ_IMPORTSF=-328.547168610049;
ZZ_IMPORTSH=-2.37054799079323;
ZZ_M_ON_SHARF=0.0063593490946998;
ZZ_M_ON_SHARH=0.084149297164759;
ZZ_M_OT_SHARF=0.00298191719568198;
ZZ_M_OT_SHARH=0.0895286056899133;
ZZ_M_O_SHARF=0.00934126629038178;
ZZ_M_O_SHARH=0.173677902854672;
ZZ_M_SHARF=0.00787651700806085;
ZZ_M_SHARH=0.19982806118916;
ZZ_N_SHARF=0.741008772713445;
ZZ_N_SHARH=0.90078198910348;
ZZ_PIE4F=12.5;
ZZ_PIE4H=12.5;
ZZ_REALEX=31.7327227026121;
ZZ_RNOMF=15.8749999999999;
ZZ_RNOMH=15.8749999999999;
ZZ_UTILITYF=-1.86610854895021;
ZZ_UTILITYH=-1.9297829736965;
ZZ_XBALF_TOT_RAT=0;
ZZ_XBALH_TOT_RAT=-7.6432037132987e-015;
Z_OF=1;
Z_OH=1;
E_ZBH=0;
E_ZUH=0;
E_ZUF=0;
E_ZEYEH=0;
E_ZEYEF=0;
E_GAMMAH=0;
E_GAMMAF=0;
E_LANDH=0;
E_LANDF=0;
E_GAH = 0;
E_GAF = 0;
E_CAPAH=0;
E_CAPAF=0;
end;
vcov = [
0.000324 0 0 0 0 0 0 0 0 0 0 0 0;
0 0.0004 0 0 0 0 0 0 0 0 0 0 0;
0 0 0.00000001 0 0 0 0 0 0 0 0 0 0;
0 0 0 0.000004 0 0 0 0 0 0 0 0 0;
0 0 0 0 0.000289 0 0 0 0 0 0 0 0;
0 0 0 0 0 0.000025 0 0 0 0 0 0 0;
0 0 0 0 0 0 0.0049 0 0 0 0 0 0;
0 0 0 0 0 0 0 0.000001 0 0 0 0 0;
0 0 0 0 0 0 0 0 0.000016 0 0 0 0;
0 0 0 0 0 0 0 0 0 0.00001225 0 0 0;
0 0 0 0 0 0 0 0 0 0 0.0000005625 0 0;
0 0 0 0 0 0 0 0 0 0 0 0.01 0;
0 0 0 0 0 0 0 0 0 0 0 0 0.0001
];
order=2;
\ No newline at end of file
dynare++/tests/judd.mod deleted 100644 → 0
View file @ 219d2bb3
var C K1 K2 L S1 S2 THETA V V1 V2;
varexo KSI;
parameters theta_ss lambda delta1 delta2 alpha1 alpha2 eta beta gamma depr1 depr2;
theta_ss=1;
lambda=0.8;
delta1=0.1;
delta2=0.05;
alpha1=0.3;
alpha2=0.15;
eta=3;
beta=0.95;
gamma=0.5;
depr1=0.1;
depr2=0.05;
model;
C = THETA*K1^alpha1*K2^alpha2*L^(1-alpha1-alpha2)-S1*K1-S2*K2;
K1 = (1-depr1+(1-0.5*delta1*S1)*S1)*K1(-1);
K2 = (1-depr2+(1-0.5*delta2*S2)*S2)*K2(-1);
THETA = THETA(-1)^lambda*theta_ss^(1-lambda)*exp(KSI);
/*
THETA = THETA(-1)*lambda+theta_ss*(1-lambda)+KSI;
*/
C^(-gamma)*THETA*K1^alpha1*K2^alpha2*L^(-alpha1-alpha2)*(1-alpha1-alpha2) = L^eta;
C^(-gamma) = beta*V1(+1)*(1-delta1*S1);
C^(-gamma) = beta*V2(+1)*(1-delta2*S2);
V1 = C^(-gamma)*(alpha1*THETA*K1^(alpha1-1)*K2^alpha2*L^(1-alpha1-alpha2)-S1)+beta*V1(+1)*(1-depr1+(1-0.5*delta1*S1)*S1);
V2 = C^(-gamma)*(alpha2*THETA*K1^alpha1*K2^(alpha2-1)*L^(1-alpha1-alpha2)-S2)+beta*V2(+1)*(1-depr2+(1-0.5*delta2*S2)*S2);
V = (C^(1-gamma)/(1-gamma)-L^(1+eta)/(1+eta)) + beta*V(+1);
end;
initval;
C= 1.33341818203972;
K1= 3.80023995548668;
K2= 3.80023995563911;
L= 0.85120255261552;
S1= 0;
S2= 0;
THETA= 1.00000000000000;
V1= 0.59202988402399;
V2= 0.59202988402399;
V= -17.6239;
end;
vcov = [ 0.001 ];
order = 6;
dynare++/tests/judd_norm.mod deleted 100644 → 0
View file @ 219d2bb3
var C K1 K2 L S1 S2 THETA V V1 V2;
varexo KSI;
parameters theta_ss lambda delta1 delta2 alpha1 alpha2 eta beta gamma depr1 depr2;
theta_ss=1;
lambda=0.5;
delta1=0.05;
delta2=0.2;
alpha1=0.3;
alpha2=0.3;
eta=3;
beta=0.95;
gamma=0.5;
depr1=0.1;
depr2=0.05;
model;
1 = (THETA*K1^alpha1*K2^alpha2*L^(1-alpha1-alpha2)-S1*K1-S2*K2)/C;
1 = (1-depr1+(1-0.5*delta1*S1)*S1)*K1(-1)/K1;
1 = (1-depr2+(1-0.5*delta2*S2)*S2)*K2(-1)/K2;
1 = THETA(-1)^lambda/THETA*theta_ss^(1-lambda)*exp(KSI);
/*
1 = (THETA(-1)*lambda+theta_ss*(1-lambda)+KSI)/THETA;
*/
C^(-gamma)*THETA*K1^alpha1*K2^alpha2*L^(-alpha1-alpha2)*(1-alpha1-alpha2)*L^(-eta)=1;
1 = beta*V1(+1)*(1-delta1*S1)*C^gamma;
1 = beta*V2(+1)*(1-delta2*S2)*C^gamma;
1 = (C^(-gamma)*(alpha1*THETA*K1^(alpha1-1)*K2^alpha2*L^(1-alpha1-alpha2)-S1)+beta*V1(+1)*(1-depr1+(1-0.5*delta1*S1)*S1))/V1;
1 = (C^(-gamma)*(alpha2*THETA*K1^alpha1*K2^(alpha2-1)*L^(1-alpha1-alpha2)-S2)+beta*V2(+1)*(1-depr2+(1-0.5*delta2*S2)*S2))/V2;
1 = (C^(1-gamma)/(1-gamma)-L^(1+eta)/(1+eta) + beta*V(+1))/V;
end;
initval;
C =1.0997055 ;
L =0.9425540 ;
S1 =0.1005051 ;
S2 =0.0500627 ;
K1 =2.9378521 ;
K2 =2.1952681 ;
THETA =1. ;
V =38.000392 ;
V1 =1.0139701 ;
V2 =1.0062981 ;
end;
vcov = [ 0.05 ];
order = 5;
dynare++/tests/kp1980_1.dyn deleted 100644 → 0
View file @ 219d2bb3
// Model from Kydland & Prescott JEDC 1980
// case 1: optimal policy, in fact, optimal control
var C G K TAU Z;
varexo EPS;
parameters eta beta alpha delta phi a rho;
eta = 2;
beta = 0.99;
alpha = 0.3;
delta = 0.10;
phi = 2.5;
a = 0.1;
rho = 0.7;
planner_objective C^(1-eta)/(1-eta) + a*G^(1-phi)/(1-phi);
planner_discount beta;
model;
K = (1-delta)*K(-1) + (exp(Z(-1))*K(-1)^alpha - C(-1) - G(-1));
G = TAU*alpha*K^alpha;
Z = rho*Z(-1) + EPS;
end;
initval;
TAU = 0.70;
K = ((delta+1/beta-1)/alpha)^(1/(alpha-1));
G = TAU*alpha*K^alpha;
C = K^alpha - delta*K - G;
Z = 0;
end;
order = 4;
vcov = [
0.01
];
dynare++/tests/kp1980_2.dyn deleted 100644 → 0
View file @ 219d2bb3
// Model from Kydland & Prescott JEDC 1980
// case 2: time inconsistent optimal policy with different policy and consumer objectives
var C G K TAU Z;
varexo EPS;
parameters eta beta alpha delta phi a rho;
eta = 2;
beta = 0.99;
alpha = 0.3;
delta = 0.10;
phi = 2.5;
a = 0.1;
rho = 0.7;
planner_objective C^(1-eta)/(1-eta) + a*G^(1-phi)/(1-phi);
planner_discount beta;
model;
K = (1-delta)*K(-1) + (exp(Z(-1))*K(-1)^alpha - C(-1) - G(-1));
G = TAU*alpha*K^alpha;
Z = rho*Z(-1) + EPS;
C^(-eta) = beta*C(+1)^(-eta)*(1-delta+exp(Z(+1))*alpha*K(+1)^(alpha-1)*(1-alpha*TAU(+1)));
end;
initval;
TAU = 0.70;
K = ((delta+1/beta-1)/(alpha*(1-alpha*TAU)))^(1/(alpha-1));
G = TAU*alpha*K^alpha;
C = K^alpha - delta*K - G;
Z = 0;
end;
order = 4;
vcov = [
0.01
];
dynare++/tests/kp1980_3.dyn deleted 100644 → 0
View file @ 219d2bb3
// Model from Kydland & Prescott JEDC 1980
// case 3: optimal policy with consistent objective (equivalent to kp1980_1.dyn)
var C G K TAU Z;
varexo EPS;
parameters eta beta alpha delta phi a rho;
eta = 2;
beta = 0.99;
alpha = 0.3;
delta = 0.10;
phi = 2.5;
a = 0.1;
rho = 0.7;
planner_objective C^(1-eta)/(1-eta) + a*G^(1-phi)/(1-phi);
planner_discount beta;
model;
K = (1-delta)*K(-1) + (exp(Z(-1))*K(-1)^alpha - C(-1) - G(-1));
G = TAU*alpha*K^alpha;
Z = rho*Z(-1) + EPS;
C^(-eta) = beta*C(+1)^(-eta)*(1-delta+exp(Z(+1))*alpha*K(+1)^(alpha-1)*(1-alpha*TAU(+1))) + beta*a*G(+1)^(-phi)*TAU(+1)*exp(Z(+1))*alpha^2*K(+1)^(alpha-1);
end;
initval;
TAU = 0.70;
K = ((delta+1/beta-1)/alpha)^(1/(alpha-1));
G = TAU*alpha*K^alpha;
C = K^alpha - delta*K - G;
Z = 0;
end;
order = 4;
vcov = [
0.01
];
dynare++/tests/lucas78.mod deleted 100644 → 0
View file @ 219d2bb3
var Y P;
varexo EXO_Y;
parameters beta gamma rho y_ss;
beta = 0.95;
gamma= 0.5;
rho = 0.9;
y_ss = 2;
model;
Y-y_ss = rho*(Y(-1)-y_ss) + EXO_Y;
Y^(-gamma)*P = beta*Y(+1)^(-gamma)*(P(+1) + Y(+1));
end;
initval;
Y = 2;
P = 38;
end;
vcov = [
10
];
order = 7;
dynare++/tests/m_0_3_0_0_0_0_0_0.mod deleted 100644 → 0
View file @ 219d2bb3
var lambda a1 a2 k1 k2 i1 i2 c1 c2 l1 l2;
varexo e e1 e2;
parameters beta delta rho sigma phi AA alpha gamma mu chi b Le tau;
beta = 0.99;
delta = 0.025;
rho = 0.95;
sigma = 0.001;
phi = 0.5;
alpha = 0.36;
AA = 0.028058361;
tau = 6.36522e-11;
gamma = 0.25;
mu = -0.2;
chi = 0.83;
b = 2.16872693993;
Le = 2.5;
model;
log(a1) = rho*log(a1(-1))+sigma*(e+e1);
log(a2) = rho*log(a2(-1))+sigma*(e+e2);
lambda = tau*c1^(-1/chi)*(c1^(1-1/chi)+b*(Le-l1)^(1-1/chi))^((1-1/gamma)/(1-1/chi)-1);
lambda = tau*c2^(-1/chi)*(c2^(1-1/chi)+b*(Le-l2)^(1-1/chi))^((1-1/gamma)/(1-1/chi)-1);
tau*(-b)*(Le-l1)^(-1/chi)*(c1^(1-1/chi)+b*(Le-l1)^(1-1/chi))^((1-1/gamma)/(1-1/chi)-1) = -lambda*a1*AA*(1-alpha)*l1^(mu-1)*(alpha*k1(-1)^mu+(1-alpha)*l1^mu)^(1/mu-1);
tau*(-b)*(Le-l2)^(-1/chi)*(c2^(1-1/chi)+b*(Le-l2)^(1-1/chi))^((1-1/gamma)/(1-1/chi)-1) = -lambda*a2*AA*(1-alpha)*l2^(mu-1)*(alpha*k2(-1)^mu+(1-alpha)*l2^mu)^(1/mu-1);
lambda*(1+phi*(i1/k1(-1)-delta)) =beta*lambda(+1)*(1+a1(+1)*AA*alpha*k1^(mu-1)*(alpha*k1^mu+(1-alpha)*l1(+1)^mu)^(1/mu-1)+phi*(1-delta+i1(+1)/k1-0.5*(i1(+1)/k1-delta))*(i1(+1)/k1-delta));
lambda*(1+phi*(i2/k2(-1)-delta)) =beta*lambda(+1)*(1+a2(+1)*AA*alpha*k2^(mu-1)*(alpha*k2^mu+(1-alpha)*l2(+1)^mu)^(1/mu-1)+phi*(1-delta+i2(+1)/k2-0.5*(i2(+1)/k2-delta))*(i2(+1)/k2-delta));
k1 = i1 + (1-delta)*k1(-1);
k2 = i2 + (1-delta)*k2(-1);
c1+i1-delta*k1(-1) + c2+i2-delta*k2(-1) = a1*AA*(alpha*k1(-1)^mu+(1-alpha)*l1^mu)^(1/mu)-(phi/2)*k1(-1)*(i1/k1(-1)-delta)^2 + a2*AA*(alpha*k2(-1)^mu+(1-alpha)*l2^mu)^(1/mu)-(phi/2)*k2(-1)*(i2/k2(-1)-delta)^2;
end;
initval;
a1 = 1;
a2 = 1;
k1 = 1;
k2 = 1;
c1 = 0.028058361;
c2 = 0.028058361;
i1 = 0.025;
i2 = 0.025;
l1 = 1;
l2 = 1;
lambda = 1;
end;
vcov = [1 0 0; 0 1 0; 0 0 1];
order = 4;
\ No newline at end of file
dynare++/tests/m_1_3_0_0_0_0_0_0.mod deleted 100644 → 0
View file @ 219d2bb3
var lambda a1 a2 k1 k2 i1 i2 c1 c2 l1 l2;
varexo e e1 e2;
parameters beta delta rho sigma phi AA alpha Le gamma1 gamma2 mu1 mu2 chi1 chi2 b1 b2 tau1 tau2;
beta = 0.99;
delta = 0.025;
rho = 0.95;
sigma = 0.01;
phi = 0.5;
alpha = 0.36;
AA = 0.028058361;
tau1 = 1.0604611e-11;
tau2 = 2.9305887e-08;
Le = 2.5;
gamma1 = 0.2;
gamma2 = 0.4;
chi1 = 0.75;
chi2 = 0.9;
mu1 = -0.3;
mu2 = 0.3;
b1 = 3.6164368;
b2 = 1.4937381;
model;
log(a1) = rho*log(a1(-1))+sigma*(e+e1);
log(a2) = rho*log(a2(-1))+sigma*(e+e2);
lambda = tau1*c1^(-1/chi1)*(c1^(1-1/chi1)+b1*(Le-l1)^(1-1/chi1))^((1-1/gamma1)/(1-1/chi1)-1);
lambda = tau2*c2^(-1/chi2)*(c2^(1-1/chi2)+b2*(Le-l2)^(1-1/chi2))^((1-1/gamma2)/(1-1/chi2)-1);
tau1*(-b1)*(Le-l1)^(-1/chi1)*(c1^(1-1/chi1)+b1*(Le-l1)^(1-1/chi1))^((1-1/gamma1)/(1-1/chi1)-1) = -lambda*a1*AA*(1-alpha)*l1^(mu1-1)*(alpha*k1(-1)^mu1+(1-alpha)*l1^mu1)^(1/mu1-1);
tau2*(-b2)*(Le-l2)^(-1/chi2)*(c2^(1-1/chi2)+b2*(Le-l2)^(1-1/chi2))^((1-1/gamma2)/(1-1/chi2)-1) = -lambda*a2*AA*(1-alpha)*l2^(mu2-1)*(alpha*k2(-1)^mu2+(1-alpha)*l2^mu2)^(1/mu2-1);
lambda*(1+phi*(i1/k1(-1)-delta)) =beta*lambda(+1)*(1+a1(+1)*AA*alpha*k1^(mu1-1)*(alpha*k1^mu1+(1-alpha)*l1(+1)^mu1)^(1/mu1-1)+phi*(1-delta+i1(+1)/k1-0.5*(i1(+1)/k1-delta))*(i1(+1)/k1-delta));
lambda*(1+phi*(i2/k2(-1)-delta)) =beta*lambda(+1)*(1+a2(+1)*AA*alpha*k2^(mu2-1)*(alpha*k2^mu2+(1-alpha)*l2(+1)^mu2)^(1/mu2-1)+phi*(1-delta+i2(+1)/k2-0.5*(i2(+1)/k2-delta))*(i2(+1)/k2-delta));
k1 = i1 + (1-delta)*k1(-1);
k2 = i2 + (1-delta)*k2(-1);
c1+i1-delta*k1(-1) + c2+i2-delta*k2(-1) = a1*AA*(alpha*k1(-1)^mu1+(1-alpha)*l1^mu1)^(1/mu1)-(phi/2)*k1(-1)*(i1/k1(-1)-delta)^2 + a2*AA*(alpha*k2(-1)^mu2+(1-alpha)*l2^mu2)^(1/mu2)-(phi/2)*k2(-1)*(i2/k2(-1)-delta)^2;
end;
initval;
a1 = 1;
a2 = 1;
k1 = 1;
k2 = 1;
c1 = 0.028058361;
c2 = 0.028058361;
i1 = 0.025;
i2 = 0.025;
l1 = 1;
l2 = 1;
lambda = 1;
end;
vcov = [1 0 0; 0 1 0; 0 0 1];
order = 4;
dynare++/tests/m_1_3_0_0_0_0_0_1.mod deleted 100644 → 0
View file @ 219d2bb3
var lambda a1 a2 a3 a4 k1 k2 k3 k4 i1 i2 i3 i4 c1 c2 c3 c4 l1 l2 l3 l4;
varexo e e1 e2 e3 e4;
parameters beta delta rho sigma phi AA alpha Le gamma1 gamma2 gamma3 gamma4 mu1 mu2 mu3 mu4 chi1 chi2 chi3 chi4 b1 b2 b3 b4 tau1 tau2 tau3 tau4;
beta = 0.99;
delta = 0.025;
rho = 0.95;
sigma = 0.001;
phi = 0.5;
alpha = 0.36;
AA = 0.028058361;
tau1 = 1.0604611e-11;
tau2 = 1.8099765e-09;
tau3 = 2.1096359e-08;
tau4 = 2.9305887e-08;
Le = 2.5;
gamma1 = 0.2;
gamma2 = 0.266666666667;
gamma3 = 0.333333333333;
gamma4 = 0.4;
chi1 = 0.75;
chi2 = 0.8;
chi3 = 0.85;
chi4 = 0.9;
mu1 = -0.3;
mu2 = -0.1;
mu3 = 0.1;
mu4 = 0.3;
b1 = 3.6164368;
b2 = 2.5958433;
b3 = 1.9373921;
b4 = 1.4937381;
model;
log(a1) = rho*log(a1(-1))+sigma*(e+e1);
log(a2) = rho*log(a2(-1))+sigma*(e+e2);
log(a3) = rho*log(a3(-1))+sigma*(e+e3);
log(a4) = rho*log(a4(-1))+sigma*(e+e4);
lambda = tau1*c1^(-1/chi1)*(c1^(1-1/chi1)+b1*(Le-l1)^(1-1/chi1))^((1-1/gamma1)/(1-1/chi1)-1);
lambda = tau2*c2^(-1/chi2)*(c2^(1-1/chi2)+b2*(Le-l2)^(1-1/chi2))^((1-1/gamma2)/(1-1/chi2)-1);
lambda = tau3*c3^(-1/chi3)*(c3^(1-1/chi3)+b3*(Le-l3)^(1-1/chi3))^((1-1/gamma3)/(1-1/chi3)-1);
lambda = tau4*c4^(-1/chi4)*(c4^(1-1/chi4)+b4*(Le-l4)^(1-1/chi4))^((1-1/gamma4)/(1-1/chi4)-1);
tau1*(-b1)*(Le-l1)^(-1/chi1)*(c1^(1-1/chi1)+b1*(Le-l1)^(1-1/chi1))^((1-1/gamma1)/(1-1/chi1)-1) = -lambda*a1*AA*(1-alpha)*l1^(mu1-1)*(alpha*k1(-1)^mu1+(1-alpha)*l1^mu1)^(1/mu1-1);
tau2*(-b2)*(Le-l2)^(-1/chi2)*(c2^(1-1/chi2)+b2*(Le-l2)^(1-1/chi2))^((1-1/gamma2)/(1-1/chi2)-1) = -lambda*a2*AA*(1-alpha)*l2^(mu2-1)*(alpha*k2(-1)^mu2+(1-alpha)*l2^mu2)^(1/mu2-1);
tau3*(-b3)*(Le-l3)^(-1/chi3)*(c3^(1-1/chi3)+b3*(Le-l3)^(1-1/chi3))^((1-1/gamma3)/(1-1/chi3)-1) = -lambda*a3*AA*(1-alpha)*l3^(mu3-1)*(alpha*k3(-1)^mu3+(1-alpha)*l3^mu3)^(1/mu3-1);
tau4*(-b4)*(Le-l4)^(-1/chi4)*(c4^(1-1/chi4)+b4*(Le-l4)^(1-1/chi4))^((1-1/gamma4)/(1-1/chi4)-1) = -lambda*a4*AA*(1-alpha)*l4^(mu4-1)*(alpha*k4(-1)^mu4+(1-alpha)*l4^mu4)^(1/mu4-1);
lambda*(1+phi*(i1/k1(-1)-delta)) =beta*lambda(+1)*(1+a1(+1)*AA*alpha*k1^(mu1-1)*(alpha*k1^mu1+(1-alpha)*l1(+1)^mu1)^(1/mu1-1)+phi*(1-delta+i1(+1)/k1-0.5*(i1(+1)/k1-delta))*(i1(+1)/k1-delta));
lambda*(1+phi*(i2/k2(-1)-delta)) =beta*lambda(+1)*(1+a2(+1)*AA*alpha*k2^(mu2-1)*(alpha*k2^mu2+(1-alpha)*l2(+1)^mu2)^(1/mu2-1)+phi*(1-delta+i2(+1)/k2-0.5*(i2(+1)/k2-delta))*(i2(+1)/k2-delta));
lambda*(1+phi*(i3/k3(-1)-delta)) =beta*lambda(+1)*(1+a3(+1)*AA*alpha*k3^(mu3-1)*(alpha*k3^mu3+(1-alpha)*l3(+1)^mu3)^(1/mu3-1)+phi*(1-delta+i3(+1)/k3-0.5*(i3(+1)/k3-delta))*(i3(+1)/k3-delta));
lambda*(1+phi*(i4/k4(-1)-delta)) =beta*lambda(+1)*(1+a4(+1)*AA*alpha*k4^(mu4-1)*(alpha*k4^mu4+(1-alpha)*l4(+1)^mu4)^(1/mu4-1)+phi*(1-delta+i4(+1)/k4-0.5*(i4(+1)/k4-delta))*(i4(+1)/k4-delta));
k1 = i1 + (1-delta)*k1(-1);
k2 = i2 + (1-delta)*k2(-1);
k3 = i3 + (1-delta)*k3(-1);
k4 = i4 + (1-delta)*k4(-1);
c1+i1-delta*k1(-1) + c2+i2-delta*k2(-1) + c3+i3-delta*k3(-1) + c4+i4-delta*k4(-1) = a1*AA*(alpha*k1(-1)^mu1+(1-alpha)*l1^mu1)^(1/mu1)-(phi/2)*k1(-1)*(i1/k1(-1)-delta)^2 + a2*AA*(alpha*k2(-1)^mu2+(1-alpha)*l2^mu2)^(1/mu2)-(phi/2)*k2(-1)*(i2/k2(-1)-delta)^2 + a3*AA*(alpha*k3(-1)^mu3+(1-alpha)*l3^mu3)^(1/mu3)-(phi/2)*k3(-1)*(i3/k3(-1)-delta)^2 + a4*AA*(alpha*k4(-1)^mu4+(1-alpha)*l4^mu4)^(1/mu4)-(phi/2)*k4(-1)*(i4/k4(-1)-delta)^2;
end;
initval;
a1 = 1;
a2 = 1;
a3 = 1;
a4 = 1;
k1 = 1;
k2 = 1;
k3 = 1;
k4 = 1;
c1 = 0.028058361;
c2 = 0.028058361;
c3 = 0.028058361;
c4 = 0.028058361;
i1 = 0.025;
i2 = 0.025;
i3 = 0.025;
i4 = 0.025;
l1 = 1;
l2 = 1;
l3 = 1;
l4 = 1;
lambda = 1;
end;
vcov = [1.0 0 0 0 0; 0 1.0 0 0 0; 0 0 1.0 0 0; 0 0 0 1.0 0; 0 0 0 0 1.0];
order = 4;
dynare++/tests/or0a.mod deleted 100644 → 0
View file @ 219d2bb3
var
C
CF
CF_STAR
CH
CH_STAR
CN
CN_STAR
CT
CT_STAR
C_STAR
E
KE
KE_STAR
L
L_STAR
P
PF
PF_STAR
PH
PH_STAR
PN
PN_STAR
PT
PT_STAR
P_STAR
W
W_STAR
Y
Y_STAR
;
varexo k k_star m m_star;
parameters epsi chi thet nu phi gam;
epsi = 0.5;
nu = 3;
chi = 1.2;
phi = 4;
thet = 3;
gam = 0.5;
model;
C = (1/chi)*(exp(m)/P)^epsi;
C_STAR = (1/chi)*(exp(m_star)/P_STAR)^epsi;
CN = (1-gam)*(P/PN)*C;
CN_STAR = (1-gam)*(P_STAR/PN_STAR)*C_STAR;
CT = gam*(P/PT)*C;
CT_STAR = gam*(P_STAR/PT_STAR)*C_STAR;
CH = 0.5*(PT/PH)*CT;
CH_STAR = 0.5*(PT_STAR/PH_STAR)*CT_STAR;
CF = 0.5*(PT/PF)*CT;
CF_STAR = 0.5*(PT_STAR/PF_STAR)*CT_STAR;
P = PT^gam*PN^(1-gam);
P_STAR = PT_STAR^gam*PN_STAR^(1-gam);
PT = sqrt(PH*PF);
PT_STAR = sqrt(PH_STAR*PF_STAR);
PH = (thet/(thet-1))*W(-1);
PF_STAR = (thet/(thet-1))*W_STAR(-1);
PN = PH;
PN_STAR = PF_STAR;
L = Y;
L_STAR = Y_STAR;
(L(+1)/(P(+1)*C(+1)))*W = (phi/(phi-1))*KE(+1)*L(+1)^nu;
(L_STAR(+1)/(P_STAR(+1)*C_STAR(+1)))*W_STAR = (phi/(phi-1))*KE_STAR(+1)*L_STAR(+1)^nu;
P*C = Y*PH;
P_STAR*C_STAR = Y_STAR*PF_STAR;
Y = CH + CH_STAR + CN;
Y_STAR = CF + CF_STAR + CN_STAR;
PT = E*PT_STAR;
KE = exp(k);
KE_STAR = exp(k_star);
end;
initval;
C = 1;
PH = 1;
P = 1;
PN = 1;
PT = 1;
L = 1;
Y = 1;
W = 1;
CF = 0.25;
CH = 0.25;
CT = 0.5;
CN = 0.5;
PF = 1;
C_STAR = 1;
PH_STAR = 1;
P_STAR = 1;
PN_STAR = 1;
PT_STAR = 1;
L_STAR = 1;
Y_STAR = 1;
W_STAR = 1;
CF_STAR = 0.25;
CH_STAR = 0.25;
CT_STAR = 0.5;
CN_STAR = 0.5;
PF_STAR = 1;
KE = 1;
KE_STAR = 1;
E = 1;
k = 0;
k_star = 0;
m = 0;
m_star = 0;
end;
vcov = [
0.01 0 -0.01 0;
0 0.01 0 -0.01;
-0.01 0 0.01 0;
0 -0.01 0 0.01
];
order=4;
dynare++/tests/portfolio.mod deleted 100644 → 0
View file @ 219d2bb3
var DOTQ Q1 Q2 X1 X2 C D1 D2;
varexo E_D1 E_D2;
parameters beta, r1, r2, gamma, d, rho1, rho2;
beta = 0.95;
r1 = 0.2;
r2 = 0.05;
gamma = 0.78;
d = 0.10;
rho1 = 0.8;
rho2 = 0.2;
model;
C + X1 + X2 = D1*Q1 + D2*Q2;
Q1+Q2 = 1;
C^(-gamma)/(1-2*r1*X1) = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r1*X1(+1))*(D1(+1)*(1-2*r1*X1(+1))+1);
C^(-gamma)/(1-2*r2*X2) = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r2*X2(+1))*(D2(+1)*(1-2*r2*X2(+1))+1);
DOTQ*Q1 = Q1(-1) + X1(-1) - r1*X1(-1)^2;
DOTQ*Q2 = Q2(-1) + X2(-1) - r2*X2(-1)^2;
D1/d = D1(-1)^rho1/(d^rho1)*exp(E_D1);
D2/d = D2(-1)^rho2/(d^rho2)*exp(E_D2);
/*
D1-d = rho1*(D1(-1)-d) + E_D1;
D2-d = rho2*(D2(-1)-d) + E_D2;
*/
end;
initval;
C =0.0441234;
D1 =0.1000000000000;
D2 =0.1000000000000;
DOTQ =1.05567;
Q1 =0.333333;
Q2 =0.666667;
X1 =0.0186255;
X2 =0.0372511;
end;
vcov = [
0.04 0;
0 0.01
];
order=5;
dynare++/tests/portfolio4.mod deleted 100644 → 0
View file @ 219d2bb3
var DOTQ Q1 Q2 Q3 Q4 X1 X2 X3 X4 C D1 D2 D3 D4 V;
varexo E_D1 E_D2 E_D3 E_D4;
parameters beta, r1, r2, r3, r4, gamma, ed1, ed2, ed3, ed4, rho1, rho2, rho3, rho4;
beta = 0.95;
r1 = 0.2;
r2 = 0.1;
r3 = 0.06;
r4 = 0.03;
gamma = 0.7;
ed1 = 0.1;
ed2 = 0.1;
ed3 = 0.1;
ed4 = 0.1;
rho1 = 0.3;
rho2 = 0.01;
rho3 = 0.6;
rho4 = 0.6;
model;
Q1+Q2+Q3+Q4 = 1;
C + X1 + X2 + X3 + X4 = D1*Q1 + D2*Q2 + D3*Q3 + D4*Q4;
DOTQ*Q1 = Q1(-1) + X1(-1) - r1*X1(-1)*X1(-1);
DOTQ*Q2 = Q2(-1) + X2(-1) - r2*X2(-1)*X2(-1);
DOTQ*Q3 = Q3(-1) + X3(-1) - r3*X3(-1)*X3(-1);
DOTQ*Q4 = Q4(-1) + X4(-1) - r4*X4(-1)*X4(-1);
C^(-gamma)/(1-2*r1*X1) = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r1*X1(+1))*(D1(+1)*(1-2*r1*X1(+1))+1);
C^(-gamma)/(1-2*r2*X2) = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r2*X2(+1))*(D2(+1)*(1-2*r2*X2(+1))+1);
C^(-gamma)/(1-2*r3*X3) = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r3*X3(+1))*(D3(+1)*(1-2*r3*X3(+1))+1);
C^(-gamma)/(1-2*r4*X4) = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r4*X4(+1))*(D4(+1)*(1-2*r4*X4(+1))+1);
V = C^(1-gamma)/(1-gamma) + beta*V(+1);
D1/ed1 = D1(-1)^rho1/(ed1^rho1)*exp(E_D1);
D2/ed2 = D2(-1)^rho2/(ed2^rho2)*exp(E_D2);
D3/ed3 = D3(-1)^rho3/(ed3^rho3)*exp(E_D3);
D4/ed4 = D4(-1)^rho4/(ed4^rho4)*exp(E_D4);
/*
D1-ed1 = rho1*(D1(-1)-ed1) + E_D1;
D2-ed2 = rho2*(D2(-1)-ed2) + E_D2;
D3-ed3 = rho3*(D3(-1)-ed3) + E_D3;
D4-ed4 = rho4*(D4(-1)-ed4) + E_D4;
*/
end;
initval;
D1 = ed1;
D2 = ed2;
D3 = ed3;
D4 = ed4;
DOTQ = 1.05; // provide a guess not larger than the true value
X1 = 1/2/r1*(1-1/D1*(1/beta*DOTQ^gamma - 1));
X2 = 1/2/r2*(1-1/D2*(1/beta*DOTQ^gamma - 1));
X3 = 1/2/r3*(1-1/D3*(1/beta*DOTQ^gamma - 1));
X4 = 1/2/r4*(1-1/D4*(1/beta*DOTQ^gamma - 1));
Q1 = 1/(DOTQ-1)*(X1 - r1*X1*X1);
Q2 = 1/(DOTQ-1)*(X2 - r2*X2*X2);
Q3 = 1/(DOTQ-1)*(X3 - r3*X3*X3);
Q4 = 1/(DOTQ-1)*(X4 - r4*X4*X4);
C = ed1*Q1 + ed2*Q2 + ed3*Q3 + ed4*Q4 - X1 - X2 - X3 - X4;
V = 1/(1-beta)*C^(1-gamma)/(1-gamma);
end;
/*
vcov = [
0.0005 0 0 0;
0 0.00025 0 0;
0 0 0.0005 0;
0 0 0 0.00025
];
*/
vcov = [
0.05 0 0 0;
0 0.025 0 0;
0 0 0.05 0;
0 0 0 0.025
];
order=5;
dynare++/tests/portfolio4_norm.mod deleted 100644 → 0
View file @ 219d2bb3
var DOTQ Q1 Q2 Q3 Q4 X1 X2 X3 X4 C D1 D2 D3 D4 V;
varexo E_D1 E_D2 E_D3 E_D4;
parameters beta, r1, r2, r3, r4, gamma, ed1, ed2, ed3, ed4, rho1, rho2, rho3, rho4;
beta = 0.95;
r1 = 0.2;
r2 = 0.1;
r3 = 0.06;
r4 = 0.03;
gamma = 0.7;
ed1 = 0.1;
ed2 = 0.1;
ed3 = 0.1;
ed4 = 0.1;
rho1 = 0.3;
rho2 = 0.01;
rho3 = 0.6;
rho4 = 0.6;
model;
1 = (C + X1 + X2 + X3 + X4)/ (D1*Q1 + D2*Q2 + D3*Q3 + D4*Q4);
1 = (Q1(-1) + X1(-1) - r1*X1(-1)*X1(-1))/(DOTQ*Q1);
1 = (Q2(-1) + X2(-1) - r2*X2(-1)*X2(-1))/(DOTQ*Q2);
1 = (Q3(-1) + X3(-1) - r3*X3(-1)*X3(-1))/(DOTQ*Q3);
1 = (Q4(-1) + X4(-1) - r4*X4(-1)*X4(-1))/(DOTQ*Q4);
Q1+Q2+Q3+Q4 = 1;
1 = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r1*X1(+1))*(D1(+1)*(1-2*r1*X1(+1))+1)/(C^(-gamma)/(1-2*r1*X1));
1 = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r2*X2(+1))*(D2(+1)*(1-2*r2*X2(+1))+1)/(C^(-gamma)/(1-2*r2*X2));
1 = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r3*X3(+1))*(D3(+1)*(1-2*r3*X3(+1))+1)/(C^(-gamma)/(1-2*r3*X3));
1 = beta*DOTQ(+1)^(-gamma)*C(+1)^(-gamma)/(1-2*r4*X4(+1))*(D4(+1)*(1-2*r4*X4(+1))+1)/(C^(-gamma)/(1-2*r4*X4));
1 = D1(-1)^rho1/D1*ed1^(1-rho1)*exp(E_D1);
1 = D2(-1)^rho2/D2*ed2^(1-rho2)*exp(E_D2);
1 = D3(-1)^rho3/D3*ed3^(1-rho3)*exp(E_D3);
1 = D4(-1)^rho4/D4*ed4^(1-rho4)*exp(E_D4);
/*
D1-ed1 = rho1*(D1(-1)-ed1) + E_D1;
D2-ed2 = rho2*(D2(-1)-ed2) + E_D2;
D3-ed3 = rho3*(D3(-1)-ed3) + E_D3;
D4-ed4 = rho4*(D4(-1)-ed4) + E_D4;
*/
V/(C^(1-gamma)/(1-gamma) + beta*V(+1)) = 1;
end;
initval;
Q1 =0.0769231;
Q2 =0.1538462;
Q3 =0.2564103;
Q4 =0.5128205;
X1 =0.0049761;
X2 =0.0099522;
X3 =0.0165871;
X4 =0.0331741;
D1 =0.1;
D2 =0.1;
D3 =0.1;
D4 =0.1;
DOTQ =1.0646251;
C =0.0353105;
V =24.450057;
end;
/*
vcov = [
0.0005 0 0 0;
0 0.00025 0 0;
0 0 0.0005 0;
0 0 0 0.00025
];
*/
vcov = [
0.05 0 0 0;
0 0.025 0 0;
0 0 0.05 0;
0 0 0 0.025
];
order=4;
dynare++/tests/psd_exo3.mod deleted 100644 → 0
View file @ 219d2bb3
var c k x;
varexo rho;
parameters a alph gam bet lamb;
alph = 0.7;
bet = 0.95;
gam = 2;
a = 1.052632;
lamb = 0.9;
model;
c^(-gam) = bet*c(+1)^(-gam)*a*exp(x(+1))*k^(-alph);
k = a*exp(x)*k(-1)^(1-alph)/(1-alph)-c;
x = lamb*x(-1)+rho;
end;
initval;
k = 1;
c = 2.508;
x = 0;
rho = 0;
end;
vcov=[0.0001];
order=6;
\ No newline at end of file
dynare++/tests/q3a2.mod deleted 100644 → 0
View file @ 219d2bb3
var c,a1,k1,a2,k2;
varexo e,e1,e2;
parameters beta, gamma, phi, delta, alpha, rho, zeta, sigma, N;
alpha = 0.36;
delta = 0.025;
phi = 2;
gamma = 1;
beta = 0.99;
rho = 0.95;
zeta = 2;
sigma = 0.005;
N = 2;
model;
c^(-gamma)*(1+phi*zeta*(k1-k1(-1))^(zeta-1)/(2*k1(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k1-k1(-1))^(zeta-1)/(2*k1(-1))+phi*(k1-k1(-1))^zeta/(2*k1(-1)^2)-delta+alpha*a1(+1)*k1^(alpha-1));
log(a1) = rho*log(a1(-1))+sigma*(e+e1);
c^(-gamma)*(1+phi*zeta*(k2-k2(-1))^(zeta-1)/(2*k2(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k2-k2(-1))^(zeta-1)/(2*k2(-1))+phi*(k2-k2(-1))^zeta/(2*k2(-1)^2)-delta+alpha*a2(+1)*k2^(alpha-1));
log(a2) = rho*log(a2(-1))+sigma*(e+e2);
N*c+k1+phi*(k1-k1(-1))^zeta/(2*k1)-(1-delta)*k1(-1)+k2+phi*(k2-k2(-1))^zeta/(2*k2)-(1-delta)*k2(-1)= a1*k1(-1)^alpha+a2*k2(-1)^alpha;
end;
initval;
c=1;
e=0;
a1=1;
e1=0;
k1=10;
a2=1;
e2=0;
k2=10;
end;
vcov=[1 0 0; 0 1 0; 0 0 1];
order = 4;
\ No newline at end of file
dynare++/tests/q3a50.mod deleted 100644 → 0
View file @ 219d2bb3
periods 5000;
var c,a1,k1,a2,k2,a3,k3,a4,k4,a5,k5,a6,k6,a7,k7,a8,k8,a9,k9,a10,k10,a11,k11,a12,k12,a13,k13,a14,k14,a15,k15,a16,k16,a17,k17,a18,k18,a19,k19,a20,k20,a21,k21,a22,k22,a23,k23,a24,k24,a25,k25,a26,k26,a27,k27,a28,k28,a29,k29,a30,k30,a31,k31,a32,k32,a33,k33,a34,k34,a35,k35,a36,k36,a37,k37,a38,k38,a39,k39,a40,k40,a41,k41,a42,k42,a43,k43,a44,k44,a45,k45,a46,k46,a47,k47,a48,k48,a49,k49,a50,k50;
varexo e,e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19,e20,e21,e22,e23,e24,e25,e26,e27,e28,e29,e30,e31,e32,e33,e34,e35,e36,e37,e38,e39,e40,e41,e42,e43,e44,e45,e46,e47,e48,e49,e50;
parameters beta, gamma, phi, delta, alpha, rho, zeta, sigma, N;
alpha = 0.36;
delta = 0.1;
phi = 2;
gamma = 1;
beta = 0.99;
rho = 0.95;
zeta = 2;
sigma = 2;
N = 50;
model;
c^(-gamma)*(1+phi*zeta*(k1-k1(-1))^(zeta-1)/(2*k1(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k1-k1(-1))^(zeta-1)/(2*k1(-1))+phi*(k1-k1(-1))^zeta/(2*k1(-1)^2)-delta+alpha*a1(+1)*k1^(alpha-1));
log(a1) = rho*log(a1(-1))+sigma*(e+e1);
c^(-gamma)*(1+phi*zeta*(k2-k2(-1))^(zeta-1)/(2*k2(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k2-k2(-1))^(zeta-1)/(2*k2(-1))+phi*(k2-k2(-1))^zeta/(2*k2(-1)^2)-delta+alpha*a2(+1)*k2^(alpha-1));
log(a2) = rho*log(a2(-1))+sigma*(e+e2);
c^(-gamma)*(1+phi*zeta*(k3-k3(-1))^(zeta-1)/(2*k3(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k3-k3(-1))^(zeta-1)/(2*k3(-1))+phi*(k3-k3(-1))^zeta/(2*k3(-1)^2)-delta+alpha*a3(+1)*k3^(alpha-1));
log(a3) = rho*log(a3(-1))+sigma*(e+e3);
c^(-gamma)*(1+phi*zeta*(k4-k4(-1))^(zeta-1)/(2*k4(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k4-k4(-1))^(zeta-1)/(2*k4(-1))+phi*(k4-k4(-1))^zeta/(2*k4(-1)^2)-delta+alpha*a4(+1)*k4^(alpha-1));
log(a4) = rho*log(a4(-1))+sigma*(e+e4);
c^(-gamma)*(1+phi*zeta*(k5-k5(-1))^(zeta-1)/(2*k5(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k5-k5(-1))^(zeta-1)/(2*k5(-1))+phi*(k5-k5(-1))^zeta/(2*k5(-1)^2)-delta+alpha*a5(+1)*k5^(alpha-1));
log(a5) = rho*log(a5(-1))+sigma*(e+e5);
c^(-gamma)*(1+phi*zeta*(k6-k6(-1))^(zeta-1)/(2*k6(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k6-k6(-1))^(zeta-1)/(2*k6(-1))+phi*(k6-k6(-1))^zeta/(2*k6(-1)^2)-delta+alpha*a6(+1)*k6^(alpha-1));
log(a6) = rho*log(a6(-1))+sigma*(e+e6);
c^(-gamma)*(1+phi*zeta*(k7-k7(-1))^(zeta-1)/(2*k7(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k7-k7(-1))^(zeta-1)/(2*k7(-1))+phi*(k7-k7(-1))^zeta/(2*k7(-1)^2)-delta+alpha*a7(+1)*k7^(alpha-1));
log(a7) = rho*log(a7(-1))+sigma*(e+e7);
c^(-gamma)*(1+phi*zeta*(k8-k8(-1))^(zeta-1)/(2*k8(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k8-k8(-1))^(zeta-1)/(2*k8(-1))+phi*(k8-k8(-1))^zeta/(2*k8(-1)^2)-delta+alpha*a8(+1)*k8^(alpha-1));
log(a8) = rho*log(a8(-1))+sigma*(e+e8);
c^(-gamma)*(1+phi*zeta*(k9-k9(-1))^(zeta-1)/(2*k9(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k9-k9(-1))^(zeta-1)/(2*k9(-1))+phi*(k9-k9(-1))^zeta/(2*k9(-1)^2)-delta+alpha*a9(+1)*k9^(alpha-1));
log(a9) = rho*log(a9(-1))+sigma*(e+e9);
c^(-gamma)*(1+phi*zeta*(k10-k10(-1))^(zeta-1)/(2*k10(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k10-k10(-1))^(zeta-1)/(2*k10(-1))+phi*(k10-k10(-1))^zeta/(2*k10(-1)^2)-delta+alpha*a10(+1)*k10^(alpha-1));
log(a10) = rho*log(a10(-1))+sigma*(e+e10);
c^(-gamma)*(1+phi*zeta*(k11-k11(-1))^(zeta-1)/(2*k11(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k11-k11(-1))^(zeta-1)/(2*k11(-1))+phi*(k11-k11(-1))^zeta/(2*k11(-1)^2)-delta+alpha*a11(+1)*k11^(alpha-1));
log(a11) = rho*log(a11(-1))+sigma*(e+e11);
c^(-gamma)*(1+phi*zeta*(k12-k12(-1))^(zeta-1)/(2*k12(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k12-k12(-1))^(zeta-1)/(2*k12(-1))+phi*(k12-k12(-1))^zeta/(2*k12(-1)^2)-delta+alpha*a12(+1)*k12^(alpha-1));
log(a12) = rho*log(a12(-1))+sigma*(e+e12);
c^(-gamma)*(1+phi*zeta*(k13-k13(-1))^(zeta-1)/(2*k13(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k13-k13(-1))^(zeta-1)/(2*k13(-1))+phi*(k13-k13(-1))^zeta/(2*k13(-1)^2)-delta+alpha*a13(+1)*k13^(alpha-1));
log(a13) = rho*log(a13(-1))+sigma*(e+e13);
c^(-gamma)*(1+phi*zeta*(k14-k14(-1))^(zeta-1)/(2*k14(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k14-k14(-1))^(zeta-1)/(2*k14(-1))+phi*(k14-k14(-1))^zeta/(2*k14(-1)^2)-delta+alpha*a14(+1)*k14^(alpha-1));
log(a14) = rho*log(a14(-1))+sigma*(e+e14);
c^(-gamma)*(1+phi*zeta*(k15-k15(-1))^(zeta-1)/(2*k15(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k15-k15(-1))^(zeta-1)/(2*k15(-1))+phi*(k15-k15(-1))^zeta/(2*k15(-1)^2)-delta+alpha*a15(+1)*k15^(alpha-1));
log(a15) = rho*log(a15(-1))+sigma*(e+e15);
c^(-gamma)*(1+phi*zeta*(k16-k16(-1))^(zeta-1)/(2*k16(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k16-k16(-1))^(zeta-1)/(2*k16(-1))+phi*(k16-k16(-1))^zeta/(2*k16(-1)^2)-delta+alpha*a16(+1)*k16^(alpha-1));
log(a16) = rho*log(a16(-1))+sigma*(e+e16);
c^(-gamma)*(1+phi*zeta*(k17-k17(-1))^(zeta-1)/(2*k17(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k17-k17(-1))^(zeta-1)/(2*k17(-1))+phi*(k17-k17(-1))^zeta/(2*k17(-1)^2)-delta+alpha*a17(+1)*k17^(alpha-1));
log(a17) = rho*log(a17(-1))+sigma*(e+e17);
c^(-gamma)*(1+phi*zeta*(k18-k18(-1))^(zeta-1)/(2*k18(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k18-k18(-1))^(zeta-1)/(2*k18(-1))+phi*(k18-k18(-1))^zeta/(2*k18(-1)^2)-delta+alpha*a18(+1)*k18^(alpha-1));
log(a18) = rho*log(a18(-1))+sigma*(e+e18);
c^(-gamma)*(1+phi*zeta*(k19-k19(-1))^(zeta-1)/(2*k19(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k19-k19(-1))^(zeta-1)/(2*k19(-1))+phi*(k19-k19(-1))^zeta/(2*k19(-1)^2)-delta+alpha*a19(+1)*k19^(alpha-1));
log(a19) = rho*log(a19(-1))+sigma*(e+e19);
c^(-gamma)*(1+phi*zeta*(k20-k20(-1))^(zeta-1)/(2*k20(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k20-k20(-1))^(zeta-1)/(2*k20(-1))+phi*(k20-k20(-1))^zeta/(2*k20(-1)^2)-delta+alpha*a20(+1)*k20^(alpha-1));
log(a20) = rho*log(a20(-1))+sigma*(e+e20);
c^(-gamma)*(1+phi*zeta*(k21-k21(-1))^(zeta-1)/(2*k21(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k21-k21(-1))^(zeta-1)/(2*k21(-1))+phi*(k21-k21(-1))^zeta/(2*k21(-1)^2)-delta+alpha*a21(+1)*k21^(alpha-1));
log(a21) = rho*log(a21(-1))+sigma*(e+e21);
c^(-gamma)*(1+phi*zeta*(k22-k22(-1))^(zeta-1)/(2*k22(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k22-k22(-1))^(zeta-1)/(2*k22(-1))+phi*(k22-k22(-1))^zeta/(2*k22(-1)^2)-delta+alpha*a22(+1)*k22^(alpha-1));
log(a22) = rho*log(a22(-1))+sigma*(e+e22);
c^(-gamma)*(1+phi*zeta*(k23-k23(-1))^(zeta-1)/(2*k23(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k23-k23(-1))^(zeta-1)/(2*k23(-1))+phi*(k23-k23(-1))^zeta/(2*k23(-1)^2)-delta+alpha*a23(+1)*k23^(alpha-1));
log(a23) = rho*log(a23(-1))+sigma*(e+e23);
c^(-gamma)*(1+phi*zeta*(k24-k24(-1))^(zeta-1)/(2*k24(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k24-k24(-1))^(zeta-1)/(2*k24(-1))+phi*(k24-k24(-1))^zeta/(2*k24(-1)^2)-delta+alpha*a24(+1)*k24^(alpha-1));
log(a24) = rho*log(a24(-1))+sigma*(e+e24);
c^(-gamma)*(1+phi*zeta*(k25-k25(-1))^(zeta-1)/(2*k25(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k25-k25(-1))^(zeta-1)/(2*k25(-1))+phi*(k25-k25(-1))^zeta/(2*k25(-1)^2)-delta+alpha*a25(+1)*k25^(alpha-1));
log(a25) = rho*log(a25(-1))+sigma*(e+e25);
c^(-gamma)*(1+phi*zeta*(k26-k26(-1))^(zeta-1)/(2*k26(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k26-k26(-1))^(zeta-1)/(2*k26(-1))+phi*(k26-k26(-1))^zeta/(2*k26(-1)^2)-delta+alpha*a26(+1)*k26^(alpha-1));
log(a26) = rho*log(a26(-1))+sigma*(e+e26);
c^(-gamma)*(1+phi*zeta*(k27-k27(-1))^(zeta-1)/(2*k27(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k27-k27(-1))^(zeta-1)/(2*k27(-1))+phi*(k27-k27(-1))^zeta/(2*k27(-1)^2)-delta+alpha*a27(+1)*k27^(alpha-1));
log(a27) = rho*log(a27(-1))+sigma*(e+e27);
c^(-gamma)*(1+phi*zeta*(k28-k28(-1))^(zeta-1)/(2*k28(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k28-k28(-1))^(zeta-1)/(2*k28(-1))+phi*(k28-k28(-1))^zeta/(2*k28(-1)^2)-delta+alpha*a28(+1)*k28^(alpha-1));
log(a28) = rho*log(a28(-1))+sigma*(e+e28);
c^(-gamma)*(1+phi*zeta*(k29-k29(-1))^(zeta-1)/(2*k29(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k29-k29(-1))^(zeta-1)/(2*k29(-1))+phi*(k29-k29(-1))^zeta/(2*k29(-1)^2)-delta+alpha*a29(+1)*k29^(alpha-1));
log(a29) = rho*log(a29(-1))+sigma*(e+e29);
c^(-gamma)*(1+phi*zeta*(k30-k30(-1))^(zeta-1)/(2*k30(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k30-k30(-1))^(zeta-1)/(2*k30(-1))+phi*(k30-k30(-1))^zeta/(2*k30(-1)^2)-delta+alpha*a30(+1)*k30^(alpha-1));
log(a30) = rho*log(a30(-1))+sigma*(e+e30);
c^(-gamma)*(1+phi*zeta*(k31-k31(-1))^(zeta-1)/(2*k31(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k31-k31(-1))^(zeta-1)/(2*k31(-1))+phi*(k31-k31(-1))^zeta/(2*k31(-1)^2)-delta+alpha*a31(+1)*k31^(alpha-1));
log(a31) = rho*log(a31(-1))+sigma*(e+e31);
c^(-gamma)*(1+phi*zeta*(k32-k32(-1))^(zeta-1)/(2*k32(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k32-k32(-1))^(zeta-1)/(2*k32(-1))+phi*(k32-k32(-1))^zeta/(2*k32(-1)^2)-delta+alpha*a32(+1)*k32^(alpha-1));
log(a32) = rho*log(a32(-1))+sigma*(e+e32);
c^(-gamma)*(1+phi*zeta*(k33-k33(-1))^(zeta-1)/(2*k33(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k33-k33(-1))^(zeta-1)/(2*k33(-1))+phi*(k33-k33(-1))^zeta/(2*k33(-1)^2)-delta+alpha*a33(+1)*k33^(alpha-1));
log(a33) = rho*log(a33(-1))+sigma*(e+e33);
c^(-gamma)*(1+phi*zeta*(k34-k34(-1))^(zeta-1)/(2*k34(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k34-k34(-1))^(zeta-1)/(2*k34(-1))+phi*(k34-k34(-1))^zeta/(2*k34(-1)^2)-delta+alpha*a34(+1)*k34^(alpha-1));
log(a34) = rho*log(a34(-1))+sigma*(e+e34);
c^(-gamma)*(1+phi*zeta*(k35-k35(-1))^(zeta-1)/(2*k35(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k35-k35(-1))^(zeta-1)/(2*k35(-1))+phi*(k35-k35(-1))^zeta/(2*k35(-1)^2)-delta+alpha*a35(+1)*k35^(alpha-1));
log(a35) = rho*log(a35(-1))+sigma*(e+e35);
c^(-gamma)*(1+phi*zeta*(k36-k36(-1))^(zeta-1)/(2*k36(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k36-k36(-1))^(zeta-1)/(2*k36(-1))+phi*(k36-k36(-1))^zeta/(2*k36(-1)^2)-delta+alpha*a36(+1)*k36^(alpha-1));
log(a36) = rho*log(a36(-1))+sigma*(e+e36);
c^(-gamma)*(1+phi*zeta*(k37-k37(-1))^(zeta-1)/(2*k37(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k37-k37(-1))^(zeta-1)/(2*k37(-1))+phi*(k37-k37(-1))^zeta/(2*k37(-1)^2)-delta+alpha*a37(+1)*k37^(alpha-1));
log(a37) = rho*log(a37(-1))+sigma*(e+e37);
c^(-gamma)*(1+phi*zeta*(k38-k38(-1))^(zeta-1)/(2*k38(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k38-k38(-1))^(zeta-1)/(2*k38(-1))+phi*(k38-k38(-1))^zeta/(2*k38(-1)^2)-delta+alpha*a38(+1)*k38^(alpha-1));
log(a38) = rho*log(a38(-1))+sigma*(e+e38);
c^(-gamma)*(1+phi*zeta*(k39-k39(-1))^(zeta-1)/(2*k39(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k39-k39(-1))^(zeta-1)/(2*k39(-1))+phi*(k39-k39(-1))^zeta/(2*k39(-1)^2)-delta+alpha*a39(+1)*k39^(alpha-1));
log(a39) = rho*log(a39(-1))+sigma*(e+e39);
c^(-gamma)*(1+phi*zeta*(k40-k40(-1))^(zeta-1)/(2*k40(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k40-k40(-1))^(zeta-1)/(2*k40(-1))+phi*(k40-k40(-1))^zeta/(2*k40(-1)^2)-delta+alpha*a40(+1)*k40^(alpha-1));
log(a40) = rho*log(a40(-1))+sigma*(e+e40);
c^(-gamma)*(1+phi*zeta*(k41-k41(-1))^(zeta-1)/(2*k41(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k41-k41(-1))^(zeta-1)/(2*k41(-1))+phi*(k41-k41(-1))^zeta/(2*k41(-1)^2)-delta+alpha*a41(+1)*k41^(alpha-1));
log(a41) = rho*log(a41(-1))+sigma*(e+e41);
c^(-gamma)*(1+phi*zeta*(k42-k42(-1))^(zeta-1)/(2*k42(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k42-k42(-1))^(zeta-1)/(2*k42(-1))+phi*(k42-k42(-1))^zeta/(2*k42(-1)^2)-delta+alpha*a42(+1)*k42^(alpha-1));
log(a42) = rho*log(a42(-1))+sigma*(e+e42);
c^(-gamma)*(1+phi*zeta*(k43-k43(-1))^(zeta-1)/(2*k43(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k43-k43(-1))^(zeta-1)/(2*k43(-1))+phi*(k43-k43(-1))^zeta/(2*k43(-1)^2)-delta+alpha*a43(+1)*k43^(alpha-1));
log(a43) = rho*log(a43(-1))+sigma*(e+e43);
c^(-gamma)*(1+phi*zeta*(k44-k44(-1))^(zeta-1)/(2*k44(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k44-k44(-1))^(zeta-1)/(2*k44(-1))+phi*(k44-k44(-1))^zeta/(2*k44(-1)^2)-delta+alpha*a44(+1)*k44^(alpha-1));
log(a44) = rho*log(a44(-1))+sigma*(e+e44);
c^(-gamma)*(1+phi*zeta*(k45-k45(-1))^(zeta-1)/(2*k45(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k45-k45(-1))^(zeta-1)/(2*k45(-1))+phi*(k45-k45(-1))^zeta/(2*k45(-1)^2)-delta+alpha*a45(+1)*k45^(alpha-1));
log(a45) = rho*log(a45(-1))+sigma*(e+e45);
c^(-gamma)*(1+phi*zeta*(k46-k46(-1))^(zeta-1)/(2*k46(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k46-k46(-1))^(zeta-1)/(2*k46(-1))+phi*(k46-k46(-1))^zeta/(2*k46(-1)^2)-delta+alpha*a46(+1)*k46^(alpha-1));
log(a46) = rho*log(a46(-1))+sigma*(e+e46);
c^(-gamma)*(1+phi*zeta*(k47-k47(-1))^(zeta-1)/(2*k47(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k47-k47(-1))^(zeta-1)/(2*k47(-1))+phi*(k47-k47(-1))^zeta/(2*k47(-1)^2)-delta+alpha*a47(+1)*k47^(alpha-1));
log(a47) = rho*log(a47(-1))+sigma*(e+e47);
c^(-gamma)*(1+phi*zeta*(k48-k48(-1))^(zeta-1)/(2*k48(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k48-k48(-1))^(zeta-1)/(2*k48(-1))+phi*(k48-k48(-1))^zeta/(2*k48(-1)^2)-delta+alpha*a48(+1)*k48^(alpha-1));
log(a48) = rho*log(a48(-1))+sigma*(e+e48);
c^(-gamma)*(1+phi*zeta*(k49-k49(-1))^(zeta-1)/(2*k49(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k49-k49(-1))^(zeta-1)/(2*k49(-1))+phi*(k49-k49(-1))^zeta/(2*k49(-1)^2)-delta+alpha*a49(+1)*k49^(alpha-1));
log(a49) = rho*log(a49(-1))+sigma*(e+e49);
c^(-gamma)*(1+phi*zeta*(k50-k50(-1))^(zeta-1)/(2*k50(-1))) = beta*c(+1)^(-gamma)*(1+phi*zeta*(k50-k50(-1))^(zeta-1)/(2*k50(-1))+phi*(k50-k50(-1))^zeta/(2*k50(-1)^2)-delta+alpha*a50(+1)*k50^(alpha-1));
log(a50) = rho*log(a50(-1))+sigma*(e+e50);
N*c+k1+phi*(k1-k1(-1))^zeta/(2*k1)-(1-delta)*k1(-1)+k2+phi*(k2-k2(-1))^zeta/(2*k2)-(1-delta)*k2(-1)+k3+phi*(k3-k3(-1))^zeta/(2*k3)-(1-delta)*k3(-1)+k4+phi*(k4-k4(-1))^zeta/(2*k4)-(1-delta)*k4(-1)+k5+phi*(k5-k5(-1))^zeta/(2*k5)-(1-delta)*k5(-1)+k6+phi*(k6-k6(-1))^zeta/(2*k6)-(1-delta)*k6(-1)+k7+phi*(k7-k7(-1))^zeta/(2*k7)-(1-delta)*k7(-1)+k8+phi*(k8-k8(-1))^zeta/(2*k8)-(1-delta)*k8(-1)+k9+phi*(k9-k9(-1))^zeta/(2*k9)-(1-delta)*k9(-1)+k10+phi*(k10-k10(-1))^zeta/(2*k10)-(1-delta)*k10(-1)+k11+phi*(k11-k11(-1))^zeta/(2*k11)-(1-delta)*k11(-1)+k12+phi*(k12-k12(-1))^zeta/(2*k12)-(1-delta)*k12(-1)+k13+phi*(k13-k13(-1))^zeta/(2*k13)-(1-delta)*k13(-1)+k14+phi*(k14-k14(-1))^zeta/(2*k14)-(1-delta)*k14(-1)+k15+phi*(k15-k15(-1))^zeta/(2*k15)-(1-delta)*k15(-1)+k16+phi*(k16-k16(-1))^zeta/(2*k16)-(1-delta)*k16(-1)+k17+phi*(k17-k17(-1))^zeta/(2*k17)-(1-delta)*k17(-1)+k18+phi*(k18-k18(-1))^zeta/(2*k18)-(1-delta)*k18(-1)+k19+phi*(k19-k19(-1))^zeta/(2*k19)-(1-delta)*k19(-1)+k20+phi*(k20-k20(-1))^zeta/(2*k20)-(1-delta)*k20(-1)+k21+phi*(k21-k21(-1))^zeta/(2*k21)-(1-delta)*k21(-1)+k22+phi*(k22-k22(-1))^zeta/(2*k22)-(1-delta)*k22(-1)+k23+phi*(k23-k23(-1))^zeta/(2*k23)-(1-delta)*k23(-1)+k24+phi*(k24-k24(-1))^zeta/(2*k24)-(1-delta)*k24(-1)+k25+phi*(k25-k25(-1))^zeta/(2*k25)-(1-delta)*k25(-1)+k26+phi*(k26-k26(-1))^zeta/(2*k26)-(1-delta)*k26(-1)+k27+phi*(k27-k27(-1))^zeta/(2*k27)-(1-delta)*k27(-1)+k28+phi*(k28-k28(-1))^zeta/(2*k28)-(1-delta)*k28(-1)+k29+phi*(k29-k29(-1))^zeta/(2*k29)-(1-delta)*k29(-1)+k30+phi*(k30-k30(-1))^zeta/(2*k30)-(1-delta)*k30(-1)+k31+phi*(k31-k31(-1))^zeta/(2*k31)-(1-delta)*k31(-1)+k32+phi*(k32-k32(-1))^zeta/(2*k32)-(1-delta)*k32(-1)+k33+phi*(k33-k33(-1))^zeta/(2*k33)-(1-delta)*k33(-1)+k34+phi*(k34-k34(-1))^zeta/(2*k34)-(1-delta)*k34(-1)+k35+phi*(k35-k35(-1))^zeta/(2*k35)-(1-delta)*k35(-1)+k36+phi*(k36-k36(-1))^zeta/(2*k36)-(1-delta)*k36(-1)+k37+phi*(k37-k37(-1))^zeta/(2*k37)-(1-delta)*k37(-1)+k38+phi*(k38-k38(-1))^zeta/(2*k38)-(1-delta)*k38(-1)+k39+phi*(k39-k39(-1))^zeta/(2*k39)-(1-delta)*k39(-1)+k40+phi*(k40-k40(-1))^zeta/(2*k40)-(1-delta)*k40(-1)+k41+phi*(k41-k41(-1))^zeta/(2*k41)-(1-delta)*k41(-1)+k42+phi*(k42-k42(-1))^zeta/(2*k42)-(1-delta)*k42(-1)+k43+phi*(k43-k43(-1))^zeta/(2*k43)-(1-delta)*k43(-1)+k44+phi*(k44-k44(-1))^zeta/(2*k44)-(1-delta)*k44(-1)+k45+phi*(k45-k45(-1))^zeta/(2*k45)-(1-delta)*k45(-1)+k46+phi*(k46-k46(-1))^zeta/(2*k46)-(1-delta)*k46(-1)+k47+phi*(k47-k47(-1))^zeta/(2*k47)-(1-delta)*k47(-1)+k48+phi*(k48-k48(-1))^zeta/(2*k48)-(1-delta)*k48(-1)+k49+phi*(k49-k49(-1))^zeta/(2*k49)-(1-delta)*k49(-1)+k50+phi*(k50-k50(-1))^zeta/(2*k50)-(1-delta)*k50(-1)= a1*k1(-1)^alpha+a2*k2(-1)^alpha+a3*k3(-1)^alpha+a4*k4(-1)^alpha+a5*k5(-1)^alpha+a6*k6(-1)^alpha+a7*k7(-1)^alpha+a8*k8(-1)^alpha+a9*k9(-1)^alpha+a10*k10(-1)^alpha+a11*k11(-1)^alpha+a12*k12(-1)^alpha+a13*k13(-1)^alpha+a14*k14(-1)^alpha+a15*k15(-1)^alpha+a16*k16(-1)^alpha+a17*k17(-1)^alpha+a18*k18(-1)^alpha+a19*k19(-1)^alpha+a20*k20(-1)^alpha+a21*k21(-1)^alpha+a22*k22(-1)^alpha+a23*k23(-1)^alpha+a24*k24(-1)^alpha+a25*k25(-1)^alpha+a26*k26(-1)^alpha+a27*k27(-1)^alpha+a28*k28(-1)^alpha+a29*k29(-1)^alpha+a30*k30(-1)^alpha+a31*k31(-1)^alpha+a32*k32(-1)^alpha+a33*k33(-1)^alpha+a34*k34(-1)^alpha+a35*k35(-1)^alpha+a36*k36(-1)^alpha+a37*k37(-1)^alpha+a38*k38(-1)^alpha+a39*k39(-1)^alpha+a40*k40(-1)^alpha+a41*k41(-1)^alpha+a42*k42(-1)^alpha+a43*k43(-1)^alpha+a44*k44(-1)^alpha+a45*k45(-1)^alpha+a46*k46(-1)^alpha+a47*k47(-1)^alpha+a48*k48(-1)^alpha+a49*k49(-1)^alpha+a50*k50(-1)^alpha;
end;
initval;
c=1;
e=0;
a1=1;
e1=0;
k1=10;
a2=1;
e2=0;
k2=10;
a3=1;
e3=0;
k3=10;
a4=1;
e4=0;
k4=10;
a5=1;
e5=0;
k5=10;
a6=1;
e6=0;
k6=10;
a7=1;
e7=0;
k7=10;
a8=1;
e8=0;
k8=10;
a9=1;
e9=0;
k9=10;
a10=1;
e10=0;
k10=10;
a11=1;
e11=0;
k11=10;
a12=1;
e12=0;
k12=10;
a13=1;
e13=0;
k13=10;
a14=1;
e14=0;
k14=10;
a15=1;
e15=0;
k15=10;
a16=1;
e16=0;
k16=10;
a17=1;
e17=0;
k17=10;
a18=1;
e18=0;
k18=10;
a19=1;
e19=0;
k19=10;
a20=1;
e20=0;
k20=10;
a21=1;
e21=0;
k21=10;
a22=1;
e22=0;
k22=10;
a23=1;
e23=0;
k23=10;
a24=1;
e24=0;
k24=10;
a25=1;
e25=0;
k25=10;
a26=1;
e26=0;
k26=10;
a27=1;
e27=0;
k27=10;
a28=1;
e28=0;
k28=10;
a29=1;
e29=0;
k29=10;
a30=1;
e30=0;
k30=10;
a31=1;
e31=0;
k31=10;
a32=1;
e32=0;
k32=10;
a33=1;
e33=0;
k33=10;
a34=1;
e34=0;
k34=10;
a35=1;
e35=0;
k35=10;
a36=1;
e36=0;
k36=10;
a37=1;
e37=0;
k37=10;
a38=1;
e38=0;
k38=10;
a39=1;
e39=0;
k39=10;
a40=1;
e40=0;
k40=10;
a41=1;
e41=0;
k41=10;
a42=1;
e42=0;
k42=10;
a43=1;
e43=0;
k43=10;
a44=1;
e44=0;
k44=10;
a45=1;
e45=0;
k45=10;
a46=1;
e46=0;
k46=10;
a47=1;
e47=0;
k47=10;
a48=1;
e48=0;
k48=10;
a49=1;
e49=0;
k49=10;
a50=1;
e50=0;
k50=10;
end;
vcov = [
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1];
order = 2;
\ No newline at end of file
dynare++/tests/sedmodel1.mod deleted 100644 → 0
View file @ 219d2bb3
var A, Disp, G, Int, L,
LStar, pi, Welf, WelfStar, x0,
Y, YGap, YStar, z1, z2, Cbar, Cequiv;
varexo eps1 eps2 eps3;
parameters alpha beta gamma rhoa rhog rho phi chi chi0 theta xi
ABar GBar KBar ZBar piBar Istar;
alpha = 0.3;
beta = 0.99;
gamma = 15;
rhoa = 0.8;
rhog = 0.7;
phi = 1.5;
chi = 1.5;
chi0 = 1;
theta = 0.333333333333;
xi = 0.75;
ABar = 4.0266;
GBar = 0.3163;
KBar = 9.489;
ZBar = .03;
piBar = 1;
rho=.8;
Istar=1.01010101010101;
model;
z1 - ((Y-G)^(1-phi) + beta *xi *piBar *pi(+1)^(1/theta) *z1(+1));
z2 - (Y *chi0 *(1-L-ZBar)^(-chi) / ((1-alpha) *A *KBar^alpha
*L^(-alpha)) + beta *xi *pi(+1)^((1+theta)/theta) *z2(+1));
x0 - (1+theta)*z2 /z1;
pi^(-1/theta) - ((1-xi) *(x0*pi)^(-1/theta) + xi *piBar^(-1/theta));
Y - (Disp^(-1) *A *KBar^alpha *L^(1-alpha));
Disp - ((1-xi) *x0^(-(1+theta)/theta)
+ xi *(piBar/pi)^(-(1+theta)/theta) *Disp(-1));
log(A/ABar) - (rhoa *log(A(-1)/ABar) + eps1);
log(G/GBar) - (rhog *log(G(-1)/GBar) + eps2);
(Y-G)^(-phi) - (beta *(Int/pi(+1)) *(Y(+1)-G(+1))^(-phi));
Welf - ((Y-G)^(1-phi) /(1-phi)
+ chi0*(1-L-ZBar)^(1-chi) /(1-chi) + beta *Welf(+1));
Cequiv = (((1-beta)*Welf-chi0*(1-LStar-ZBar)^(1-chi) /(1-chi))*(1-phi))^(1/(1-phi));
(1-alpha) *A *KBar^alpha *LStar^(-alpha)
- (1+theta) *YStar *(YStar-G)^(phi-1) *chi0
*(1-LStar-ZBar)^(-chi);
YStar - A *KBar^alpha *LStar^(1-alpha);
YGap - (log(Y/YStar));
WelfStar - ((YStar-G)^(1-phi) /(1-phi)
+ chi0*(1-LStar-ZBar)^(1-chi) /(1-chi) + beta *WelfStar(+1));
Int = (Int(-1)^rho)*((Istar*(pi/piBar)^gamma)^(1-rho))*exp(eps3);
Cbar=(1/100)*((1-phi)*((1-beta)*WelfStar-chi0*(1-LStar-ZBar)^(1-chi)/(1-chi)))^(1/(1-phi));
end;
initval;
A= 4.022;
Disp= 1;
G= 0.3;
Int= 1.0101;
L= 0.22;
LStar= 0.22;
pi= 1;
Welf= -359;
WelfStar= -359;
x0= 1;
Y= 2.8;
YGap= 0;
YStar= 2.8;
z1= 2.5;
z2= 1.8;
Cbar= 0.024;
Cequiv = 0.024;
end;
vcov = [0.001 0 0 ; 0 0.001 0; 0 0 0.001];
order=4;