function Result = Function15_02_SS(x,PAR,EXO)

alpha     = PAR.alpha;
epsilon   = PAR.epsilon;
delta     = PAR.delta;
mu        = PAR.mu;
rho       = PAR.rho;
Z_0       = PAR.Z_0;
t_L       = EXO.t_L;
t_C       = EXO.t_C;
G         = EXO.G;

% Computes the steady-state equilibrium
% x(1)     is W
% x(2)     is R
% x(3)     is C
% x(4)     is K
% x(5)     is Y
% x(6)     is I
% x(7)     is L
%
% Result(1) =  W * L - (1 - alpha) * Y ;
% Result(2) =  (R - rho) * (1 + t_C) * C - epsilon * mu * (rho + mu) * K ; 
% Result(3) =  Y - C - delta * K - G ;
% Result(4) =  Y - Z_0 * K^alpha * L^(1-alpha) ;
% Result(5) =  R + delta - alpha * Y / K ;
% Result(6) =  I - delta * K ;
% Result(7) =  W * (1 - t_L) * (1 - L) - ((1 - epsilon)/epsilon) * (1+t_C) * C ;

    Result = zeros(size(x));
    Result(1) =  x(1) * x(7) - (1 - alpha) * x(5) ;
    Result(2) = (x(2) - rho) * (1 + t_C) * x(3) - epsilon * mu * (rho + mu) * x(4) ; 
    Result(3) =  x(5) - x(3) - delta * x(4) - G ;
    Result(4) =  x(5) - Z_0 * x(4)^alpha * x(7)^(1-alpha) ;
    Result(5) =  (x(2) + delta) * x(4) - alpha * x(5) ;
    Result(6) =  x(6) - delta * x(4) ;
    Result(7) =  x(1) * (1 - t_L) * (1 - x(7)) - ((1-epsilon)/epsilon) * (1+t_C) * x(3) ;
    
end
