解決 R 中的 ODE

請注意，必須以與狀態變數規範相同的順序返回變化率。例如洛倫茲模型這意味著，在函式 Lorenz 命令中

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

與狀態變數的定義具有相同的順序

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