The development of numerical tools, including modern symplectic methods, are presented as they pertain to the identification of short-term chaos and long term integrations of the orbits of many astronomical architectures such as stellar triples, planets in binaries, and single stars that host multiple exoplanets.