Romerg, mentioned by @rafael.guerra above, is surprisingly accurate for integrating smooth functions but it requires the values being integrated be evenly sampled over the independent variable.
You can always fall back to trapezoidal integration if the points are not evenly sampled, but it’s low accuracy.
I"m not sure how the quadrature weighting method would compare in performance and accuracy. Need to try it out.