A finite element method is developed for approximating the solution of the Dirichlet problem for the biharmonic operator, as a canonical example of a higher order elliptic boundary value problem. The solution is approximated by special choices of classes of discontinuous functions, piecewise polynomial functions, by virtue of a special variational formulation of the boundary value problem. The ...