function F = solvess2b(x,mc1_ss,gamma,aa,chi,raf1_ss,yf_ss,af1_ss,eta)
%Solve for the steady state variables for full model with money


        c1_temp   = mc1_ss/(chi*x(1)^eta);

F=[ x(4)*x(1) - (1-gamma)*x(4)^(-aa)*c1_temp-x(2)^(-aa)*yf_ss;
    

    x(3)*af1_ss+gamma*x(3)^(1-aa)*c1_temp   - x(3)*x(2)*x(2)^(-aa)*yf_ss-x(3)*(1+raf1_ss)*af1_ss;
    
    
    1-x(2)*x(3)/x(4)
    1 - (((1-gamma)*x(4)^(1-aa)+ gamma*x(3)^(1-aa)))] ;




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