守恒方程的数值解，该程序很简单，没有处理边界条件-Numerical solution of the conservation equations, the program is very simple, does not deal with the boundary conditions