"foam1D.m" is the main program;
"decomt.m" is a subroutine for LU decomposition;
"solvet.m" is a subroutine for backward substitution;
"plots.m" is the code for graphical output.