import matplotlib.pyplot as plt

from math import exp
# Modelisation de l'hydrolyse de l'anhydride acetique dans des conditions
# adiabatiques
## Systeme etudie
# 1 mol d'anhydride acetique (102 g) et 10 mol d'eau (180 g) sont introduits
# dans un recipient calorifuge sous agitation magnetique.
# On cherche a determiner l'evolution de la temperature en fonction du temps
## Donnees (variables globales et unités)
T =             #
Cp =            #
Delta_rH =      #
dt =            #
V =             #
## Constante cinetique (loi d'Arrhenius)
def k(T):
    # Renvoie la constante cinetique a la temperature T (min-1)
    return
TT = [T]    # Initialisation de la liste des
ksi = 0     # Avancement de la reaction
t = 0       # Initialisation du temps
tt = [t]    # Initialisation de la liste des valeurs de temps
c =         # Concentration initiale d'anhydride (mol/L)

## Methode d'Euler - lorsque le temps augmente de dt...

while ksi <  :
    # On s'arrete lorsque 99.9% de l'anhydride acetique est hydrolysé
    dksi =           # increment de l'avancement - réaction d'ordre 1 (eau en large excès) (bilan de matiere)
    dT =             # increment de temperature - réacteur adiabatique (bilan energetique)
    # Incrementation des variables
    ksi +=    #
    T +=      #
    t +=      #
    c -=
    # Stockage des variables
    TT.append(T)
    tt.append(t)

## Trace du graphe

plt.title('      ')
plt.xlabel('      ')
plt.ylabel('        ')
plt.plot(  ,  )
plt.show()
plt.savefig('Temperature.png')