Solving the KdV hierarchy with self-consistent sources by inverse scattering method

The evolution of the eigenfunctions in the Lax representation of the KdV hierarchy with self-consistent sources possesses singularity. By proposing a method to treat the singularity to determine the evolution of scattering data, the KdV hierarchy with self-consistent sources is integrated by the inverse scattering method. The soliton solutions of these equations are obtained. It is shown that the insertion of a source may cause the variation of the speed of soliton. This approach can be applied to other (1+1)-dimensional soliton hierarchies.