Summary: A mixed finite element discretization is applied to Richard’s equation, a nonlinear, possibly degenerate parabolic partial differential equation modeling water flow through porous medium. The equation is considered in its pressure formulation and includes both variably and fully saturated flow regime. Characteristic of such problems is the lack in regularity of the solution. To handle this, we use a time-integrated scheme. We analyze the scheme and present error estimates showing its convergence.