cython

1
def f(x):
2
    return x ** 2 - x
3

4

5
def integrate_f(a, b, N):
6
    s = 0
7
    dx = (b - a) / N
8
    for i in range(N):
9
        s += f(a + i * dx)
10
    return s * dx
11