我阅读了带有源项的双曲守恒定律的 ADER WENO 有限体积方案。我想在 Matlab 编码中实现这个方法。在计算数值通量的过程中，我在编码方面遇到了一些问题。请帮我Matlab编码

$$u_t+f(u)_x=S(u,x,t),\; u(0,x)=u_0(x),\; x\in \mathbb{R}, \;t>0,\;u(t,-\infty)=u(t,\infty)=0.$$