A Galerkin-type finite element method is employed to solve the quasilinear partial differential equations of transient seepage in saturated-unsaturated porous media. The resulting computer program is capable of handling nonuniform flow regions having complex boundaries and arbitrary degrees of local anisotropy. Flow can take place in a vertical plane, in a horizontal plane, or in a three-dimensional system with radial symmetry. An arbitrary number of seepage faces can be considered simultaneously, and the positions of the exit points on these boundaries are adjusted automatically during each time step. Two examples, one of seepage through an earth dam with a sloping core and horizontal drainage blanket, and the other of seepage through a layered medium cut by a complex topography, are included. These examples indicate that the classical concept of a free surface is not always applicable when dealing with transient seepage through soils.