  We propose two finite-difference algorithms of fourth order of accuracy for
solving the initial problem of the Zakharov-Shabat system. Both schemes have
the exponential form and conserve quadratic invariant of Zakharov-Shabat
system. The second scheme contains the spectral parameter in exponent only and
allows to apply the fast computational algorithm.
