Implementation of the finite difference scheme (explicit) for 1D heat equation: heat_explicit_no_thrills.m (simple .m script file) OR heat_explicit.m (with fancy user-defined functions)
Implementation of the theta method for the 1D heat equation: heat_theta_method.m (See also PDF report version, Quicktime animantion) includes explicit Euler scheme (theta =0), implicit Euler scheme (theta =1), Crank-Nicolson (theta=1/2)