import numpy as np
import matplotlib.pyplot as plt

# données du problème
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 = 5 # nombres d'intervalles ds l'espace
T_int = 20
T_ext = 5
K = [Instruction 1] # diffusivité th

# discrétisation de  l'espace et du tps
dx = [Instruction 2.1]
dt = [Instruction 2.2]
Temp = np.zeros((N_t + 1, N_x + 1))

# initialisation de la température
# conditions initiales
Temp[0,0] = [Instruction 3.1]

for i in range (1,N_x + 1):
    [Instruction 3.2]

# conditions aux limites
for n in range (1,N_t + 1):
    [Instruction 3.3]
    [Instruction 3.4]

# calcul des températures aux différents instants
for n in [Instruction 4.1]:
    for i in [Instruction 4.2]:      # ni le 1er ni le dernier point
        [Instruction 4.3]

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