A stable, direct solver for the gradient equation

A new finite element discretization of the equation is introduced. This discretization gives rise to an invertible system that can be solved directly, requiring a number of operations proportional to the number of unknowns. We prove an optimal error estimate, and furthermore show that the method is stable with respect to perturbations of the right-hand side . We discuss a number of applications related to the Stokes equations.