G(x,t) = (4πκt)^(−d/2) exp(−|x|²/4κt). A delta-function initial condition spreads as a Gaussian; multiple sources superpose. Click canvas to add heat sources.