∂u/∂t = D∇²u — heat/diffusion equation on a 2D domain. Watch concentration spread from initial conditions. Click to add sources. Choose boundary shapes (Dirichlet = 0 or periodic). Explicit finite-difference scheme with CFL stability control.