function F = solvess8b (x,Gamma1,gamma,aa,raf1_ss,yf_ss,af1_ss,alpha,h1_ss,i1_ss,ce1_ss,delta1_b,kb1_ss)
%Solve for the steady state variables for full model with money


 c1_temp   = (x(1)*Gamma1^(-alpha)*h1_ss-x(2)^(-aa)*yf_ss-delta1_b*kb1_ss)/(x(1)^(-aa)*(1-gamma))-i1_ss-ce1_ss;
 
        
F=[  (1+raf1_ss)*af1_ss+gamma*x(3)^(1-aa)*(c1_temp+i1_ss+ce1_ss) - x(3)*x(2)*x(2)^(-aa)*yf_ss-af1_ss;  
     1 - (((1-gamma)*x(1)^(1-aa)+ gamma*x(3)^(1-aa)));
     1-x(3)*x(2)/x(1)] ;

%         ph1_ss=solxx1(1);
%         px1_ss=solxx1(2);
%         pm1_ss=solxx1(3)
%         RER1_ss=pm1_ss ;

