r - 在R中找到有理函数的局部最大值

Question

如何在 R 中找到诸如 f(x) = 500*t/(2*t^2 + 9) 之类的有理(多项式)函数的局部最大值？我尝试了以下方法:

t <- seq(0,10000,by=0.000001)
f <- 500*t/(2*t^2 + 9)
max(f)

但是我很快就耗尽了内存。

Answer 1

所以，这就是我的做法。感谢 G. 和 alistaire。

f <- function(x){
  z <- 500*x/(2*x^2+9)
  return(z)
}

optimize(f,c(0,1000),maximum=TRUE)

Output:
$`maximum`
[1] 2.121332

$objective
[1] 58.92557