有没有办法在MATLAB中输出ode求解器的中间计算？

% 定义输出函数
function status = outputFcn(t, y, flag)
    persistent intermediate_results;

    % 在第一个时间步骤上初始化中间结果
    if isempty(intermediate_results)
        intermediate_results = [t, y];
    else
        intermediate_results = [intermediate_results; [t, y]];
    end

    % 在某些特定时间步骤上进行处理
    if strcmp(flag, 'init')
        % 处理初始时间步骤
        disp('Init step');
    elseif strcmp(flag, 'done')
        % 处理最后一个时间步骤
        disp('Final step');
    else
        % 处理中间时间步骤
        disp('Intermediate step');
    end

    status = 0; % 继续求解器的运行
end

% 定义ODE函数
function dydt = odefun(t, y)
    dydt = -2 * t * y;
end

% 设定求解器选项和参数
options = odeset('OutputFcn', @outputFcn);
tspan = [0 1];
y0 = 1;

% 调用 ode 求解器
[t, y] = ode45(@odefun, tspan, y0, options);