All modeling is a gross simplification of the real system. So how much complexity should be in a model. The partial differential equation for advection (flow), dispersion (spreading), and decay (breaking down) is

used as an illustrative example on how the same problem can be solved with different levels of complexity. Three methods will be discussed for solving the equation with the associated drawbacks and benefits. Also the experimentation with a numerical approach hints at a pure analytical solution.