PolSturm(P) = Number of real roots in real polynomial P:
P is a polynomial with real coefficients and must be square free for the function to work. An alarm is given if P is not real or if P is not square free.
P = PolMonic(P)
a = 0
for (int r = 1; r <= P.m; r++) a = a + Abs(p[r])
return PolSturmI(p, −a, a)
Return to Matrix and Polynomial Computations
Return to Harry's Home Page