Estimate of the Lyapunov timescale as a function of the number of particles. Here the horizontal axis is not linear, but in ln(ln(N)) to illustrate the scaling proposed in Goodman et al. (1993). The different symbols and colors represent different calculations (see legend). The vertical bars, plotted for Newton’s Hermite only, show the root-mean-square of the dispersion in the series of solutions. The error bars in the results obtained with Brutus are statistically indistinguishable from the presented bars.