This paper considers matrices arising from the use of finite element techniques in least-squares approximation and in elliptic partial differential equations; it studies their properties of numerical stability, and in particular, it establishes bounds for their inverses with respect to the uniform norm.