A theory of mixed finite element approximations of a class of linear self-adjoint boundary value problems is developed which involves splitting a problem of the type T*Tu equals f into a pair of canonical equations, Tu equals v, T*v equals f. Consistency and stability of Galerkin/finite element approximations of the pair is described, convergence criteria are established, and error estimates are derived for a number of important cases.
