数值方法期末作业有限体积WENO+北大高彬

高同学讲得很细也很清楚，报告时的自信从容、有条不紊、对所做格式熟悉到信手拈来。考虑到班里的很多同学对偏微分方程数值数值解不熟悉，建议高同学能够把自己代码的注释进一步的完善，使得其他同学以后再遇到这样的问题的时候能够方便地参考高同学的代码。

---下面是我看了高同学的代码之后抛砖引玉，增加的一些注释。 ----

在北太天元上运行得到数值结果如上图所示。

另外，可以做一些小改可以使得计算得更快，

把 WENO函数 修改的可以处理向量参数，

实际山就是把 ^2 改成 .^2, 把 * 改成 .*

然后把下面的for循环

%for i=1:N %重构

%uL(i+1)=WENO(uL2(i),uL1(i),u(i),uR1(i),uR2(i));

%uR(i+1)=WENO(uR3(i),uR2(i),uR1(i),u(i),uL1(i));

%end

改成

uL(2:end)=WENO(uL2,uL1,u,uR1,uR2);

uR(2:end)=WENO(uR3,uR2,uR1,u,uL1);