如何跟踪 ode 求解器的显式状态？

问：

在我的导数函数中，我有一个导数列表，这些导数依赖于当前时间步长的值。但是，我也有一些中间状态对这个导数有贡献。例如：这是我意思的一个简单的例子。我还想跟踪每个时间点的这个状态变量的中间值。这怎么能做到？这些状态没有导数。现在我的（伪代码）模式如下：dy/dtSTATE(i)=2-RATE(i);

partition t manually
for 1:length(t)-1
  compute states(i)
  rates(i+1) = ode(dydt with updated states(i), [t(i), t(i+1)], rates(i))
end
plot rates
plot states 

通常，你会把所有状态都放在导数函数中。但我也想跟踪它以进行绘图+分析。