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