Computational Physics, Physics, Heat Equation, Python, Math.
matplotlib.animation import FuncAnimation plate_length = 50 max_iter_time = 1000 alpha = 2.0 delta_x = 1 # Calculated params delta_t = (delta_x ** 2)/(4 * alpha) gamma = (alpha * delta_t) / (delta_x ** 2) # Initialize solution: the grid of u(k, i, j) u = np.empty((max_iter_time, plate_length,