I figured it out myself. In file rootisolation.py of Sympy, only the square-free part of the polynomial is considered for the calculation of the Sturm chain. The function dup_sqf_part(f, K) is the first thing called.
1 Like