numeerinenintegrointi.py

# Numeerinen integrointi Pythonilla
# Juhani Kaukoranta 5.2.2019
import scipy.integrate as integrate
import scipy
# f on yksikköympyrän x²+y² = 1  kaari f = y = sqrt(1-x**2)

f = lambda x: scipy.sqrt(1-x**2)

# y integroinnin tulos, err on tarkkuus
y,err = integrate.quad(f,-1,1)
print("yksikköpuoliympyrän ala integroimalla ala -1 --> 1 = ",y," virhe = ",err)
print("yksikköympyrän ala integroimalla -1 --> 1 = ",2*y)
print("teoreettinen arvo = ",scipy.pi)