function ff=findss_nopol(guess)

global  alfab alfag deltab deltag omega ksi rhoL deltax R Abar erow rhoy alfayb alfayg d0 d1 d2

Kb_Lb=guess(1);
Kg_Lg=guess(2);
Lb=guess(3);
Lg=guess(4);
emis=guess(5);

L = (Lb^(1+rhoL) + Lg^(1+rhoL))^(1/(1+rhoL));

X = (erow+emis)/(1-deltax);
dX = d2*(X)^2+d1*(X)+d0;

Yb = (1-dX)*Abar*Kb_Lb^(alfab)*Lb;
Yg = (1-dX)*Abar*Kg_Lg^(alfag)*Lg;

Y = ((alfayb^(1/rhoy))*Yb^((rhoy-1)/rhoy)+(alfayg^(1/rhoy))*Yg^((rhoy-1)/rhoy))^(rhoy/(rhoy-1));

pg=(alfayg*Y/Yg)^(1/rhoy);
pb=(alfayb*Y/Yb)^(1/rhoy);

Rkg = R;
Rkb = Rkg;

ff1 = omega*(L^(ksi-rhoL))*Lg^rhoL - (1-alfag)*pg*(1-dX)*Abar*(Kg_Lg)^alfag;

ff2 = omega*(L^(ksi-rhoL))*Lb^rhoL - (1-alfab)*pb*(1-dX)*Abar*(Kb_Lb)^alfab;

ff3 = Kg_Lg - (pg*(1-dX)*Abar*alfag/(Rkg-(1-deltag)))^(1/(1-alfag));

ff4 = Kb_Lb - (pb*(1-dX)*Abar*alfab/(Rkb-(1-deltab)))^(1/(1-alfab));

ff5 = emis-Yb;

ff=[ff1; ff2; ff3; ff4; ff5];


