matlab *.m files to solve u_t + (f(u))_x = 0.
These programs are for the equation u_t + (f(u))_x = 0.
the naive code
which gets the speed wrong because the equation isn't in
conservation form. Here's
where I show that it's not working as desired.
the demos for the conservative schemes.
You need the conservative
time-stepping scheme; the
initial data files:
shock2.m; the numerical flux functions for
Lax-Wendroff; and the files that find the