There is some confusion here as to what your variable of integration is to be. In the way you define ‘h’ it would appear to be some variable that varies from 0.308 to 0.478 in ten discrete steps. But how does that relate to the two vectors ‘r’ and ’T’ with their ’n’ values? You would need for 'n' to equal eleven to be compatible with the ten in the definition of ‘h’.
Also the term “+4*sum(f(r(2:2:end),T(2:2:end)))” looks incorrect. It should be “+4*sum(f(r(2:2:end-1),T(2:2:end-1)))”.
