求助 激波管编程

baiyang1985

激波管问题：设有一充满完全理想气体的无穷长直管，绝热指数1.4。初始时刻物理量分布为：x<=0时，(ρ(x,0),u(x,0)，p(x,0))=(1.0,0.75,1.0)
x>0时，(ρ(x,0),u(x,0)，p(x,0))=(0.125,0.,0.1)
用一阶FVS格式，一阶Roe格式或NND格式计算t=0.2时的数值解。（计算域取(-0.5,0.5)）。希望好心人帮帮忙，谢谢了。