The eigenvalue problem for 3x3 octonionic Hermitian matrices contains somesurprises, which we have reported elsewhere. In particular, the eigenvaluesneed not be real, there are 6 rather than 3 real eigenvalues, and thecorresponding eigenvectors are not orthogonal in the usual sense. Thenonassociativity of the octonions makes computations tricky, and all of theseresults were first obtained via brute force (but exact) Mathematicacomputations. Some of them, such as the computation of real eigenvalues, havesubsequently been implemented more elegantly; others have not. We describe herethe use of Mathematica in analyzing this problem, and in particular its use inproving a generalized orthogonality property for which no other proof is known.