function z = dsge_steadystate_helperz(z0, rhoz)

options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);

z = fsolve(@(z)) (1-rhoz)*ln(z_ss)+rhoz*z(-1),z0,options);
