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