  Employing a particularly suitable higher order symplectic integration
algorithm, we integrate the 1-$d$ nonlinear Schr\"odinger equation numerically
for solitons moving in external potentials. In particular, we study the
scattering off an interface separating two regions of constant potential. We
find that the soliton can break up into two solitons, eventually accompanied by
radiation of non-solitary waves. Reflection coefficients and inelasticities are
computed as functions of the height of the potential step and of its steepness.
