### matlab *.m files to solve u_t + (f(u))_x = 0.

These programs are for the equation u_t + (f(u))_x = 0.

Here's
the naive code
which gets the speed wrong because the equation isn't in
conservation form. Here's
the code
where I show that it's not working as desired.

Here are
the demos for the conservative schemes.
You need the conservative
time-stepping scheme; the
initial data files:
ID.m,
shock1.m, and
shock2.m; the numerical flux functions for
explicit upwind,
Lax-Friedrichs, and
Lax-Wendroff; and the files that find the
location
of the
shock region.