function HW9_14b %Simpson

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

%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

%n=6, s=0.636297500790914