function F=ss_dynare_forum(x)
  
 
%% Variables
c=x(1);
w=x(2);
tau=x(3);
l_a=x(4);
y_a=x(5);
x_a=x(6);
k_a=x(7);
b_a=x(8);
c_a=x(9);
nw_a=x(10);
lamb2_a=x(11);
p_a=x(12);
l_b=x(13);
y_b=x(14);
x_b=x(15);
k_b=x(16);
b_b=x(17);
c_b=x(18);
nw_b=x(19);
lamb2_b=x(20);
p_b=x(21);
mc=x(22);
r1=x(23);
r2=x(24);
i=x(25);
PED=x(26);
R=x(27);
y=x(28);
k=x(29);
l=x(30);
R_b=x(31); 
   
  
% Parameters
phi=1.45;
beta=0.985;          
alpha=0.33;       
beta_E=0.97;         
delta=0.025;      
a=0.5;                 
b=0.5;                 
rho=-10;                                
m_a=0.2;                                                           
m_b=0.2;                      
N=5;                                       
R_s=0.25;                
%% Steady state system

F = [% (01) - Determinacao de R
     R - (1/beta); 
     
     % (02) - Determinacao de w  --> nao consegui isolar w q da pau no sistema --> w - ((phi * c)/((1-l)*(1-tau)));
     phi - (((1-l)*w*(1-tau))/c);
     
     % (03) - Determinacao de l_a
     l_a - (((1-alpha)*y_a)/(x_a*w));
     
     % (04) - Determinacao de x_a
     x_a - (1/p_a);
     
     % (05) - Determinacao de p_a 
     p_a - ((1/(beta_E*alpha*y_a))*((k_a)*( 1 - beta_E*(1-delta) - c_a*lamb2_a*(1-delta)*(m_a/R_s))));     
     
     % (06) - Determinacao de b_a
     b_a - ((k_a) - (beta_E*nw_a));
     
     % (07) - Determinacao de k_a
     k_a - ((R_s*b_a)/(m_a*(1-delta)));
     
     % (08) - Determinacao de y_a
     y_a - ((k_a^alpha)*(l_a^(1-alpha))); 
     
     % (09) - Determinacao de c_a
     c_a - ((1-beta_E)*nw_a);
     
     % (10) - Determinacao de nw_a
     nw_a - ((alpha*y_a*p_a) + ((1-delta)*k_a) - (R_s*b_a));
     
     % (11) - Determinacao de lamb2_a
     lamb2_a - ((1-beta_E*R_s)/c_a);
     
     % (12) - Determinacao de l_b
     l_b - (((1-alpha)*y_b)/(x_b*w));
     
     % (13) - Determinacao de x_b
     x_b - (1/p_b);
     
     % (14) - Determinacao de p_b  
     p_b - ((1/(beta_E*alpha*y_b))*((k_b)*( 1 - beta_E*(1-delta) - c_b*lamb2_b*(1-delta)*(m_b/R_b))));

     % (15) - Determinacao de b_b
     b_b -  (k_b - (beta_E*nw_b));
     
     % (16) - Determinacao de k_b
     k_b - ((R_b*b_b)/(m_b*(1-delta)));     
     
     % (17) - Determinacao de y_b
     y_b - ((k_b^alpha)*(l_b^(1-alpha)));
     
     % (18) - Determinacao de c_b
     c_b - ((1-beta_E)*nw_b);
     
     % (19) - Determinacao de nw_b
     nw_b - ((alpha*y_b*p_b) + ((1-delta)*k_b) - (R_b*b_b));
     
     % (20) - Determinacao de lamb2_b
     lamb2_b - ((1-beta_E*R_b)/c_b);
     
     % (21) - Determinacao de mc
     mc - (p_a/((a+b*r1)^(1/rho)) + p_b/((a*r2+ b)^(1/rho))); 
     
     % (22) - Determinacao de r1
     r1 - (((p_b*a)/(p_a*b))^((rho)/((rho)-1)));
     
     % (23) - Determinacao de r2
     r2 - (((p_a*b)/(p_b*a))^((rho)/((rho)-1)));
     
     % (24) - Determinacao de i
     i - (k*delta);
     
     % (25) - Determinacao de PED
     PED - (1 + (1/(1-alpha))*((m_b * (1 - delta))/(R_b * beta * alpha * p_b * ((k_b^(alpha-1)) * (l_b^(1-alpha))))));
     
     % (26) - Determinacao de R_b
     R_b - ((R)/(1 - (1/N) * (1/PED)));
     
     % (27) - Determinacao de tau
     tau - ((b_a*(R-R_s))/(w*(l_a+l_b)));
     
     % (28) - Determinacao de c
     c -  (-y + c_a + c_b + i);
     
     % (29) - Determinacao de y
     y - ((a*((y_a^(rho))) + b*((y_b^(rho))))^(1/(rho)));
     
     % (30) - Determinacao de k
     k - (k_a + k_b);
     
     % (31) - Determinacao de l
     l - (l_a + l_b)
     ];
end
    
%% Now put this on command window:
  
%Initial guess  -->>  x0_all = [0.5;1.1;1;0.15;0.45;1;2;1.5;0.25;1.65;1;1;0.15;0.45;1;2;1.5;0.25;1.65;1;1;0.9;1;1;0.1;20;1;1;0.7;4;1.2];
%options = optimoptions('fsolve','Display','iter');
%Call fsolve -->> [x,fval] = fsolve(@ss_dynare_forum,x0_all,options)