infl=Pis;
R=infl/betal;
Lam=betal;
Rr=Rrs;
R=Rs;
A=As;
s=ss;
pmF=s;
Rf=sp+R;
Qf=1/(Rf-kappa);
Rb=Delta*sp+Rs;
Qb=1/(Rb-kappa);
Rhm=R*exp(-zeta*s);
vpr=vprs;
pF=pFs;
pFr=pFrs;
pH=pHs;
X=Xs;
EX=X*(pH/s)^(-v);
hm=hms;
YmF=(EX-Rhm*s*hm)/pmF;
YF=vpr*YmF;
x1r=x1rs;
x2r=x2rs;
YhH=gamma/(1-gamma)*YF;
YH=YhH+EX;
Y=Ys;
pr=prs;
vp=vps;
x1=(pr/pH)^(1-epsilon)*YH / (1-Phip*inflH^(1-epsilon)*Lam*(pr/pH)^(1-epsilon));
x2=x1*(epsilon-1)/epsilon;
pm=x2*pH/YH*((pr/pH)^epsilon-Phip*inflH^(-epsilon)*Lam);
Ym=vp*YH;
theta1=Lam/(Qf*infl*(1-kappa*Lam/infl))-1;
theta2=theta1+1;
K=alpha*Lam*pm*Ym/(theta2*(1-Lam*(1-delta)));
L=(Ym/(A*K^alpha))^(1/(1-alpha));
w=(1-alpha)*pm*Ym/L;
I=delta*K;
f=I/(Qf*(1-kappa/infl));
laml=lamls;
lam1=lam1s;
lam2=lam1*(Qb-betab*(1+kappa*Qb)/infl) / (Qb*(1-kappa*betab/infl));
Cl=1/laml;
Cb=1/lam1;
qh=gammah/(laml*(1-betal))+gammah/(lam1*(1-betab)-xi*lam2);
Hl=gammah/(qh*laml*(1-betal));
Hb=gammah/(qh*(lam1*(1-betab)-xi*lam2));
b=xi*qh*Hb/(Qb*(1-kappa/infl));
Lb=Lbs;
Ll=Lls;
n=ns;
d=Qf*f+Qb*b-n-s*hm;
phi=phis;
Omega=1-sigma+sigma*theta*phi;
lam=lams;


ys = [ log(Cl);
    log(Hl);
    log(Ll);
    log(laml);
    log(Cb);
    log(Hb);
    log(b);
    log(Lb);
    log(lam1);
    log(lam2);
    log(Rf);
    log(Rb);
    log(Rhm);
    log(Lam);
    log(Omega);
    log(n);
    log(f);
    log(d);
    log(lam);
    log(phi);
    log(Qb);
    log(Y);
    log(YhH);
    log(pH);
    log(Ym);
    log(theta1);
    log(theta2);
    log(L);
    log(K);
    log(I);
    log(pr);
    log(pm);
    log(YH);
    log(x1);
    log(x2);
    log(vp);
    log(pFr);
    log(pF);
    log(YF);
    log(x1r);
    log(x2r);
    log(vpr);
    log(R);
    log(pmF);
    log(EX);
    log(X);
    log(YmF);
    log(hm);
    log(Rr);
    log(infl);
    log(w);
    log(qh);
    log(Qf);
    log(s);
    log(A)];
    