Another simple program demostrates how to apply WENO-LF schemes of 3rd and 5th orders to reconstruct fluxes for 1D Buckley-Leverett equation. // Journal of Computational Physics, 126, pp.202-228 (1996).
Also Method of lines (MOL) have used for time-discretization as well as Lax-Friedrichs flux splitting procedure.
Feel free to ask any questions.
Some mistakes are corrected
Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.