function F = steady_state_solver(x,para)

%variables conventions
% x(1) ---> R
% x(2) ---> Rk
% x(3) ---> I
% x(4) ---> K
% x(5) ---> Wreal
% x(6) ---> C
% x(7) ---> N
% x(8) ---> Yh
% x(9) ---> Yh*=Yf
% x(10)---> Y
% x(11)---> Ph/P
% x(12)---> Ph*/P*=Pf/P
% x(13)---> MC
% x(14)---> G
% x(15)---> RER

%setting parameters
beta=para(1,1);
delta=para(1,2);
b=para(1,3);
ew=para(1,4);
eta=para(1,5);
alpha=para(1,6);
ep=para(1,7);
omega=para(1,8);
theta=para(1,9);
percG=para(1,10);
uf=para(1,11); % Assuming ef=ep=11
m=1;

F = [x(1)-(1/beta); %HH1
     x(2)*m-(1/beta)+(1-delta); %HH2
     x(3)-delta*x(4); %HH4
     (x(5)/(x(6)*(1-b)))-(ew/(ew-1))*x(7)^eta; %Wage-setting
     x(7)/(x(4)*m)-((1-alpha)/alpha)*(x(2)/x(5)); %Capital-labour ratio
     x(8)+x(9)-(x(4)*m)^alpha*x(7)^(1-alpha);
     x(11)-(ep/(ep-1))*x(13);
     x(15)*x(12)-(ep/(ep-1))*x(13);
     x(8)/x(9)-(omega/(1-omega))*(x(11)/x(12))^(-theta);
     x(10)-(omega^(1/theta)*x(8)^((theta-1)/theta)+(1-omega)^(1/theta)*x(9)^((theta-1)/theta))^(theta/(theta-1));
     x(10)-x(6)-x(3)-x(14); % Market clearing
     x(14)-percG*x(10); %Gvt expenditure
     1-uf*(omega*x(11)^(1-theta)+(1-omega)*x(12)^(1-theta))^(1/(1-theta));
     x(13)-(x(2)^(alpha)*x(5)^(1-alpha))/((1-alpha)^(1-alpha)*alpha^alpha); %Marginal cost
     x(15)*x(12)*x(9)-x(12)*x(9)];