function out = IM(l,w)

out = quad(@MGrnd,w*exp(-l),w);
