Efficiently and accurately simulating multi-dimensional M-coupled nonlinear Schrödinger equations with fourth-order time integrators and Fourier spectral method.