Trapezoid method for integral calculations