infl=Pis;
Lam=beta;
Rd=Pis/Lam;
RF=Rd+sp;
RB=Delta*(RF-Rd)+Rd;
Q=1/(RF-kappa);
QB=1/(RB-kappa);
pr=(    (1-phip*Pis^((gammap-1)*(1-epsilon)))       /    (1-phip)       )^(1/(1-epsilon));
pm=pr*(epsilon-1)/epsilon;
L=Ls;
u=1-L;
A=1;
K=L*(   (1-Lam*(1-delta))      /     (alpha*pm*A*Q*Pis*(1-kappa*Lam/Pis))     )^(1/(alpha-1));
Ym=A*K^alpha*L^(1-alpha);
I=delta*K;
vp=(1-phip)*pr^(-epsilon)/(1-phip*Pis^((1-gammap)*epsilon));
Y=Ym/vp;
x1=pm*Y/(1-phip*Lam);
x2=Y/(1-phip*Lam);
fm=pk*delta*K/(Q*(1-kappa/Pis));
v1=Lam/(Q*(Pis-Lam*kappa));
v2=v1-1;
Ih=I;
G=Gs;
C=Y-I-G;
mu=(1-bh*beta)/(C*(1-bh));
Qgamma=Qgammas;
Gamma=deltaN*L/Qgamma;
U=1-(1-deltaN)*L;
M=mubar*U^mum*Gamma^(1-mum);
QU=mubar*U^(mum-1)*Gamma^(1-mum);
Theta=xi/Qgamma;
w=(1-alpha)*pm*Y/L-Theta+Lam*Theta*(1-deltaN);
RTR=Rd;
Rre=RTR;
fcb=fcbs;
bcb=bcbs;
re=Q*fcb+QB*bcb;
f=fm-fcb;
b=bG-bcb;
n=((RF-Rd)*Q*f+(RB-Rd)*QB*b+(Rre-Rd)*re)     /       (Pis-Rd);
phi=(Q*f+Delta*QB*b)/n;
theta=thetas;
Omega=1-sigma+sigma*theta+sigma*theta*phi;
lam=(Lam*Omega*(RF-Rd))     /     (theta*Pis-Lam*Omega*(RF-Rd));
omega=-(1+lam)*Lam*Omega*(Rre-Rd)/Pis;
T=G+bG/Pis-QB*bG*(1-kappa/Pis);
d=Q*f+QB*b+re-n;
lev=re/d;
div=n-n*Rd/Pis+Y-w*L-fm/Pis+fm*Q*(1-kappa/Pis)-xi*Gamma-I;
pk=1;


ys = [log(infl);
    log(Lam);
    log(Rd);
    log(RF);
    log(RB);
    log(Q);
    log(QB);
    log(pr);
    log(pm);
    log(L);
    log(u);
    log(K);
    log(Ym);
    log(I);
    log(pk);
    vp;
    log(Y);
    x1;
    x2;
    log(fm);
    v1;
    v2;
    log(Ih);
    log(G);
    log(C);
    log(mu);
    log(Qgammas);
    log(Gamma);
    log(U);
    log(M);
    log(QU);
    log(Theta);
    log(w);
    log(RTR);
    log(Rre);
    fcb;
    bcb;
    re;
    log(f);
    log(b);
    log(n);
    log(phi);
    log(thetas);
    Omega;
    lam;
    omega;
    T;
    log(d);
    lev;
    log(A);
    div];

  