Numerical Simulation of the Buckley-Leverett Problem