We develop a Fourier–Galerkin spectral technique for computing the solutions of type of interacting localized waves. To this end, a special complete orthonormal system of functions in L2(−∞,∞) is used and a time-stepping algorithm implementing the spectral method is developed. The rate...