function F = SS_k(xss)
global beta a1 phi delta tau T eta omega alpha fi nu psi1 A1 A2

Y_1 = xss(1);
Y_2 = xss(2);
K_1 = xss(3);
K_2 = xss(4);
L_1 = xss(5);
L_2 = xss(6);
N = xss(7);
P_1 = xss(8);
P_2 = xss(9);
W = xss(10);
R_fd = xss(11);

% Equacoes preliminares
    % Eq. 3.6
    R_n = 1/beta;
    % Eq. 3.7
    Y = ((1-a1)^(1/phi)*Y_1^((phi-1)/phi) + a1^(1/phi)*Y_2^((phi-1)/phi))^(phi/(phi-1));
    % Eq. 3.9 
    P = ((1-a1)*P_1^(1-phi) + a1*P_2^(1-phi))^(1/(1-phi));
    % Eq. 3.31
    Q = - P;
    % Eq. 3.30
    R = Q*(beta*(1-delta) -1)*(1/beta);
    % Eq. 3.32
    R_f = R_n;
    % Eq. 3.37
    L = L_1 + L_2;
    % Eq. 3.38
    K = K_1 + K_2;
    % Eq. 3.26
    I = delta*K;
    % Eq. 3.36
    C = Y - I;
    
% Sistema de equacoes otimizadas
F = [% Restricao orcamentaria 
     W*L*(1-tau) + N*(R_n - 1) - C*P + T;
     % Oferta de trabalho
    (C^eta)*(L^omega)*P - W*(1-tau);
    % Fc de producao 1
    Y_1 - A1*(K_1^alpha)*(L_1^(1-alpha));
    % Fc de producao 2
    Y_2 - A2*(K_2^alpha)*(L_2^(1-alpha));
    % Rendimento de K_1
    ((1-nu)*R_f*R + nu*R_fd*R)/P_1 - ((fi - 1)/fi)*(alpha)*(Y_1/K_1);
    % Rendimento de K_2
    (R*R_f)/P_2 - alpha*((fi - 1)/fi)*(Y_2/K_2);
    % Rendimento de L_1
    W/P_1 - ((fi - 1)/fi)*(1 - alpha)*(Y_1/L_1);
    % Rendimento de L_2
    W/P_2 - ((fi - 1)/fi)*(1 - alpha)*(Y_2/L_2);
    % Preco igual a markup e custo marginal 1
    1 - ((fi-1)/fi)*(1/alpha)^alpha*(1/(1-alpha))^(1-alpha)*((1-nu)*R_f*R+nu*R_fd*R)^(alpha)*W^(1-alpha);
        %((fi-1)/fi)*(((1-nu)*R_f*W + nu*R_fd*W)/(1-alpha))*(((1-alpha)/alpha)*(R/((1-nu)*R_f*W + nu*R_fd*W)))^alpha - 1;
    % Preco igual a markup e custo marginal 2
        %((fi-1)/fi)*((R_f*W)/(1-alpha))*(((1-alpha)/alpha)*(R/(R_f*W)))^alpha - 1;
    1 - ((fi-1)/fi)*(1/alpha)^alpha*(1/(1-alpha))^(1-alpha)*(R_f*R)^(alpha)*(W)^(1-alpha);
    % BNDES
    nu*(R_f - R_fd)*R*K_1 - tau*W*L;
    % Capital 1
    R/P - psi1;
    % SF
    R*(L_1*(1-nu) + L_2) - N;
    % Composicao
    Y_1/Y_2 - ((1-a1)/a1)*(P_2/P_1)^phi;    
    ];
