由四位同学完成。

原理 步骤顺序 网格化 设定初值 求解网格面两侧的值 计算网格中波速 计算网格面上流的近似值 更新$U$ 反算$P$ 消除真空 重复循环