function y = F(x)

betaa = 0.98;
varphi = 0.5;  
sigma = 2;
omega_T = 1.455;
omega_N = 1.455;
omega_C = 1.455;
chi = 0.35;
delta = 0.025; 
alpha_T = 0.35;
alpha_N = 0.25;
alpha_C = 0.35;
gammaa_T = 0.05;
phi_T = 4.6;
phi_N = 9.3;
phi_C =10.3;
v_C = -0.014; 
v_D = 0.077;
rho_C = 0.88;
rho_T = 0.87;
rho_N = 0.84;
psi_D = 0.077;
psi_C = -0.014;

reer = x(1);
KL_N = x(2);
Lss_N = x(3);
Pss_N = x(4);

Dss = 0.12;
Rss = 1/betaa;
muss_T = (1 - betaa*(1-delta))/betaa;
muss_N = muss_T;
muss_C = muss_T;
TBss = (Rss-1)*Dss;
Pss_C = 0.69;

% Commdity sector: {Kss_C/Lss_C and Wss_C}
KL_C = (alpha_C*Pss_C/muss_C)^(1/(1-alpha_C));
Wss_C = (1-alpha_C)*Pss_C*(KL_C^alpha_C);

% Tradable sector: {Kss_T/Lss_T and Wss_T}
KL_T = ((alpha_T/muss_T)*(((gammaa_T*muss_T)/(Pss_C*alpha_T))^gammaa_T))^(1/(1 - alpha_T - gammaa_T));
Wss_T = (1 - gammaa_T - alpha_T)*((gammaa_T*muss_T/Pss_C*alpha_T)^gammaa_T)*KL_T^(alpha_T + gammaa_T); 

% Nontradable sector: {solve numerically for reer, Kss_N/Lss_N, Lss_N and Pss_N}
y(1) = reer - (chi^varphi + ((1-chi)^varphi)*(Pss_N^(1-varphi)))^(1/(1 - varphi));
y(2) = muss_N - alpha_N*Pss_N*(KL_N^(alpha_N - 1));
y(3) = reer*(Lss_N^(omega_N - 1)) - (1-alpha_N)*Pss_N*(KL_N^alpha_N);
y(4) = Pss_N - ((1-chi)/chi)*(((((gammaa_T*muss_T)/(Pss_C*alpha_T))^gammaa_T)*(KL_T^(alpha_T + gammaa_T))*((Wss_T/reer)^(1/(omega_T - 1))) - delta*KL_T*((Wss_T/reer)^(1/(omega_T - 1))) - delta*KL_N*Lss_N - delta*KL_C*((Wss_C/reer)^(1/(omega_C - 1))) - (TBss - Pss_C*(((KL_C^alpha_C)*(Wss_C/reer)^(1/(omega_C-1))) - (gammaa_T/Pss_C)*((((gammaa_T*muss_T)/(Pss_C*alpha_T))^gammaa_T)*(KL_T^(alpha_T + gammaa_T))*((Wss_T/reer)^(1/(omega_T - 1)))))))/((KL_N^alpha_N)*Lss_N))^(1/varphi);

end