function Result = Function19_01_SS(x,PAR,EXO)

alpha     = PAR.alpha;
epsilon_c = PAR.epsilon_c;
epsilon_l = PAR.epsilon_l;
epsilon_m = PAR.epsilon_m;
delta     = PAR.delta;
theta     = PAR.theta;
rho       = PAR.rho;
Omega_0   = PAR.Omega_0;
G         = EXO.G;
M_0       = EXO.M_0;
Lbar      = PAR.Lbar;

% Computes the steady-state equilibrium
% x(1)     is Y
% x(2)     is K
% x(3)     is L
% x(4)     is C
% x(5)     is w
% x(6)     is R
% x(7)     is mc
% x(8)     is P

%     Result(1) = Y0 - C0 - delta * K0 - G ;
%     Result(2) = R0 - rho;
%     Result(3) = M_0 / P0 - (epsilon_m/epsilon_c) * C0 * (1 + 1/R0) ;
%     Result(4) = 1 - L0 - (epsilon_l/epsilon_c) * C0 / w0;
%     Result(5) = w0 - (1 - alpha) * mc0 * Y0 / (L0 - Lbar);   
%     Result(6) = rho + delta - alpha * mc0 * Y0 / K0;
%     Result(7) = mc0 - (theta - 1) / theta;
%     Result(8) = Y0 - Omega_0 * K0^(alpha) * (L0 - Lbar)^(1-alpha); 

    Result = zeros(size(x));
    
    Result(1) = x(1) - x(4) - delta * x(2) - G ;
    Result(2) = x(6) - rho ;
    Result(3) = M_0 / x(8) - (epsilon_m/epsilon_c) * x(4) * (1 + 1/x(6)) ;
    Result(4) = 1 - x(3) - (epsilon_l/epsilon_c) * x(4) / x(5);
    Result(5) = x(5) - (1 - alpha) * x(7) * x(1) / (x(3) - Lbar);   
    Result(6) = rho + delta - alpha * x(7) * x(1) / x(2);
    Result(7) = x(7) - (theta - 1) / theta;
    Result(8) = x(1) - Omega_0 * x(2)^(alpha) * (x(3) - Lbar)^(1-alpha); 
    
end

