该代码给出了激波管的欧拉一维非定常黎曼问题的精确解。 解析解是通过牛顿-拉夫森法和特征方程来计算的。 将输入设置为左右气体变量后，可以轻松检查不连续性产生的波； 程序在输出中返回解决方案的类型，例如可以是 RCS（意味着对第一个系列产生了扩展，对第二个系列产生了冲击）。 然后用户有机会获得最重要的变量（声速、压力、速度）与激波管长度坐标的关系图。 注意：所有的交互都被忽略了。 输入：
p1：激波管左侧的气体压力rho1：激波管左侧的气体密度u1：激波管左侧粒子的速度p4：激波管右侧的气体压力rho4：激波管右侧的气体密度u4：激波管右侧的粒子速度tol：溶液的耐受性
弗吉尼亚诺塔罗
