function a = dsge_steadystate_helpera(a0)
  
  options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
 a = fsolve (@(ln(a)) rhoa*ln(a(-1))+epsilona, a0,options);
