methode trapeze

def fonction_a_integrer(x):
    return x**2  # Remplacez cette fonction par celle que vous souhaitez intégrer

def methode_trapezes(a, b, n):
    # Calculer la largeur de chaque sous-intervalle
    h = (b - a) / n

    # Calculer la somme des hauteurs des trapèzes
    somme = 0.5 * (fonction_a_integrer(a) + fonction_a_integrer(b))

    # Ajouter la contribution de chaque trapèze
    for i in range(1, n):
        x_i = a + i * h
        somme += fonction_a_integrer(x_i)

    # Multiplier par la largeur du sous-intervalle
    resultat = h * somme
    return resultat

# Définir les bornes d'intégration
a, b = 0, 1

# Définir le nombre de sous-intervalles
n = 1000

# Calculer l'intégrale à l'aide de la méthode des trapèzes
resultat_integration = methode_trapezes(a, b, n)

# Afficher le résultat
print(f"Résultat de l'intégrale par la méthode des trapèzes : {resultat_integration}")