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