def riemann(f,a,b,n):
    s=0 # s contient la somme des aires des rectangles
    for i in range(n):
        r=((b-a)/n)*f(a+(i/n)*(b-a)) # r contient l'aire d'un rectangle
        s=s+r
    return s

def f(x):
    return (1-x**2)**0.5 # c'est la fonction arc de cercle


print(2*riemann(f,-1,1,100)) # affiche l'approximation de pi pour 100 rectangles