Add trapezoid numerical integration