function HW9_24b %Simpson

%Input data
re = 0.308 ;
ro = 0.478;
n = 10;
h = (ro-re)/n;
r = re:h:ro;
T = zeros(1,11);

T(1) = 640;    T(2) = 794;     T(3) = 885;   T(4) = 943;
T(5) = 1034; T(6) = 1064;   T(7) = 1114; T(8) = 1152;
T(9) = 1204; T(10) = 1222; T(11) = 1239;

f = T.*r;

%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'

%Ans
int_r = (ro^2-re^2)/2;
Ans = s/int_r

end

%s =  70.4641046666667
%Ans =  1054.69397794741