permet le calcul de la valeur approché de l’intégrale de f sur [a;b] en utilisant la méthode des trapèzes, avec a : borne inférieure, b: borne supérieure, n: augmente la précision Dans ce cas, ₀¹∫ x² fait bien 0.333… = 1/3
def trapezes(f,a,b,n): s=0 h=(b-a)/n x=a for i in range(n): xh=x+h s+=h*(f(x)+f(xh))/2 x=xh return s def f(x): return x**2 a=0 b=1 n=1000 print(trapezes(f,a,b,n))