Solution numérique du problème de conduction thermique

Il est nécessaire de résoudre le problème de conduction thermique sur un segment

Lors de la résolution, j'ai utilisé un schéma explicite

N = 10 # nombre maximal de pas en xK = 10 # nombre maximal de pas en tl = 1 # valeur de x à la frontière droiteh = l / N # pas de la grille en xT = 1 # valeur maximale du temps t à la frontière droitet = T / K # pas de la grille en temps# définissons la grille x_i = np.arange(0, N, h) # valeurs aux nœuds en xt_j = np.arange(0, K, t) # valeurs aux nœuds en tr_j = len(t_j) # nombre de nœuds en tr_i = len(x_i) # nombre de nœuds en xw_h_t = np.zero