cython

1
from sin_of_square cimport Function, SinOfSquareFunction
2

3
def integrate(Function f, double a, double b, int N):
4
    cdef int i
5
    cdef double s, dx
6
    if f is None:
7
        raise ValueError("f cannot be None")
8

9
    s = 0
10
    dx = (b - a) / N
11

12
    for i in range(N):
13
        s += f.evaluate(a + i * dx)
14

15
    return s * dx
16

17
print(integrate(SinOfSquareFunction(), 0, 1, 10000))
18