Я решаю для корней целевой функции, используя пакет 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. Есть ли простой способ сделать это?