They describe how Richardson extrapolation can be used to improve the accuracy of the solutions from 𝒪 (ɛ) to 𝒪(ɛ2), but it appears that this has not been done for their published results.