This note documents R code for the CANDECOMP and TUCKER generalizations of the Singular Value Decomposition to multiway arrays. Alternating Least Squares algorithms are used to generate a convergent sequence of low-rank approximations.