import numpy as np
import matplotlib.pyplot as plt

## Parametres du système et de la réaction A + B = C (A : anhydride ; B : eau ; C : acide acétique)
# Quantités de matière initiales en mol
nA0 = 
nB0 = 
# Température initiale en K
T0 = 
# Enthlapie standard de réaction en J/mol
deltarH0 = 
# Capacité thermique du réacteur en J/K
Créac = 250 
# Capacités thermiques molaires en J/(mol.K)
CmA = 189.7 
CmB = 75.4
CmC = 119.3

# Capacité thermique totale en J/K
def Ctot(ksi):
    return()

# Facteur pré-exponentiel en 1/min et énergie d'activitation en J/mol
KA = 1.68e7
EA = 50.5e3
# Constante des gaz parfaits en J/(mol.K)
R = 8.314

# Loi d'Arrhénius
def k(T):
    return()

## Résolution numérique
# Nombre d’iterations, instant final et pas de temps en min
n = 1001
tf = 
dt = tf/(n-1)

# Listes des instants
t = 
# Création de listes pour l'avancement et pour la température
av = [0 for i in range(n)]
Temp = [0 for i in range(n)]

# Initialisation
av[0] = 
Temp[0] = 
# Resolution par la methode d’Euler
for i in range(n-1):
    dav =                             # variation de l'avancement sur l'intervalle [ti,ti+1]
    av[i+1] = 
    Temp[i+1] = 

# Representation graphique de la température en fonction du temps