Функция для полиномов произвольного порядка (предпочтительно символический метод)
Я нашел полиномиальные коэффициенты из моих данных:
R <- c(0.256,0.512,0.768,1.024,1.28,1.437,1.594,1.72,1.846,1.972,2.098,2.4029)
Ic <- c(1.78,1.71,1.57,1.44,1.25,1.02,0.87,0.68,0.54,0.38,0.26,0.17)
NN <- 3
ft <- lm(Ic ~ poly(R, NN, raw = TRUE))
pc <- coef(ft)
Так что я могу создать полиномиальную функцию:
f1 <- function(x) pc[1] + pc[2] * x + pc[3] * x ^ 2 + pc[4] * x ^ 3
И, например, взять производную:
g1 <- Deriv(f1)
Как создать универсальную функцию, чтобы ее не нужно было переписывать для каждой новой степени полиномаNN
?