S U M M A R YSymmetric tensors are typically encountered during investigations associated with stress and strain analysis and, thus, they are of particular interest to geophysicists and geodesists. Furthermore, symmetric tensors are studied using eigentheory analysis which provides the decomposition of the tensor on its principal components (n independent eigenvalues and the corresponding eigenvectors). In this paper, an analytical expression of the covariance matrix of the eigenvalues and eigenvectors of an n-D symmetric tensor is derived based on the principles of linear algebra and differential calculus. Through numerical tests, the proposed formulation is proven to give realistic uncertainty estimates of the determined eigenparameters. The methodology also reveals the significant impact on uncertainty assessments when the parameter dependencies between principal components are neglected.