  We propose a new high-precision algorithm for solving the initial problem for
the Zakharov-Shabat system. This method has the fourth order of accuracy and is
a generalization of the second order Boffetta-Osborne scheme. It is allowed by
our method to solve more effectively the Zakharov-Shabat spectral problem for
continuous and discrete spectra.
