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