Matrix techniques for various types of diagonalizations of matrices and three-dimensional arrays are implemented in R using Jacobi plane rotations.