import numpy as np
import matplotlib.pyplot as plt
# données du pb
Lambda = 0.037
cp = 1500
rho = 1.325
L = 1  # épaisseur isolant
t_max = 20000  # temps de fin d'intégration en secondes
N_t = 100 # nombres d'intervalles ds le tps
N_x = 10 # nombres d'intervalles ds l'espace
T_int = 20
T_ext = 5
K = Lambda/cp/rho # diffusivité th

# discrétisation de  l'espace et du tps


Temp = np.zeros((N_t + 1, N_x + 1))

# initialisation de la température
# conditions initiales

# conditions aux limites


# calcul des températures aux différents instants







plt.figure(figsize=(8,8))
for n in [0, 30, 60, 90] :
    plt.plot(list(range(N_x+1)), Temp[n, :], label=str(n))
plt.legend()
plt.show()