r - 是否可以在 R 中求解代数方程？

Question

我想找到以下解决方案:

-x^3+6*x^2+51*x+44=0

但是用 R. 有可能吗？

我找到了 Ryacas 包，但似乎没有人能够使它工作。

可能听起来微不足道，但我无法找到一种简单的方法来做到这一点......

你有替代方案吗？

谢谢你们!

Answer 1

您可以使用 polynom包裹:

library(polynom)
p <- polynomial(c(44,51,6,-1))
# 44 + 51*x + 6*x^2 - x^3 
solve(p)
# [1] -4 -1 11

但是你可以简单地使用函数 polyroot来自 base包裹:

polyroot(c(44,51,6,-1))
# [1] -1+0i -4+0i 11+0i

如果你用 Re 保留实部:

Re(polyroot(c(44,51,6,-1)))
# [1] -1 -4 11