function HW9_24c %Midpoint

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

T(1) = 794; T(2) = 943; T(3) = 1064; T(4) = 1152; T(5) = 1222;

f = T.*r;

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

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

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

end

%s =  70.3794900000000
%Ans =  1053.42748091603