import numpy as np
import matplotlib.pyplot as plt
# données du pb
Lambda = # à remplir 
cp = 
mu = 
L =  # longueur du barreau
t_max = 20000  # temps de fin d'intégration en secondes
N_t =  # nombres d'intervalles ds le tps
N_x =  # nombres d'intervalles ds l'espace
T_0 = 20
T_1 = 80
K = # [Instruction 1]

# 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]:      
      #[Instruction 4.3]
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()