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.