Fourier expansions of elliptic motion functions in multiples of the true, eccentric, elliptic and mean anomalies are computed numerically by means of the fast Fourier transform. Both Hansen-like coefficients and their derivatives with respect to eccentricity of the orbit are considered. General behavior of the coefficients and the efficiency (compactness) of the expansions are investigated for various values of eccentricity of the orbit.