找到R中最大实数三次函数根的最快方法

Question

什么是在R中找到三次函数最大根的最快方法？

    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])
  }