def trapezoidal(f, a, b, n): h = (b - a) / n s = 0.5 * (f(a) + f(b)) for i in range(1, n): s += f(a + i*h) return h * s