  In this paper, we study the numerical solution of Manakov systems by using a
spectrally accurate Fourier decomposition in space, coupled with a spectrally
accurate time integration. This latter relies on the use of spectral
Hamiltonian boundary Value Methods. The used approach allows to conserve all
the physical invariants of the systems. Some numerical tests are reported.
