function F = steadystatehelperfiscal_calibration(x,alpha,alphau,A,B,beta,delta,tauC,tauSS,tauL,tauK)

    %Y = x(1), Ym = x(2), Yu = x(3), A = x(...), B = x(...), 
    %K = x(4), Km = x(5), Ku = x(6), L = x(7), Lm = x(8), Lu = x(9), 
    %I = x(10), Im = x(11), Iu = x(12), C = x(13), rm = x(14), ru = x(15), 
    %wm = x(16), wu = x(17), G = x(18), Tc = x(19), Tss = x(20), 
    %Tl = x(21), Tk = x(22), TE = x(23), B0 = x(24), B1 = x(25), p = x(26),
    %s = x(27), xi = x(28), phi = x(29), sigma = x(30)
    
       F(1) = x(14) - (1/beta - 1 + delta) / (1-tauK); 
       F(2) = x(15) - (1/beta - 1 + delta); 
       F(3) = x(16) - A * (1-alpha) * (1-tauC) / (1+tauSS-tauC) * ((alpha * A * (1-tauC) / x(14))^(1/(1-alpha)))^alpha; 
       F(4) = x(17) - B * (1-alphau) * ((alphau * B * (1-x(26)*tauC*x(27)) / x(15))^(1/(1-alphau)))^alphau; 
       F(5) = (1-tauL) * x(16) - x(24)*(x(8)+x(9))^x(28) / (x(13)^(-x(30))); 
       F(6) = x(17) - (x(24)*(x(8)+x(9))^x(28) / (x(13)^(-x(30))) + x(25) * x(9)^x(29) / (x(13)^(-x(30)))); 
       F(7) = x(2) - A * ((alpha * A * (1-tauC) / x(14))^(1/(1-alpha)))^alpha * x(8); 
       F(8) = x(3) - B * ((alphau * B * (1-x(26)*tauC*x(27)) / x(15))^(1/(1-alphau)))^alphau * x(9); 
       F(9) = x(1) - (x(2) + x(3));
       F(10) = x(5) - (alpha * A * (1-tauC) / x(14))^(1/(1-alpha)) * x(8); 
       F(11) = x(6) - (alphau * B * (1-x(26)*tauC*x(27)) / x(15))^(1/(1-alphau)) * x(9); 
       F(12) = x(4) - (x(5) + x(6));
       F(13) = x(7) - (x(8) + x(9));
       F(14) = x(11) - delta * x(5); 
       F(15) = x(12) - delta * x(6); 
       F(16) = x(10) - (x(11) + x(12));
       F(17) = x(13) + x(10) - ((1-tauL) * x(16) * x(8) + (1-tauK) * x(14) * x(5) + x(17) * x(9) + x(15) * x(6)); 
       F(18) = x(1) - (x(13) + x(10) + x(18)); 
       F(19) = x(19) - (tauC * ((x(2) - x(16)*x(8)) + x(26)*x(27) * (x(3) - x(17)*x(9)))); 
       F(20) = x(20) - tauSS * x(16) * x(8); 
       F(21) = x(21) - tauL * x(16) * x(8); 
       F(22) = x(22) - tauK * x(14) * x(5);
       F(23) = x(23) - ((1-x(26))*tauC*(x(3)-x(17)*x(9)) + tauSS*x(17)*x(9) + tauL*x(17)*x(9) + tauK*x(15)*x(6)); 
       F(24) = x(3) / x(1) - 0.448757372602233;
       F(25) = x(9) / x(7) - 0.405352297070370;
       F(26) = x(13) / x(1) - 0.540091818948535;
       F(27) = x(10) / x(1) - 0.179062103298792;
       F(28) = x(18) / x(1) - 0.280846077752673;
       F(29) = x(23) / x(1) - 0.244603762598837;
       F(30) = x(7) - 0.986424426012063;

end