def trapezes(f,a,b,n): h=(b-a)/n L=0 for i in range(n) : L+=(f(a+i*h)+f(a+i*h+h))/2 return L*h