blob: ad4c8540bd006adc7f21d8251a888b2b1c215d53 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
from sin_of_square cimport Function, SinOfSquareFunction
def integrate(Function f, double a, double b, int N):
cdef int i
cdef double s, dx
if f is None:
raise ValueError("f cannot be None")
s = 0
dx = (b - a) / N
for i in range(N):
s += f.evaluate(a + i * dx)
return s * dx
print(integrate(SinOfSquareFunction(), 0, 1, 10000))
|
