matlab欧拉法和改进欧拉法代码

下面是使用Matlab实现欧拉法和改进欧拉法的代码示例：

欧拉法：

function [t, y] = euler(f, tspan, y0, h)
    t = tspan(1):h:tspan(2);
    y = zeros(size(t));
    y(1) = y0;
    for i = 1:length(t)-1
        y(i+1) = y(i) + h * feval(f, t(i), y(i));
    end
end

其中，f为微分方程的右端函数，tspan为时间区间，y0为初始值，h为步长。

改进欧拉法：

function [t, y] = improved_euler(f, tspan, y0, h)
    t = tspan(1):h:tspan(2);
    y = zeros(size(t));
    y(1) = y0;
    for i = 1:length(t)-1
        k1 = feval(f, t(i), y(i));
        k2 = feval(f, t(i+1), y(i) + h*k1);
        y(i+1) = y(i) + h/2 * (k1 + k2);
    end
end

其中，k1和k2为欧拉法中的斜率，feval函数用于计算函数f的值。