Rotation matrices that minimize or maximize the sum of the squared distances between corresponding atoms for two structures are found using a constrained least-squares procedure solved analytically as an eigenvalue problem in quaternion parameters.