In this paper, a numerical solution is presented for one-dimensional unsaturated flows in the subsurface. Water flow in the subsurface, however, is highly nonlinear and in most cases, exact analytical solutions are impossible. The method of reference-operators has been used to formulate a discrete model of the continuum physical system. Many of the standard finite difference methods and also th...