Digital series are applied to the solution of linear elastostatic problems using a mixed finite element formulation. The hierarchical nature of the digital series suggests the use of adaptive, macro-elements. The properties of the digital functions are exploited to establish closed-form solutions for the integral expressions of the structural matrices. Implementation of numerical integration procedures can be thus avoided. The structure of the governing system is well suited to the use of efficient algorithms for the solution of large, sparse systems of linear equations. Determination and graphic representation of the stress and displacement fields are easily implementable because both fields are directly interpolated and manipulation of linear combinations of digital functions are extremely simple to perform.