解决 R 中的 ODE

请注意，必须以与状态变量规范相同的顺序返回变化率。例如洛伦兹模型这意味着，在函数 Lorenz 命令中

return(list(c(dX, dY, dZ)))

与状态变量的定义具有相同的顺序

yini <- c(X = 1, Y = 1, Z = 1)