descr = 'obtaining positions of the Sun and planets in au' sun = astropy.coordinates.get_body_barycentric ('sun', t) mercury = astropy.coordinates.get_body ...