function fct = integrand(x);

   x(x==0)=NaN;
   fct = log(1+x)./x./(1+x);
   fct(isnan(fct))=1;