C'est une méthode de Newton-Cotes : On subdivise l'intervalle d'intégration [a ;b] en n sous intervalles de longueur égale h.
h=(b-a)/n
Sur chaque sous intervalle [xi;xi+1] on remplace f par son polynôme d'interpolation de degré 1 (sa corde) passant par les points (xi,f(xi)) et
(xi+1, f(xi+1)).
Algorithme : Intégrale f(x)dx entre [a,b] = h[ (f(a)+f(b))/2 + SOMME de i=1 à n-1 de f(a+i*h)]