Numerische Lösung des Wärmeleitungsproblems

Es ist notwendig, das Wärmeleitungsproblem auf einem Intervall zu lösen.

Bei der Lösung wurde ein explizites Schema verwendet.

N = 10 # maximale Anzahl der Schritte in xK = 10 # maximale Anzahl der Schritte in tl = 1 # Wert von x am rechten Randh = l / N # Schrittweite des Gitters in xT = 1 # maximaler Zeitwert t am rechten Randt = T / K # Schrittweite des Gitters in der Zeit# definiere das Gitterx_i = np.arange(0, N, h) # Werte an den Knoten in xt_j = np.arange(0, K, t) # Werte an den Knoten in tr_j = len(t_j) # Anzahl der Knoten in tr_i = len(x_i) # Anzahl der Knoten in xw_h_t = np.zero