Jacobi Solver for the Unsteady Heat Equation

Although this is a trivial exercise, many people find puzzling the correct procedure to quantify the error produced by a finite difference approximation. The following scripts are a simple demonstration on how to measure the norm errors for 1d, 2d and 3d finite difference solvers. Here, by using the classical heat equation with a Jacobi scheme we simply demonstrate the computation of the L1, L2 and Linf norm error, for each case.