import numpy as np
import matplotlib.pyplot as plt

# Données relatives au problème A COMPLETER
R = #J/mol/K
n0 = #mol
Delta_rH0 = #J/mol
Ea = #J/mol/K
A = #min-1
T0 = #K
Ceau = #J/g/K
meau = #g

# Fonction constante de vitesse k(Temp)
def k(Temp):
    A COMPLETER

# Pas de temps
dt = A COMPLETER # min

# Méthode d'Euler
t_i = A COMPLETER # min instant initial
T_i = A COMPLETER # K température initale
xi_i = A COMPLETER # mol avancement initial

ls_t = [t_i] # liste des instants
ls_T = [T_i] # liste des températures
ls_xi = [xi_i] # liste des avancements

xi = xi_i # initialisation de xi
t = t_i # initialisation de t
T = T_i # initialisation de T

while xi < 0.99*n0 :
    t = t + dt
    ls_t.append(t)
    dxi = A COMPLETER
    xi = xi + dxi
    ls_xi.append(xi)
    dT = A COMPLETER
    T = T + dT
    ls_T.append(T)

# Tracé de T = f(t)
A COMPLETER