function F = ss_matheus(x)

w = x(1);
l = x(2);
k = x(3);
cE = x(4);
yw = x(5);
PED = x(6);
Rb = x(7);

x = 5/6;  %% Mark up
phi_l = 1.45;  %% elasticidade do labor
beta = 0.4;        
alpha = 0.33;      
betaE = 0.97;        
delta = 0.025;                  
mk = 0.3;  
N = 6;
%% Steady state system
  
F = [% (02) - Determinacao de w  
     phi_l - (1/( ((k)^(alpha))*(l^(1-alpha)) - cE - (delta*k)))*(1-l)*w;  
      
     % (03) - Determinacao de l
     l -  (1-alpha)*(k^(alpha))  / (w*x)^(1/alpha);
      
     % (05) - Determinacao de k
     1 -  betaE * ( ( alpha * (k)^(alpha)*(l)^(1-alpha) )/k*x )  + betaE *(1-delta)  + ((1 - betaE * 1/beta) * mk * (1-delta)/Rb) ;
	 
     % (09) - Determinacao de cE
     cE - (1 - betaE)*(alpha/x)*( ((k)^(alpha))*(l^(1-alpha)) ) + (1-delta)*k - (1/(beta))*(mk *k *(1-delta)/ Rb);
      
   	 % (17) - Determinacao de yw
     yw - (( (k)^(alpha) ) * ( l^(1-alpha)));
        
     % (24) - Determinacao de PED
     PED - (1 + ( 1 / (1-alpha) ) * ( mk /(1/(beta))) * ( 1 /( betaE * (alpha * ( (k)^(alpha-1) ) * ( l^(1-alpha) ) ) / x) )); 
      
     % Determinação do Rb
     Rb - (((1/(beta)))/(1 - (1/N) * (1/(PED))));
     ];
end

%Initial guess  -->>  x0_7 = [1;1;1;1;1;1;1];
%options = optimoptions('fsolve','Display','iter');
%Call fsolve -->> [x,fval] = fsolve(@ss_matheus,x0_6,options)