clear;
global beta eta sigma lambda theta gamma psi tau alpha delta_a delta_b zeta epsilon varphi nu eta_1 omega gamma_p kappa_pi kappa_y
    beta  = .989;
    sigma = 1.5;
    eta   = 2;
    alpha = .39;
    epsilon= 4.16;
    delta_a = .0204;
    delta_b = .0376;
    zeta = 7.2;
    tau = .2;
    psi = .3;
    nu = 10;
    eta_1 = 3.42;
    lambda = .2800;
    theta = .9750;
    omega = .00009;
    gamma = .74;
    gamma_p = .33;
    kappa_pi = 2.43;
    kappa_y = .5;
 
xss0 = [19 2.3 11 2 1.6 .002 1.02 1.016 1.02 2.3 .001 -1.13 .35 2.40 2.40 .6 1.9 7 6];
options = optimset('Display','off');
[xss,Fval,exitflag] = fsolve(@SS_DK,xss0,options);
exitflag

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;

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

W = L^omega*C^eta;

