我正在使用rootSolve包解决目标函数的根。这是代码,它给我一个显示根的图表并打印解决方案:
library(rootSolve)
P_func <- function (x) (0.11)*(-10^6/(10^4-(x/0.03))+20)*(10^4-(x/0.03)^2/(10^6))
curve(P_func(x), 0, 8000)
abline(h = 0, lty = 3)
All <- uniroot.all(P_func, c(0,8000))
points(All, y = rep(0,length(All)), pch = 16, cex = 2)
All
我想抓住根并将其插入另一个功能:
W_func <- function (x) (10^4-x/0.03)
x的值由我为P_func找到的每个解决方案定义。有一种简单的方法可以做到这一点吗?
答案 0 :(得分:0)
是的确有一种简单的方法。
All
是一个向量,您的W_func
接受一个向量作为参数。所以只需使用W_func(All)
。