cumsum doesn't sound very accurate. Should we still get the integrate.quad solution somewhere?

Also don't we need to correct for the distance between points.