On several problems in the numerical integration of celestial orbits