function HW9_14c %Midpoint %Input data a = 1 ; b = 2; n = 128; h = (b-a)/(n+2); x = a+h:2*h:b-h; f = x.*log(x); %weight=(1,1,1,1,1,...,1,1,1,1) w(1:n/2+1)=1; %Sum s = (2*h)*w*f' end %n=128, s=0.636287525399937