Solution (source code)

= Solution

Euler's equations are $I_1\dot\omega_1=(I_2-I_3)\omega_2\omega_3$ and cyclic permutations. Setting two components to zero gives uniform rotation about any principal axis.