function z = dsge_steadystate_helperz(z0,rhoz,lnz_ss)
options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
z = fsolve(@(z)(1-rhoz)*lnz_ss+rhoz*z,z0,options);
