function F = SS_DK(xss)
global beta eta sigma lambda theta gamma psi tau alpha delta_a delta_b zeta epsilon  omega 

K = xss(1);
N = xss(2);
B = xss(3);
PHI = xss(4);
ETA = xss(5);
NU = xss(6);
X = xss(7);
R_f = xss(8);
Z = xss(9);
Ne = xss(10);
Nn = xss(11);
G = xss(12);
lucro = xss(13);
Ym = xss(14);
Y = xss(15);
L = xss(16);
C = xss(17);
DD = xss(18);
NN = xss(19);

S_C = 1;
S_L = 1;
A = 1;
S_K = 1;
R = 1/beta;
In = 0.001;
U = 1;
DELTA = delta_a + delta_b*U^(1+zeta)/(1+zeta);
PSI = psi;
LAMBDA = 1;
Pi = 1;
Pistar = Pi;
D = 1;
Pm = (epsilon - 1)/epsilon;
in = 1/beta - 1;
Q = 1;
W = L^sigma*C^eta;


S = K;
Sg = PSI*S;
Bg = Sg;
Sp = S - Sg;
I = DELTA*K;

F = [Sp - N - B;
    Sp - PHI*N;
    PHI*(lambda - NU) - ETA;
    NU*(1 - beta*theta*X) - (1-theta)*beta*(R_f - R);
    X - Z;
    ETA*(1 - theta*beta*Z) - (1-theta)*beta*R;
    Z - (R_f - R)*PHI - R;
    N - Ne - Nn;
    Ne - theta*Z*N;
    Nn - omega*Sp;
    G + tau*PSI*K - (R_f - R)*Bg;
    lucro - alpha*Pm*Ym + DELTA*K
    K*(R_f - 1) - lucro;
    W*L - (1 - alpha)*Pm*Ym;
    delta_b*U^(1+zeta)*K - alpha*Pm*Ym;
    Ym - K^alpha*L^(1-alpha);
    Y - C - I - (R_f - R)*Bg;
    DD*(epsilon - 1) - epsilon*NN;
    NN*(1-gamma*beta) - Pm*Y;
    DD*(1-gamma*beta) - Ym;
    Y - Ym
    ];
