function a = dsge_steadystate_helpera(a0,gamm,c,lambd,m)
options = optimset('Display','Final','TolX',1e-10,'TolFun',1e-10);
a = fsolve(@(a)lambd*c^(1/gamm)*(c^((gamm-1)/gamm)+e^(1/gamm)*(m^((gamm-1)/gamm))),a0,options);
