A mathematical method for studying the effects of small perturbations on the motion of celestial objects.