function [F]=ss_sol_mod(steady, OMEGA, ETA, TAU_C, R_SS, VARPHI, CHI,BIG_M, XI, GAMMA, G_SS, PAI_SS, PHI, KAPPA, RHO, BETA, MC_SS, PSI, TAU_W)

NF_SS=steady(1);      
NI_SS=steady(2);
CF_SS=steady(3);
C_SS=steady(4);
PI_SS=steady(5);
LAMBDA_SS=steady(6);
THETA_SS=steady(7);
V_SS=steady(8);
WF_SS=steady(9);
TLS_SS=steady(10);

F(1)=C_SS-(OMEGA*CF_SS^ETA+(1-OMEGA)*NI_SS^ETA)^(1/ETA);
F(2)=(1/C_SS)*OMEGA*(C_SS/CF_SS)^(1-ETA)-(1+TAU_C)*LAMBDA_SS;
F(3)=(1/C_SS)*(1-OMEGA)*(C_SS/NI_SS)^(1-ETA)-PI_SS*LAMBDA_SS*(2-1/R_SS);
F(4)=(1/LAMBDA_SS)*CHI*NI_SS^VARPHI-PI_SS+BIG_M*(THETA_SS^(1-XI))*(1-TAU_W)*WF_SS+(1-BIG_M*(THETA_SS^(1-XI)))*GAMMA*WF_SS;
F(5)=G_SS+GAMMA*WF_SS*(1-NI_SS-NF_SS)-TAU_W*WF_SS*NF_SS-TAU_C*CF_SS-((PAI_SS-1)/PAI_SS)*PI_SS*NI_SS-TLS_SS;
F(6)=THETA_SS-V_SS/(1-(1-RHO)*NF_SS-NI_SS);
F(7)=NF_SS*(1-(PHI/2)*(PAI_SS-1)^2)-CF_SS-G_SS-KAPPA*V_SS;
F(8)=WF_SS-PSI*MC_SS-GAMMA*WF_SS*(1-PSI)/(1-TAU_W)-PSI*BETA*(1-RHO)*KAPPA*THETA_SS;
F(9)=NF_SS*RHO-V_SS*BIG_M*(THETA_SS^-XI);
F(10)=(KAPPA/BIG_M)*(THETA_SS^XI)*(1-BETA*(1-RHO))-MC_SS+WF_SS;

