I am encoutering a problem when trying to perform the numerical integration of a function I have obtained from interpolation that returns a list. The problem is just with NIntegrate, because Plot works well.
dataX = Table[{n, {n, 2 n}}, {n, 0, 10}];
dataF = Interpolation[dataX, InterpolationOrder -> 1];
dataFi[x_?NumericQ, i_] := dataF[x][[i]];
Plot[dataF[x], {x, 0, 10}]
NIntegrate[dataF[x], {x, 0, 10}] (*Returns NIntegrate::inum*)
NIntegrate[Table[dataFi[x, i], {i, 1, 2}], {x, 0, 10}] (*Returns result with several Part::partw warnings*)
Table[NIntegrate[dataFi[x, i], {x, 0, 10}], {i, 1, 2}] (*Returns result and one Part::partw warning*)
I think I have done this in the past with previous versions and it was working, now with 13.3 it does not. Is it a bug or is this intended? Possible workarounds?