We applied a new numerical technique, developed recently by one of the authors and based on a method found by Zadunaisky, to a classical problem of Celestial Mechanics. The problem consists in the discovery of a massive object whose existence is assumed on the base of systematic deviatiosn form the orbits of known planets.