function F = solvess5b(x,mc1_ss,gamma,aa,chi,raf1_ss,yf_ss,af1_ss,eta,Gamma1,delta1,alpha,rk1_ss,q1_ss,rl1_ss,l_y)
%Solve for the steady state variables for full model with money


c_temp=(1-alpha)*mc1_ss/(chi*x(1)^eta*Gamma1^alpha);
i_temp=delta1*Gamma1^-1*x(1);
ce_temp=(rk1_ss+(1-delta1)*q1_ss)*Gamma1^-1*x(1) - q1_ss*Gamma1^-1*x(1)-rl1_ss*x(4)*x(1)*Gamma1^-alpha*l_y;
         %(rk1+(1-delta1)*q1)*k1(-1)-q1*k1 = ce1

F=[ x(4)*x(1)*Gamma1^-alpha - ( (1-gamma)*x(4)^-aa*(c_temp +i_temp+ce_temp) + x(2)^(-aa)*yf_ss);
    gamma*x(3)^(1-aa)*  (c_temp+i_temp+ce_temp) - (x(3)*x(2)*x(2)^(-aa)*yf_ss+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);
%         ph1_ss=solxx1(4);