A numerical alternate linear-nonlinear method for solving two-dimensional unsteady infiltration of water in unsaturated soils is presented. The results compare favorably with previously published results obtained by conventional numerical methods on similar examples. (AA)