# -*- coding: utf-8 -*-
"""
    Concours Centrale-Supélec
    Oral de Physique-Chimie 2 MP
"""
import numpy as np
import matplotlib.pyplot as plt

#exemple d'application numérique
e=1 #épaisseur en mm
T=1 #durée de l'expérience en s
dt=0.01 #en s
Diff=.5 #diffusivité en mm2/s
   
tps=np.arange(0,T,dt)
# A compléter... l'expression de u

 
plt.plot(tps,u,'b') #solution exacte"
plt.plot(tps,u_approx,'r') #solution approchée
plt.axis([0,T,0,1])
plt.xlabel('$temps (s)$',fontsize=14)
plt.ylabel('$Température \, réduite$',fontsize=14)
plt.legend(('$Courbe \, analytique$','$Courbe \, approchée$'))
plt.show()