cython

1
from cython.cimports.sin_of_square import Function, SinOfSquareFunction
2

3
def integrate(f: Function, a: float, b: float, N: cython.int):
4
    i: cython.int
5

6
    if f is None:
7
        raise ValueError("f cannot be None")
8

9
    s: float = 0
10
    dx: float = (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