2

の値について解きたい方程式がありますm

m^2+x*m+1=0

私はこれらのルーツを使用して見つけることができることを知っていますpolyroot(c(1,x,1))。

x多項式のすべての実根がの値ごとにプロットされるように、これらの根をRにプロットして変化させたいと思いますx。

多項式はこれよりも複雑になり、2次式または3次式の適用が除外される可能性があることに注意してください。

私がこれをどのように達成できるかについて何か考えはありますか？

ありがとう！

4

1 に答える 1

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)

ここに画像の説明を入力してください

于 2012-08-07T16:03:09.210 に答える