Institute of Applied Astronomy, 197042 St.Petersburg, Russia
Using an original numerical technique we construct intermediate orbit for general planetary theory in the form of multivariate Fourier series with numerical coefficients, the arguments of the Fourier series being linear functions of time. We illustrate the structure and efficiency of the derived series giving various statistical properties of the coefficients.
We investigate the ability of the recently proposed elliptic function approach to compactify the Fourier series representing the intermediate orbit. Our results confirm that for univariate Fourier series the elliptic function approach is quite efficient and allows one to compactify the series substantially. However, for multivariate Fourier series the elliptic function approach does not give any advantages.