% File name: financial_friction_num_eval.m 
% File generated by anal_deriv_print2f.m Date: 11-Sep-2009

nfx=[[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-yy/yyback*gback]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-c/cback*gback]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-ivv/ivvback*gback]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[yy/yyback*gback]
[c/cback*gback]
[ivv/ivvback*gback]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-PHI*g*(g-G)]
[a*k^ALFA*ALFA*(g*h)^(1-ALFA)]
[-(1-DELTA)*k]
[0]
[(1-ALFA)*a*g^(1-ALFA)*(k/h)^ALFA*ALFA]
[nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*PHI*g]
[0]
[0]
[0]
[PHI*(g-G)*g/yy]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[d+d^2*g/(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)^2*PSSI*exp(d-DBAR)]
[0]
[0]
[BETTA/(g^GAMA)*PSSI*d*exp(d-DBAR)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)]
[0]
[-nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*(1+PHI*(g-G))]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[RHONU]
[0]
[0]
[d*g/(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)^2*mu*exp(mu-1)]
[0]
[0]
[BETTA/(g^GAMA)*mu*exp(mu-1)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[RHOMU]
[0]
[s]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-s/yy]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[RHOS]
[PHI*g*(g-G)-d*g/(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)]
[a*k^ALFA*(g*h)^(1-ALFA)*(1-ALFA)]
[k*g]
[-BETTA/(g^GAMA)*(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*GAMA]
[(1-ALFA)^2*a*g^(1-ALFA)*(k/h)^ALFA]
[-nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*PHI*g-BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*(1-DELTA+ALFA*a*(g*h/k)^(1-ALFA)+PHI*g*(g-G)-1/2*PHI*(g-G)^2)*GAMA]
[0]
[0]
[RHOG]
[-PHI*(g-G)*g/yy]
[0]
[0]
[0]
[0]
[0]
[0]
[g]
[0]
[0]
[0]
[0]
[a*k^ALFA*(g*h)^(1-ALFA)]
[0]
[0]
[(1-ALFA)*a*g^(1-ALFA)*(k/h)^ALFA]
[0]
[0]
[RHOA]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]];
nfxp=[[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-yyback]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-cback]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-ivvback]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-gback]
[0]
[0]
[0]
[PHI*g*(g-G)]
[0]
[k*g]
[0]
[0]
[-nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*PHI*g+BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*ALFA*a*(g*h/k)^(1-ALFA)*(-1+ALFA)]
[k]
[0]
[0]
[-PHI*(g-G)*g/yy]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-d*g/(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[BETTA/(g^GAMA)*(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)]
[0]
[BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*(1-DELTA+ALFA*a*(g*h/k)^(1-ALFA)+PHI*g*(g-G)-1/2*PHI*(g-G)^2)]
[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]
[-1]
[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]
[BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*(ALFA*a*(g*h/k)^(1-ALFA)*(1-ALFA)+PHI*g^2)]
[0]
[0]
[-1]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*ALFA*a*(g*h/k)^(1-ALFA)]
[0]
[-1]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]];
nfy=[[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-gy]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-gc]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-giv]
[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]
[a*k^ALFA*(g*h)^(1-ALFA)*(1-ALFA)]
[0]
[-nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*GAMA*THETA*h^OMEGA/(c-THETA/OMEGA*h^OMEGA)]
[-THETA*h^(OMEGA-1)*(OMEGA-1)-(1-ALFA)*a*g^(1-ALFA)*(k/h)^ALFA*ALFA]
[-nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*GAMA*THETA*h^OMEGA/(c-THETA/OMEGA*h^OMEGA)*(1+PHI*(g-G))]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-yy]
[-yy]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[1-(yy-c-s-ivv-1/2*PHI*(g-G)^2)/yy]
[yy/yyback*gback]
[0]
[0]
[yy]
[0]
[0]
[0]
[0]
[0]
[0]
[c]
[0]
[0]
[nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*GAMA*c/(c-THETA/OMEGA*h^OMEGA)]
[0]
[nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*GAMA*c/(c-THETA/OMEGA*h^OMEGA)*(1+PHI*(g-G))]
[0]
[0]
[0]
[-c/yy]
[0]
[c/cback*gback]
[0]
[0]
[c]
[0]
[0]
[0]
[0]
[0]
[ivv]
[0]
[-ivv]
[0]
[0]
[0]
[0]
[0]
[0]
[-ivv/yy]
[0]
[0]
[ivv/ivvback*gback]
[0]
[0]
[ivv]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*PHI*g^2]
[-k1]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]];
nfyp=[[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[BETTA/(g^GAMA)*(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*GAMA*THETA*h^OMEGA/(c-THETA/OMEGA*h^OMEGA)]
[0]
[BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*GAMA*THETA*h^OMEGA/(c-THETA/OMEGA*h^OMEGA)*(1-DELTA+ALFA*a*(g*h/k)^(1-ALFA)+PHI*g*(g-G)-1/2*PHI*(g-G)^2)+BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*ALFA*a*(g*h/k)^(1-ALFA)*(1-ALFA)]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[-BETTA/(g^GAMA)*(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*GAMA*c/(c-THETA/OMEGA*h^OMEGA)]
[0]
[-BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*GAMA*c/(c-THETA/OMEGA*h^OMEGA)*(1-DELTA+ALFA*a*(g*h/k)^(1-ALFA)+PHI*g*(g-G)-1/2*PHI*(g-G)^2)]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*PHI*g^2]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]];
nf=[[d-yy+c+s+ivv+1/2*PHI*(g-G)^2-d*g/(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)]
[-yy+a*k^ALFA*(g*h)^(1-ALFA)]
[-ivv+k*g-(1-DELTA)*k]
[-nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)+BETTA/(g^GAMA)*(RSTAR+PSSI*(exp(d-DBAR)-1)+exp(mu-1)-1)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)]
[-THETA*h^(OMEGA-1)+(1-ALFA)*a*g^(1-ALFA)*(k/h)^ALFA]
[-nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*(1+PHI*(g-G))+BETTA/(g^GAMA)*nu*(c-THETA/OMEGA*h^OMEGA)^(-GAMA)*(1-DELTA+ALFA*a*(g*h/k)^(1-ALFA)+PHI*g*(g-G)-1/2*PHI*(g-G)^2)]
[-k1+k]
[-log(a)+RHOA*log(a)]
[-log(g/G)+RHOG*log(g/G)]
[-tby+(yy-c-s-ivv-1/2*PHI*(g-G)^2)/yy]
[-gy+yy/yyback*gback]
[-gc+c/cback*gback]
[-giv+ivv/ivvback*gback]
[-yyback+yy]
[-cback+c]
[-ivvback+ivv]
[-gback+g]
[-log(nu)+RHONU*log(nu)]
[-log(mu)+RHOMU*log(mu)]
[-log(s/S)+RHOS*log(s/S)]];
nETASHOCK=[[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[SIGMAG]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[SIGMAA]
[0]
[0]
[0]
[0]
[0]
[0]
[SIGMANU]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[SIGMAMU]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[0]
[SIGMAS]
[0]
[0]];
nfx=reshape(nfx,[20  11]);
nfxp=reshape(nfxp,[20  11]);
nfy=reshape(nfy,[20   9]);
nfyp=reshape(nfyp,[20   9]);
nf=reshape(nf,[20   1]);

nETASHOCK=reshape(nETASHOCK,[11   5]);