Generalized perturbation equations of celestial mechanics

Generalized perturbation equations of celestial mechanics in terms of orbital elements are derived. The most general case is considered: Keplerian motion of two bodies caused by gravitational forces between them is disturbed by disturbing acceleration acting on each of the bodies separately and by changes of masses of these bodies. It is also pointed out why derivation presented in Klaka (1992a) is completely physically correct only for constant masses.