Do some searches on the Finite Volume method and integration approaches. Lots of work is done in CFD/CSD coming up with stable, conservative approaches. It might be a little bit overkill than say, the trapezoid rule, but if you're looking for high-accuracy, that'd be the way to go.