where u is the dependent variable, f is the source term, and ∇² is the Laplacian operator.
% Create the mesh x = linspace(0, L, N+1); matlab codes for finite element analysis m files hot
The heat equation is:
∂u/∂t = α∇²u