function ff=findssFF_prud(guess)

global d0 d1 d2  betta gama kappa alfab alfag deltab deltag taubar theta1 theta2 ...
       rhoL deltax zetta taug taub R omega erow ksi rhoy alfayb alfayg Abar

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));

Kb = Kb_Lb*Lb;
Kg = Kg_Lg*Lg;

Qg = 1;
Qb = 1;

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

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

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);

mu = 1-emis/Yb;

Rkb = ((1-dX)*(Kb_Lb)^(alfab-1))*(alfab*(pb-theta1*mu^theta2-taubar*(1-mu)))/Qb + (1-deltab);

Rkg = Rkb- (taub-taug)*R;

p  =  Rkb-(1+taub)*R;


ff1 = Rkg -(pg*(1-dX)*alfag*(Kg_Lg)^(alfag-1) + Qg*(1-deltag))/Qg;

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

ff3 = omega*(L^(ksi-rhoL))*Lb^rhoL - ...
    (1-alfab)*(1-dX)*Abar*(Kb_Lb)^alfab*(pb-theta1*mu^theta2-taubar*(1-mu));

ff4 = taubar-theta1*theta2*mu^(theta2-1);

ff5 = (1-gama)*betta*(betta*p+zetta)*(gama*p+zetta) - kappa*(betta-gama)*((1-betta)*gama*p+(1-gama)*zetta);

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


