In this paper, the variational iteration method (VIM) and the Adomian decomposition method (ADM) are applied to solve numerically the focusing Manakov systems of coupled nonlinear Schrödinger equations. The accuracy of the methods are verified by ensuring that the conserved quantities remain almost constant. The results show that VIM is much easier, more convenient, more stable and efficient than ADM.