# -*- 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
u=#A COMPLETER
u_approx=np.arange(0,T,dt)
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$'),4)
plt.show()