I agree in general, and it would have been nice to see some plots of some of quantities that should be conserved, e.g. how well does his simulation preserve the initial energy and angular momentum?
Having said that, it is an adaptive RK scheme, and it seems to work pretty well, the article shows that the results match some reference data quite well, and it even captures some quite subtle effects, most notably the influence of the moon.