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

h1 h2 h3 h4 print("???") $$aaa^aa$$