Hey there!  I&#x27;m a grad student in my last year and my research is in orbital dynamics.  I&#x27;ve recently started getting into the rebound package written by Hanno Rein:<p><a href="https:&#x2F;&#x2F;github.com&#x2F;hannorein&#x2F;rebound" rel="nofollow">https:&#x2F;&#x2F;github.com&#x2F;hannorein&#x2F;rebound</a><p>There&#x27;s some tutorials that come with it.  It can also read data from JPL&#x27;s ephemerides to calculate orbits of various planets and minor bodies in the Solar System.  It has a couple of different integrators (some symplectic, some not), but one of them (IAS15) is new and (as far as I know) the best for long-term integrations of Solar-System-like systems.<p>Also if you&#x27;re interested in studying the long-term evolution of hierarchical triples, I wrote a Python package to study that:<p><a href="https:&#x2F;&#x2F;github.com&#x2F;joe-antognini&#x2F;kozai" rel="nofollow">https:&#x2F;&#x2F;github.com&#x2F;joe-antognini&#x2F;kozai</a><p>As far as star charts, you can do things like that with AstroPy.  If you want the coordinates of stars, look into astropy.vo, and then if you want to plot them, you can do that with some of the functions in astropy.wcs.