我有一个方程式,我想为m
m^2+x*m+1=0
我知道我可以使用polyroot(c(1,x,1))
找到这些根。
我想在R中将这些根绘制成不同的x
,以便为x
的每个值绘制多项式的所有实根。
注意,多项式可能比这更复杂,排除了二次或三次公式的应用。
有关如何实现这一目标的任何想法?
谢谢!
答案 0 :(得分:8)
我认为这样做会有所帮助。只需将X
设置为您要评估的任何值。
X <- seq(0, 10, length=21)
roots <- sapply(X, function(x) polyroot(c(1,x,1)))
roots[abs(Im(roots)) > 1e-10] <- NA
matplot(X, t(roots), pch=1)