Convergence of a finite element method based on the dual variational formulation

An "equilibrium model" with piecewise linear polynomials on triangular clements applied to the solution of a mixed boundary value problem for a second order elliptic equation is studied. The procedure is proved to be second order correct in $h$ (the maximal side in the triangulation) provided the exact solution is sufficiently smooth.