I recently wrote a Sage / Python program that uses
JPL Horizons data to create interactive 3D plots of trajectories of Solar System bodies. You can specify any target & observation centre that Horizons knows about: (1,169,191 asteroids, 3,778 comets, 211 planetary satellites {includes satellites of Earth and dwarf planet Pluto}, 8 planets, the Sun, L1, L2, select spacecraft, and system barycenters).