We are revisiting the problem of solving a discrete nonlinear Schrödinger equation by the inverse scattering transform method, by use of the recently developed ExactMPF package within MAPLE Software. ExactMPF allows for an exact Wiener–Hopf factorization of matrix polynomials regardless of the partial indices of the matrix. The package can be widely used in various problems, where Wiener–Hopf factorization as one of the effective mathematical tools is required, as its code has already been disclosed. The analysis presented in this paper contains not only numerical examples of its use, but is also supported by appropriate and accurate a priori estimations. The procedure itself guarantees that the ExactMPF package produces all computations arithmetically exactly, and a detailed numerical analysis of various aspects of the computational algorithm and approximation strategies is provided in the case of a finite initial impulse.