function e = dsge_steadystate_helper(e0,m,c,r,gamm)
options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
 e = fsolve(@(e)(m/c)*(1-1/r)^gamm,e0,options);
