The paper aims to prove the orbital stability of peakon solutions in the integrable Degasperis-Procesi (DP) equation,
which resembles another integrable model, Camassa-Holm (CS) equation,
although being essentially different from it. Nevertheless, the simplest peakon solutions are identical in both equations,
with arbitrary positive constant . Unlike peakons, smooth solutions to these equations feature an infinite propagation velocity. The DP equation also has exact solutions in the form of “shock peakons”,
which can be generated by a collision of a symmetric pair of peakons, set initially at a finite distance. The stability of peakons within the framework of the CS equation was established previously. The difference in the proof of the peakon stability in the DP equation, which is the main result of this paper, is related to the different form of its dynamical invariants. In particular, the second invariant, , where is defined by relation , includes the nonlocal density. The final result is that a slightly perturbed peakon evolves into a solution close to another peakon, which may have a different velocity.