function tau = dsge_steadystate_helper(tau0,m,pi)
  
  options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
 tau = fsolve (@(tau) tau=(m/m(-1))*pi, tau0,options);
