This note describes briefly a technique, easily implemented with most compu
ter algebra systems, for the purpose of computing the two-body expansions (
in powers of the eccentricity and in Fourier series of the mean-anomaly) of
a large class of functions of the distance, the true anomaly and/or the ec
centric anomaly.