  A symbolic computation technique is developed to calculate adiabatic
evolution equations for parameters of the perturbed DNLS/MNLS solitons obtained
by the recently developed direct perturbation theory [X.-J. Chen and J. Yang,
Phys. Rev. E {\bf 65}, 066608(2002)]. Effects of the intrapulse Raman
scattering, third-order group velocity dispersion, and narrow-banded filters
with nonlinear gain on the MNLS solitons are studied as examples.
