
D.12.8.13 sturmhaseq
Procedure from library rootsur.lib (see rootsur_lib).
 Usage:
 sturmhaseq(P); P poly.
 Return:
 list: the nonzero polynomials of the SturmHabicht sequence of P
 Assume:
 P is a univariate polynomial.
 Theory:
 The SturmHabicht sequence (also subresultant sequence) is closely
related to the Sturm sequence, but behaves better with respect to
the size of the coefficients. It is defined via subresultants.
See: Basu, Pollack, Roy, Algorithms in Real Algebraic Geometry,
Springer, 2003.
Example:
 LIB "rootsur.lib";
ring r = 0,x,dp;
poly p = x5x4+x3/2;
list l = sturmhaseq(p);
l;
==> [1]:
==> 132949/16
==> [2]:
==> 25x332
==> [3]:
==> 16x2+42x24
==> [4]:
==> 4x320x+73/2
==> [5]:
==> 5x44x3+1
==> [6]:
==> x5x4+x3/2
 See also:
sturm;
sturmha;
sturmseq.
