
D.12.8.11 sturmseq
Procedure from library rootsur.lib (see rootsur_lib).
 Usage:
 sturmseq(p); p poly
 Return:
 list: a Sturm sequence of p
 Assume:
 p is a univariate polynomial with rational coefficients
 Theory:
 The Sturm sequence of p (also called remainder sequence) is the
sequence beginning with p, p' and goes on with the negative part of
the remainder of the two previous polynomials, until the remainder
is zero.
See: Basu, Pollack, Roy, Algorithms in Real Algebraic Geometry,
Springer, 2003.
Example:
 LIB "rootsur.lib";
ring r = 0,(z,x),dp;
poly p = x53x4+12x3+7x153;
sturmseq(p);
==> [1]:
==> x53x4+12x3+7x153
==> [2]:
==> x412/5x3+36/5x2+7/5
==> [3]:
==> x39/7x25/3x+317/7
==> [4]:
==> x2756/151x+2433/151
==> [5]:
==> x514191/177889
==> [6]:
==> 1
 See also:
sturm;
sturmhaseq.
