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