this thesis discusses the triangular finite element solution to second order elliptic boundary value problems. the barycentric coordinate system, which some engineers call the areal coordinate system, is used throughout in this thesis. some fundamental parts of vector calculus are developed in this coordinate system, and are applied to the triangular finite element method.