"""
Thème 4 : calcul intégral
"""

###Q1


# Faites un dessin à la main

# rem : pas de conclusion sur d'éventuelles considérations
# de sous-estimations de l'aire par les rectangles.

## Q2

def Riemann(f,a,b,n):
    h = (b-a)/n   # largeur du rectangle/ pas de la méthode des rectangles
    S  = 0        # initialisation du calcul de somme
    for k in range(n):
        S+=h*f(a+k*h)
    return S

## Q3a)

def f(x):
    return x
a,b=0,1

for n in [10,20,30,40,100,150,1000]:
    An = Riemann(f,a,b,n)
    print(f"Avec {n} rectangles, l'aire estimée est : {An}")