用matlab写一个二阶TVD型Runge-Kutta算法

好的，下面给出一个Matlab代码实现二阶TVD型Runge-Kutta算法：

function [t, y] = rk2_tvd(fun, tspan, y0, h)
% fun: 要求解的ODE函数，形式为 dy/dt = fun(t, y)
% tspan: 求解时间区间，形式为 [t0, tf]
% y0: 初始值
% h: 步长
t0 = tspan(1); tf = tspan(2);
t = t0:h:tf;
y = zeros(length(y0), length(t));
y(:,1) = y0;
for n = 1:length(t)-1
    tn = t(n);
    yn = y(:,n);
    k1 = h * feval(fun, tn, yn);
    k2 = h * feval(fun, tn + h/2, yn + k1/2);
    y1 = yn + k2;
    k3 = h * feval(fun, tn + h/2, (3*yn + y1)/4);
    y(:,n+1) = (yn + y1 + 2*k3)/3;
end
end

其中，输入参数fun是要求解的ODE函数，形式为 dy/dt = fun(t, y)，tspan是求解时间区间，形式为 [t0, tf]，y0是初始值，h是步长。输出参数t是时间序列，y是对应的函数值序列。在这个函数中，我们用了二阶TVD型Runge-Kutta方法来求解ODE。该方法与传统的二阶Runge-Kutta方法类似，但在计算k3时，使用了TVD限制器，使得算法具有更好的稳定性和精度。