import numpy as np
from scipy.optimize import bisect

# Parametres ===================================
# Constante des gaz parfaits en J/(mol.K)
R = 8.314
# Enthalpie standard de réaction en J/mol
deltarH0 = -114.7e3 
# Entropie standard de réaction en J/(mol.K)
deltarS0 = -245.9
# Pression standard en bar
P0 = 1
# Pression totale en bar 
P = 1
# Température en K
T = 373

def ZeroDicho(f,a,b,eps):
    """Recherche un zero de la fonction f entre a et b par dichotomie, et renvoie une valeur approchée à eps près
    La fonction f est supposée continue, et vérifie f(a) * f(b) <= 0"""








# Fonction à annuler pour déterminer le rendement r
def f(r):
    return()

# Résolution par dichotomie puis avec la fonction bisect
r_d=ZeroDicho()
r_b=bisect() 

# Affichage des résultats
print("rendement par dichotomie :", r_d)
print("rendement par bisect :", r_b)