cython

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

4

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