function HW9_7b %Simpson

%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,4,2,4,2,...,4,2,4,1)
w(1:n+1)=2;
for j = 1:n/2
	w(2*j) = w(2*j)*2;
end
w(1) = w(n+1) = 1;

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

end

%s =  3.10933712650887