急求！！WENO格式计算某些算例为什么会震荡的很厉害？

r1b2i3

我在用matlab代码WENO-JS5格式计算一维欧拉方程算例，初值为(rhoL,uL,pL)=(0.445,0.698,0.3528) if 0<x<0.5; (rhoR,uR,pR)=（0.5，0，0.571），时间为0.14 s,时间上采用3阶RK方法. 计算后密度图像为
请问是甚么原因？是边界条件不对吗？还是其他原因？但是我在看论文时，论文上面没有震荡。谁有正确的代码？顺便问一下，有二维双马赫反射问题或二维欧拉方程的matlab代码吗？

pswpswpsw

你肯定代码有问题。。自己DEBUG吧。