在Julia

Question

我有一堆耦合的非线性ODE，其系统可以总结如下：

其中F，G，H和I是使用所有变量和某些参数p的非线性函数，这些参数也取决于某些变量值。

我可以使用for或while循环或一些自制函数来解决它，但是我想当我想在Julia方面变得更好时，我会尝试使用DifferentialEquation.jl，我认为它将允许更快或更有效的计算，而且我现在已经阅读了不同的示例以及文档的一部分。

我遇到的问题如下：

• 所有变量的大小都不同，我已经看到了计算向量演化的示例（我可以将它们全部写成一个大向量，但是检索矩阵会很麻烦，而且这意味着将所有内容存储在每个时间步），但大小变量没有不同。
• 我将需要在一个循环之后存储所有这些变量（可以通过使用回调定义该结束），但是我只需要在每个时间步存储U [1]和U [2]。

我还没有找到任何类似于我的例子。有人可以帮忙吗？

提前谢谢！