function HW9_24a %Trapezoidal

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

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

end

%s =  70.5064120000000
%Ans =  1055.32722646310