function HW9_7a %Trapezoidal

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

%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

%s =  3.15947567425982