Transient 1D diffusion in k and t domains

The main program has two functions.

First, it finds the transfer function F(k) in the wavenumber domain for specified diffusivity kappa, time step dt, and grid spacing dx.

Then code solves the evolution equation in the time domain on a limited domain of length Lx, over time period Lt, using specified boundary conditions and initial conditions.

Users can compare predictions of stability and accuracy from the wavenumber analysis against the actual evolution in the time domain.

The m-file solve.m (Same Old Solver) solves the resulting system of linear equations at every time step.

Main program
Same Old Solver