如果我需要计算两个变量函数的泰勒级数近似值，那么在R中是否可以？我感兴趣的函数被提升到变量x1和x2之和的幂：e ^（x1 + x2）
由于taylor库中有pracma function，我尝试这样的事情：
func <- function(x) exp(1)^(x[1]+x[2])
vars <- c(0, 0) # evaluated at point (0,0)
taylor(func, vars, 2)
但R说Error in fderiv(f, x0, i) : Function 'f' must first be vectorized: Vectorize(f).
那么，我做错了什么，并且完全可以用于两个变量的函数？