求助关于ode将每一步结果导出

guai乖Maggie

在用ode45进行计算的时候，可不可以将每一步的结果导出到另外一个矩阵中？
不是全部计算完了之后再导出，而是计算一步导出一个值，将这个值放到另外一个矩阵里面
就拿一个简单的常微分方程作例子：
function dy = rigid(t,y)
dy = zeros(3,1);    % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
怎么讲每一步计算出的y导入到一个矩阵X中呢？

linlin820

仔细修改ode45.m文件好像能实现