欧拉法matlab程序代码

以下是欧拉法的 Matlab 程序代码：

function [t, y] = euler(f, tspan, y0, h)
% 欧拉法求解常微分方程
% f: 函数句柄，t 和 y 为自变量和因变量，返回值为 y' 的值
% tspan: 时间区间，格式为 [t0, tf]
% y0: 初值
% h: 步长
t0 = tspan(1);
tf = tspan(2);
t = t0:h:tf;
y = zeros(size(t));
y(1) = y0;
for i = 1:length(t)-1
y(i+1) = y(i) + h * f(t(i), y(i));
end
end

希望对你有所帮助！