Home Online Manual
Top
Back: sturm
Forward: sturmha
FastBack:
FastForward:
Up: rootsur_lib
Top: Singular Manual
Contents: Table of Contents
Index: Index
About: About this document

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 = x5-3x4+12x3+7x-153;
sturmseq(p);
==> [1]:
==>    x5-3x4+12x3+7x-153
==> [2]:
==>    x4-12/5x3+36/5x2+7/5
==> [3]:
==>    -x3-9/7x2-5/3x+317/7
==> [4]:
==>    -x2-756/151x+2433/151
==> [5]:
==>    x-514191/177889
==> [6]:
==>    1
See also: sturm; sturmhaseq.