perturbed trapezoid formula