clear;clc;a=4;b=2;M=3;f=@(z)quadl(@(x)1./((1+(z-x)/a).*(1+x/b)).^M,0,z);quadl(@(z)arrayfun(@(z)log(1+z)*f(z),z),0,1e8)
ans =
3.1092
