function f=onefoc(z,param);
%ONEFOC  Steady-state first order conditions for model described in onesect.f.

%        Ellen McGrattan, 6-23-97
%        Revised, ERM, 9-24-97

a     = param(1);
alpk  = param(2);
alph  = param(3);
beta  = param(4);
delk  = param(5);
delh  = param(6);
eta   = param(7);
psi   = param(8);
sigma = param(9);
zeta  = param(10);
tauxk = param(11);
tauxh = param(12);
tauk  = param(13);
tauh  = param(14);
tauc  = param(15);
z1    = param(16);
z2    = param(17);
alp   = alpk+alph;
sigz  = sigma-zeta*(1-sigma);
betn  = beta*(1+eta/100)^sigz;
tauxk1= 1+tauxk;
tauxh1= 1+tauxh;
tauk1 = 1-tauk;
tauh1 = 1-tauh;
tauc1 = 1+tauc;
delk1 = 1-delk;
delh1 = 1-delh;

k     = z(1);
l     = z(2);
gam   = z(3);
xh    = gam-delh1;
xk    = (gam-delk1)*k;

gs    = gam^sigz;
g2    = (z2*a*k^alpk*l^alph)^(1/alp);
F     = a*k^alpk*l^alph*g2^(1-alpk-alph);
g1    = z1*F;
c     = F-xk-xh-g1-g2;
F1    = alpk * F/k;
F2    = alph * F/l;

f(1)  = psi*c - (1-l)* F2 * tauh1/tauc1;
f(2)  = gs-betn*(delk1+F1*tauk1/tauxk1);
f(3)  = gs-betn*(delh1+F2*l*tauh1/tauxh1);
f     = f(:);
