Numerical methods are among the most powerful tools available for solving problems in subsurface hydrology. Our interest in these methods developed because of our need to solve complicated natural hydrologic problems without resorting to simplifications that are mathematically necessary but physically unrealistic.