A fast and efficient numerical integration algortihm is presented for the problem of the secular evolution of the spin axis by reducing the equations of motion to the Hamiltonian for with Lie-Poisson bracket.