High-order IMEX-spectral schemes for computing the dynamics of systems of nonlinear Schrödinger/Gross-Pitaevskii equations