Hi Philip, thanks for this very clear article & code. I used your code to simulate the 8 figure solution of three points with mass 1/3 and found it is stable for only 1 period, then it diverges into chaos. Is this normal (due to physics?) or is this due to numerical issues in plain python/numpy? (I used the initial conditions from Eq8 from https://arxiv.org/pdf/1707.06462.pdf at the bottom of page 2. thanks, Peter.