一般多项式的R-Picard方法

我目前正在用R写一个程序，用Picard的方法求一般多项式差分方程的解。 要了解其背后的数学知识（因为这里没有数学模式）：

从那以后，我一直在尝试使用Ryacas软件包进行集成。然而，我在尝试使用表达式和集成函数的组合时遇到了麻烦

library(Ryacas)

degrees = 3
a = c(3,5,4,6)
x0 = -1
maxIterations(10)


iteration = vector('expression', length = maxIterations)
iteration[1] = x0

for(i in 2:maxIterations){

  for(i in 1:degrees){

  exp1 = expression( a[i] * iteration[i-1] ^ i) 

  }

iteration[i] = x0 + Integrate(exp1, t)
}

但这导致了

粘贴时出错（（“，…，”）： 无法将类型“closure”强制为“character”类型的向量

和

exp1=表达式（一个[j]*迭代[i-1]^j）

而不是我试图实现的实际表达式。我是否可以确保R将其作为一个真正的表达式来读取（例如，I=2的

3*（x0）^j

）

提前谢谢

编辑： 我还找到了

Subst（）

函数，目前正在尝试使用它查看是否有任何问题可以解决。现在，我主要在努力为a的m个系数建立一个表达式，因为我找不到一种方法来创建例如

expression（）

命令中的for循环