function HW9_14a %Trapezoidal

%Input data
a = 1 ;
b = 2;
n = 91;
h = (b-a)/n;
x = a:h:b;
f = x.*log(x);

%weight=(1,2,2,2,2,...,2,2,2,1)
w(1:n+1)=2;
w(1) = w(n+1) = 1;

%Sum
s = (h/2)*w*f'

end

%n=91, s=0.636301336381243
%n=92, s=0.636301185569557