Modifier la fonction etat final du cours pour qu'elle affiche un…

Avancement et bilans de matière · Problèmes et outils numériques · niveau 3 (difficile)

Énoncé

Modifier la fonction etat_final du cours pour qu'elle affiche un message d'alerte lorsque le mélange initial est stœchiométrique (à $1 0^{- 9}$ près). La tester sur la synthèse de l'ammoniac stœchiométrique.

Corrigé


def etat_final_v2(n0, coeffs_r, coeffs_p):
    candidats = [n0[i] / coeffs_r[i] for i in range(len(n0))]
    x_max = min(candidats)
    if max(candidats) - x_max < 1e-9:
        print("Mélange stoechiométrique : tous les réactifs s'épuisent.")
    restes   = [n0[i] - coeffs_r[i] * x_max for i in range(len(n0))]
    produits = [c * x_max for c in coeffs_p]
    return x_max, restes, produits

print(etat_final_v2([2.0, 6.0], [1, 3], [2]))
# Mélange stoechiométrique : tous les réactifs s'épuisent.
# (2.0, [0.0, 0.0], [4.0])

Le test « tous les candidats égaux » (max $-$ min $< ε$ ) traduit exactement la condition $\frac{n _{0} ( A )}{a} = \frac{n _{0} ( B )}{b}$ — à epsilon près, réflexe des flottants.