The aim of this work is the numerical study of a nonlinear equation, which models the water flow in a partially saturated underground porous medium under the surface. We propose a discretization of this equation that combines Euler’s implicit scheme in time and spectral methods in space. We prove optimal error estimates between the continuous and discrete solutions. Some numerical experiments confirm the interest of this approach. We present numerical experiments which are in perfect coherence with the analysis.