cython

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

4

5
def integrate_f(a: cython.double, b: cython.double, N: cython.int):
6
    i: cython.int
7
    s: cython.double
8
    dx: cython.double
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