The input file to use, with typical numerical values, is:
The exact solution can be found from converting to a coordinate system moving with the average value of u and then solving the characteristic solution using the Newton-Raphson method. A function that does so is:
The following program solves the nonconservative problem:
The above programs produce a plot file burgnc.plt that can be plotted on Unix using:
| ~dommelen/bin/xyplot burgnc.plt 2 ' ' ' ' ' ' 2 '0,100,0,100,.4' |
This creates a printable plot file dump.ps. Replace the first "2" by a 1 to get a screen plot if you are using X-windows. Instead, feel free to change the program to make burgnc.plt consistent with your own graphics program.