from Eq. (92)
with flux at
Naive scheme: stencil diagram and (disastrous) test result.
The oscillations already seen in Fig. 4
After some time the numerical result does not have the faintest
resemblance to true solution.