default(realprecision,38);
T=x^3-6*x^2+11*x-6;
polrootsreal(T)
polrootsreal(T, [-oo,3/2])
polrootsreal(T, [3/2,6])
polrootsreal(T, [-oo,+oo])
polrootsreal(T, [2,3])
polrootsreal(T, [1,2])
polsturm(T, [-oo,3/2])
polsturm(T, [3/2,6])
polsturm(T, [-oo,+oo])
polsturm(T, [2,3])
polsturm(T, [1,2])
polsturm(T, [1,+oo])
polsturm(T, 2,3)
polsturm(T, 2.,3)
polsturm(T,,2)
polrootsreal(x^10 + 23*x^9 + 19*x^8 + 18*x^7 + 39*x^6 + 41*x^5 + 46*x^4 + 24*x^3 - 4*x^2 + 2*x + 42)
polrootsreal(polchebyshev(9))
polrootsreal(polchebyshev(10))
polrootsreal(x^0)
polrootsreal(1)
polrootsreal(0)
polrootsreal(Pol(0))
polrootsreal(Mod(1,2))

polroots(T*x+0.)
polroots(1)
polrootsreal(T,[1,1])
polrootsreal(T,[0,0])
polsturm(T,[1,1])
polsturm(T,[2,1])

U=(x^2-1)*(x-2);
polsturm(U)
polsturm(U,[-oo,1])
polsturm(U,[-1,+oo])

polrootsreal(x^3-2)
polrootsreal(x^3+2)
\\#1605
polsturm(33*x^2-4*x-1)
polrootsreal(4*x)
polsturm(-4*x)
polsturm((x^4-2)^2)
