r - Rで3次関数の最大実根を見つける最速の方法

Question

Rで3次関数の最大根を見つける最速の方法は何ですか?

    a x^3 + b x^2 + c x + d = 0

Answer 1

基本機能に問題はありますpolyrootか？

説明

実数または複素数の多項式の零点を見つけます。

立方体の例

 polyroot(c(1,3,3,1))


 # [1] -1+0i -1+0i -1-0i

これは、多項式の最大非複素根を見つける関数です...

maxReal <- function(params){
  x <- polyroot(params)
  reals <- sapply(x, function(i) isTRUE(all.equal(Im(i),0)))
  max(Re(x)[reals])
  }