We already have the possibility to get positions of celestial bodies with very high precision via JPL ephemerides but you would have to write and integrate your force model yourself.