Si t'as n rectangles, tu calcules les deux sommes suivantes :
Sn = (b-a)/n * somme(k=0 à n-1, f(a + k(b-a)/n)) et
Tn = (b-a)/n * somme(k=1 à n, f(a + k(b-a)/n)), qui correspondant respectivement aux sommes des aires des rectangles inférieurs et supérieurs. La valeur de l'intégrale est comprise entre ces 2 sommes.
Si tu veux une précision à 10^-p près, tu fais une boucle "tant que Tn-Sn > 10^-p".