function v = dsge_steadystate_helperv(v0)
  
  options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
 v = fsolve (@(ln(v)) rhov*ln(v(-1))+epsilonv, v0,options);
