﻿功能:求多项式Sturm变号数

格式:
PolySturm(f,a,b)
PolySturm(f,a)
PolySturm(f)

f : 一元多项式符号变量
a : 指定点的Sturm变号数
b : 计算Sturm变号数区间上限

说明:
1. 本函数主要计算一元多项式变号数,方便确定实数根区间
2. PolySturm(f,a): 表示计算多项式f在a点处的Sturm变号数
3. PolySturm(f,a,b): 表示多项式f在区间[a,b]之间变号数之差
4. PolySturm(f):在实数区间内实根个数
5. 等式
$$
\begin{cases}
\mathbf{PolySturm}(f,a,b)&=|\mathbf{PolySturm}(f,a) -\mathbf{PolySturm}(f,b)|
\\
\mathbf{PolySturm}(f) &=\mathbf{PolySturm}(f,-\inf, \inf)
\end{cases}
$$

例子:
f=@x^5+3*x^7-8*x^2+4*x^3-7*x^5-123.45;
n = PolySturm(f,10)//执行得到如下结果
n =
[ 3.00000000000000 ]
