Veamos Secuencias de Sturm . Aquí es más fácil utilizar un sistema de álgebra computacional, daré programas de Maxima al final.
Comenzamos con $P_0=P=x^8-x^7+x^2-x+15$ y $P_1=P'$ , entonces hacemos la división polinómica hasta que demos con un polinomio constante. Obtenemos: $$P_0=x^8-x^7+x^2-x+15$$ $$P_1=8x^7-7x^6+2x-1$$ $$P_2=\frac{7}{64}x^6-\frac{3}{4}x^2+\frac{27}{32}x-\frac{959}{64}$$ $$P_3=-\frac{384}{7}x^3+\frac{768}{7}x^2-1152x+960$$ $$P_4=-\frac{1877}{64}x^2-\frac{6571}{64}x+\frac{32501}{256}$$ $$P_5=\frac{8617090464}{3523129}x-\frac{7984345440}{3523129}$$ $$P_6=-\frac{13706250843331201741}{2062618001798881536}$$
Ahora el número de raíces entre $a$ y $b$ es $\sigma(a)-\sigma(b)$ , donde $\sigma(x)$ es el número de cambios de signo en $P_0(x), P_1(x),..., P_6(x)$ .
Necesitamos un límite superior para las raíces (o en su lugar calcular $\sigma(-\infty)-\sigma(\infty)$ ). Utilicemos El límite de Cauchy que está aquí $16$ . Entonces $\sigma(-16) = \sigma(16) = 3$ por lo que la diferencia es cero y no hay raíces reales en $[-16, 16]$ por lo que no hay una raíz real en $\mathbb{R}$ .
Aquí hay un código de Maxima para calcular todo esto.
sturmseq(p) := block(\[u:\[p,diff(p,x)\],q,r,n:2\],
while hipow(u\[n\],x)>0 do (
\[q,r\]:divide(u\[n-1\],u\[n\]),
u:endcons(expand(-r),u),
n:n+1),
u)$
changes(u,a) := block(\[s:subst(x=a,u\[1\]),n:length(u),i,j:0,t\],
for i from 2 thru n do (
t:subst(x=a,u\[i\]),
if (t>0 and s<0) or (t<0 and s>0) then j:j+1,
s:t),
j)$
p:x^8-x^7+x^2-x+15$
s:sturmseq(p);
changes(s,-16);
changes(s,16);