Some diff. eq. systems constructed naively from IBP contain spurious poles that look like A/(x^n+eps*x^m+...). Being spurious these disappear in the epsilon form, but dragging them though calculations mean that the polynomial needs to be solved, which introduces roots of high powers into intermediate stages, or even worse -- no closed form solution can be found, which stops the reduction altogether (Maple signals this with the RootOf function; Maxima usually just hangs).
We need to find a way to transform these poles out without the need for solving the polynomial. This can be done by a Moser-like transformation that balances these poles with e.g. infinity. The infinity pole then can be dealt with in the usual way.
See issues #12 and #16 for comments from people who've run into about this situation.
Some diff. eq. systems constructed naively from IBP contain spurious poles that look like
A/(x^n+eps*x^m+...). Being spurious these disappear in the epsilon form, but dragging them though calculations mean that the polynomial needs to be solved, which introduces roots of high powers into intermediate stages, or even worse -- no closed form solution can be found, which stops the reduction altogether (Maple signals this with theRootOffunction; Maxima usually just hangs).We need to find a way to transform these poles out without the need for solving the polynomial. This can be done by a Moser-like transformation that balances these poles with e.g. infinity. The infinity pole then can be dealt with in the usual way.
See issues #12 and #16 for comments from people who've run into about this situation.