Computational Physics, Physics, Heat Equation, Python, Math.
ℎ² which gives us the ratio 𝑟=𝑘/ℎ², and the simplified form will look look like below: For 1-Dimension import numpy as np from matplotlib import pyplot as plt L = 100 x = np.linspace(0, 1, L) u = np.zeros(L) #Initial Condition u[L//2] = 1 #Boundary Conditions #This is not necessary to declare since