function x = dsge_steadystate_helperx(x0)
  
  options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
 x = fsolve (@(ln(x)) rhox*ln(x(-1))+epsilonx, x0,options);
