Согласно PARI/GP, в соответствующем расширении многочлен разваливается на квадратные множители, корни которых находятся без проблем по стандартной формуле.
Код:
? nffactor( nfinit( y^8+y^7-7*y^6-6*y^5+15*y^4+10*y^3-10*y^2-4*y+1), (x^8-x^7-7*x^6+6*x^5+43/4*x^4-23/4*x^3-3/2*x^2+4*x-13/4)*4 )
%1 =
[x^2 + Mod(y^4 - 4*y^2 + y + 2, y^8 + y^7 - 7*y^6 - 6*y^5 + 15*y^4 + 10*y^3 - 10*y^2 - 4*y + 1)*x + Mod(-3/2*y^7 - 3/2*y^6 + 10*y^5 + 8*y^4 - 19*y^3 - 19/2*y^2 + 17/2*y, y^8 + y^7 - 7*y^6 - 6*y^5 + 15*y^4 + 10*y^3 - 10*y^2 - 4*y + 1) 1]
[x^2 + Mod(y^5 - 4*y^3 + 2*y, y^8 + y^7 - 7*y^6 - 6*y^5 + 15*y^4 + 10*y^3 - 10*y^2 - 4*y + 1)*x + Mod(1/2*y^7 + 1/2*y^6 - 4*y^5 - 4*y^4 + 9*y^3 + 17/2*y^2 - 11/2*y - 7/2, y^8 + y^7 - 7*y^6 - 6*y^5 + 15*y^4 + 10*y^3 - 10*y^2 - 4*y + 1) 1]
[x^2 + Mod(-y^7 - y^6 + 6*y^5 + 5*y^4 - 10*y^3 - 5*y^2 + 4*y - 1, y^8 + y^7 - 7*y^6 - 6*y^5 + 15*y^4 + 10*y^3 - 10*y^2 - 4*y + 1)*x + Mod(1/2*y^7 + 1/2*y^6 - 4*y^5 - 2*y^4 + 9*y^3 + 1/2*y^2 - 7/2*y + 1, y^8 + y^7 - 7*y^6 - 6*y^5 + 15*y^4 + 10*y^3 - 10*y^2 - 4*y + 1) 1]
[x^2 + Mod(y^7 + y^6 - 7*y^5 - 6*y^4 + 14*y^3 + 9*y^2 - 7*y - 2, y^8 + y^7 - 7*y^6 - 6*y^5 + 15*y^4 + 10*y^3 - 10*y^2 - 4*y + 1)*x + Mod(1/2*y^7 + 1/2*y^6 - 2*y^5 - 2*y^4 + y^3 + 1/2*y^2 + 1/2*y + 3/2, y^8 + y^7 - 7*y^6 - 6*y^5 + 15*y^4 + 10*y^3 - 10*y^2 - 4*y + 1) 1]
P.S. Ну или можно напрямую решить уравнение в косинусах, как
мы это делали раньше.