A hybrid finite analytic method, a numerical technique, was used to solve lateral subsurface storm flow described by the extended Dupuit-Forchheimer equation. In subregions the local linearized one-dimensional Dupuit-Forchheimer equation was solved analytically in space and discretized in time by a simple difference formula. This is called a hybrid finite analytic approach. The system of algebraic equations obtained by the hybrid finite analytic approach can approximately preserve the overall nonlinear effect because the term (∂H/∂X)2 and the coefficient of diffusion are treated as constants only in the local regions. A four-point numerical formula can provide results that are stable and acceptably accurate without complex calculation and without small time steps. The comparison between steady state profiles of the water table provided by this study with profiles provided by previous investigations shows that the method performs very well in the analysis of lateral subsurface storm flow. ¿ American Geophysical Union 1994