Я не знаю ни одного активного проекта R / octave, но если вы только что нашли корни для данного полинома, вы можете использовать один из пакетов polynom или PolynomF :
Вот пример с P (x) = 6 + 5 * x + 4 * x ^ 2 + 3 * x ^ 3 + 2 * x ^ 4 + x ^ 5.
в октаве,
octave[2] > p = 1:6;
octave[3] > roots(p)
ans =
0.55169 + 1.25335i
0.55169 - 1.25335i
-1.49180 + 0.00000i
-0.80579 + 1.22290i
-0.80579 - 1.22290i
В R,
> library(polynom)
> p <- polynomial(6:1)
> pz <- solve(p)
> pz
[1] -1.491798+0.000000i -0.805786-1.222905i -0.805786+1.222905i
[4] 0.551685-1.253349i 0.551685+1.253349i