It is shown how delta shock waves which consist of Dirac delta distributions and classical shocks can be used to construct non-monotone solutions of the Buckley–Leverett equation. These solutions are interpreted using a ...