  The current implementation of Thiele rational interpolation in Maple (the
ThieleInterpolation routine) breaks down when the points are not well-ordered.
In this article, it is shown how this breakdown can be avoided by ordering the
interpolation points in an adaptive way.
