integrate1.pyx 177 Bytes
Newer Older
1 2
# cython: language_level=3

3 4 5
def f(x):
    return x**2-x

6

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