function HW9_7c %Midpoint

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

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

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

end

%s =  3.00906003100695