如何存储ode求解器的中间变量？

Question

来自ODE求解器的部分代码是：

function rdot = f(t, r)

rdot(1) = r(2);
calculation expressions
...
...
rdot(2) =  some expression;

rdot(3) = r(4);

calculation expressions;
...
...
rdot(4) =  some expression;
rdot = rdot';

我希望为每个t的输入值存储rdot（1）和rdot（4）的值。稍后，我需要它们用于其他两个变量f1和f2的其他计算。所以我想我需要对两个新变量进行数组存储rdot（3）和rdot（4）。但是我该如何实现呢？

或者我应该做

function rdot = f(t, r, f1, f2)

在此功能代码的同一部分中对f1和f2的计算吗？但这可能会影响ode求解器？