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