Revision history [back]
 | 1 | initial version |
We tried a bit more and it turns out that dividing by a small irreducible polynomial will actually work:
sturm(((2*k+1) *a + 1)/(k^2+1))
=> [((2*k + 1)/(k^2 + 1))*a + 1/(k^2 + 1), (2*k + 1)/(k^2 + 1)]
The spurious factor cannot change the number of sign changes, so that's a workaround.
