function [Y,Ym,Yu,K,Km,Ku,L,Lm,Lu,I,Im,Iu,C,rm,ru,wm,wu,G,Tc,Tss,Tl,Tk,TE,B0,B1,p,s,xi,phi,sigma] = outputhelperfiscal_calibration(x,alpha,alphau,A,B,beta,delta,tauC,tauSS,tauL,tauK,Y0,Ym0,Yu0,K0,Km0,Ku0,L0,Lm0,Lu0,I0,Im0,Iu0,C0,rm0,ru0,wm0,wu0,G0,Tc0,Tss0,Tl0,Tk0,TE0,B00,B10,p0,s0,xi0,phi0,sigma0) 
    
    fun = @(x)steadystatehelperfiscal_calibration(x,alpha,alphau,A,B,beta,delta,tauC,tauSS,tauL,tauK);
    x0 = [Y0,Ym0,Yu0,K0,Km0,Ku0,L0,Lm0,Lu0,I0,Im0,Iu0,C0,rm0,ru0,wm0,wu0,G0,Tc0,Tss0,Tl0,Tk0,TE0,B00,B10,p0,s0,xi0,phi0,sigma0];
    options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10); 
    
    x = fsolve(fun,x0,options);
    
    Y = x(1); 
    Ym = x(2);  
    Yu = x(3);
    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);
    
end