function out = IM(l,w)
if l==0
    out = 0;
else
    MRes = (exp(0.5.*l)-1)./sqrt(w);
    out = quad(@MGrnd1,w*exp(-l),w)+MRes;
end
