addpath /Applications/Dynare/4.3.2/matlab

var n c hf hm ef em k lambda r wf wm y s;
parameters rho phi gammaf gammam nu theta alpha beta kappa eta cs efs ems;

rho = 0.97;
theta = 0.4;
kappa = 0.33;
gammaf = 0.6;
gammam = 0.8;
alpha = 0.3;
beta = 0.3;
nu = 0.15;
eta = 0.3;


cs=61.3340;
efs=239.5615;
ems=225.4019;

model;
(1/c) = lambda;
(rho/c(+1)) = (1/(1+r))*lambda;
theta*phi*(n/2)^gammaf*(1/ef) = lambda*(n/2);
theta*phi*(n/2)^gammam*(1/em) = lambda*(n/2);
(gammaf*phi/2)*(n/2)^(gammaf-1)*log(ef^theta*hf(-1)^alpha*hm(-1)^beta)+(gammam*phi/2)*(n/2)^(gammam-1)*log(em^theta*hf(-1)^alpha*hm(-1)^beta) = lambda*(nu*wf*hf(-1)+(1/2)*ef+(1/2)*em);
c+(n/2)*ef+(n/2)*em+s = wf*hf(-1)*(1-nu*n)+wm*hm(-1);
hf = ef^theta*hf(-1)^alpha*hm(-1)^beta;
hm = em^theta*hf(-1)^alpha*hm(-1)^beta;
r = kappa*k(-1)^(kappa-1);
wf = eta*(1-kappa)*k(-1)^kappa;
wm = (1-eta)*(1-kappa)*k(-1)^kappa;
y = k(-1)^kappa;
k = s/(1+(n/2));
end;

initval;


hf = 235.2227;
hm = 229.5596;
k = 34.2432;



n = 1.4748;
c = 61.3340;
s= 59.4956;
ef = 239.5615;
em = 225.4019;
lambda = 0.0163;
r = 0.0309;
wf = 0.6451;
wm = 1.5052;
y = 3.2093; 

end;



steady;