我的R for循环总结了一个系列有什么问题？

Question

这是我的函数循环:

answer = function(a,n) {
  for (k in 0:n) {
    x =+ (a^k)/factorial(k)
    }
  return(x)
  }

answer(1,2) 应该返回2.5,因为它是计算值

1 ^ 0/0!+ 1 ^ 1/1!+ 1 ^ 2/2!= 1 + 1 + 0.5 = 2.5

但我明白了

answer(1,2)
#[1] 0.5

看起来它无法累积所有三个术语,每次只存储最新值.+=不起作用所以我用过,=+但它仍然不对.谢谢.

Answer 1

answer = function(a,n) {
  x <- 0  ## initialize the accumulator
  for (k in 0:n) {
    x <- x + (a^k)/factorial(k)  ## note how to accumulate value in R
    }
    return(x)
  }

answer(1, 2)
#[1] 2.5

有"矢量化"解决方案:

answer = function(a,n) {
  x <- a ^ (0:n) / factorial(0:n)
  return(sum(x))
  }

在这种情况下,您不需要初始化任何内容.R将分配背后存储器<-和sum.

您正在使用泰勒展开来近似exp(a).请参阅此主题的问答.您可能需要特别注意我的答案中提到的"数值收敛"问题.