basic_astroquery.py

from astroquery.jplhorizons import Horizons

pluto        = 999
jupiter      = 599
mu69         = 486958
new_horizons = -98

paths = []
for body_name in [jupiter, pluto, new_horizons, mu69]:
    obj = \
        Horizons(
            id=body_name,
            id_type="majorbody",
            epochs={
                'start':'2006-1-20',
                'stop':'2018-12-28',
                'step':'10d'
            }
        )
    paths.append([[i["x"], i["y"], i["z"]] for i in obj.vectors()])

print paths