File: integrate.pyx

package info (click to toggle)
cython 0.29.21-3
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 14,848 kB
  • sloc: python: 67,537; ansic: 14,505; cpp: 1,340; xml: 1,031; makefile: 367; lisp: 206; sed: 11; sh: 7
file content (14 lines) | stat: -rw-r--r-- 388 bytes parent folder | download | duplicates (4)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
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))