我有一个方程式,我想解决的价值 m
m^2+x*m+1=0
Run Code Online (Sandbox Code Playgroud)
我知道我可以找到这些根源polyroot(c(1,x,1)).
我想在R中绘制这些根的变化x,以便为每个值绘制多项式的所有实根x.
注意,多项式可能变得比这更复杂,排除了二次或三次公式的应用.
有关如何实现这一目标的任何想法?
谢谢!
我认为这样做会有所帮助.只需设置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)
Run Code Online (Sandbox Code Playgroud)

| 归档时间: |
|
| 查看次数: |
780 次 |
| 最近记录: |