2013-10-10 43 views
0

我有amatrix m其包括三次多項式的係數的不同值,並且我想找到這些三次函數的最大實根。我怎樣才能做到這一點R中?應用一個矩陣的每一行用於polyroot中的R

 m = matrix(c(1, 3, 4, 6, 7, 8, 9, 2, -3, -5, 6, 2), nrow = 3) 
     colnames(m) <- c("d", "c", "b", "a") 

我想用polyroot爲m矩陣中的每一行,我怎麼能做到這一點的R'

回答

1

這是不是僅僅是:

apply(m , 1, polyroot) 

...記憶,當然,該結果將在列的順序。

apply(apply(m, 1,polyroot) , 2, function(x) ifelse(abs(Im(x)) < 0.00001, x,NA)) 
      [,1]   [,2]   [,3] 
[1,]   NA -0.416256-0i -0.4167669+0i 
[2,]   NA   NA   NA 
[3,] 2.34748+0i   NA   NA 
+0

謝謝,但我怎麼才能得到最大的真正的根? – rose

+0

嗯,這是不是原來的問題了,是嗎? –

+0

不,soory,我忘了提到它。 – rose

相關問題