function res = f_Yashiv_2022_hours_test(xstar)
 nt = xstar(1);
 ut = xstar(2);
 it = xstar(3);
 ct = xstar(4);
 pit = xstar(5);

 xt = xstar(6);
 mt = xstar(7);
 ft = xstar(8);
 at = xstar(9);
 kt = xstar(10);

 netyt = xstar(11);
% gt = xstar(12);
 dgmt = xstar(12);
 dgnt = xstar(13);
 dgkt = xstar(14);

 invt = xstar(15);
 mct = xstar(16);
 fkt = xstar(17);
 wrt = xstar(18);
 mpt = xstar(19);

 fnt = xstar(20);
 ht = xstar(21);
 mcwt = xstar(22);
 mohnt = xstar(23);

beta = 0.99;   
deltan = 0.126;  
deltak = 0.024;    
alpha = 0.66;   
epsilon = 11;     
gamma = 0.44;  
chi = 1; 
phi = 0.5;      
zeta = 120;    
e = 1;


res(1) = 1 - nt-ut;

res(2) = 1/it - beta;

%res(3) = xt - 0.7214;
res(3) = deltan*nt-xt/(1-xt)*ut;

res(4) = mt - xt/(1-xt)*ut; 

res(5) = ft - at*(nt*ht)^alpha*kt^(1-alpha); 

res(6) = netyt - ft+(e/2*(deltan)^2)*ft;

%res(7) = gt - (e/2*(deltan)^2)*ft;

res(7) = dgmt -  ft/nt*e*(deltan);

res(8) = dgnt + e*ft/nt*(deltan)^2 - alpha*(e/2*(deltan)^2)*ft/nt;

res(9) = dgkt - (1-alpha)*(e/2*(deltan)^2)*ft/kt;

res(10) = fkt - (1-alpha)*ft/kt;

res(11) = deltak*kt - invt;

%res(13) = kt - 22;
res(12) = it - (mct*(fkt-dgkt)+(1-deltak));

%res(14) = nt - 0.9511;
res(13) = mct*dgmt - mct*(fnt-dgnt)+wrt*ht-(1-deltan)*1/it*(mct*dgmt);

res(14) = ct+invt - netyt;

res(15) = (1-epsilon)/zeta + epsilon/zeta*mct;

res(16) = wrt*ht - gamma*mct*(fnt-dgnt)-(1-gamma)*(chi/(1+phi)*ct*ht^(1+phi)+xt/(1-xt)*gamma/(1-gamma)*(mct*dgmt));

res(17) = pit - 1;

res(18) = at - 1;

res(19) = mpt - 1;

res(20) = fnt - alpha*ft/nt;

res(21) = ht - (mct/(chi*ct)*at*alpha*kt^(1-alpha)*(alpha*nt^(alpha-1)-(alpha-2)*nt^(alpha-3)*(e/2)*mt^2))^(1/(phi-alpha+1));

res(22) = mcwt - mct*dgmt;

res(23) = mohnt - at*alpha*kt^(1-alpha)*(alpha*nt^(alpha-1)-((alpha-2)*nt^(alpha-3))*(e/2)*mt^2);

end

