MATLAB中用欧拉法求解微分方程组

欧拉法是一种常用的数值求解微分方程的方法，可以在MATLAB中使用。下面是一个使用欧拉法求解微分方程组的示例：

% 定义微分方程组
function dydt = myODE(t, y)
    dydt = zeros(2,1);
    dydt(1) = 2*(1-0.01*y(1))*y(1);
    dydt(2) = -0.5*y(2);
end

% 设置初始条件和时间步长
t0 = 0;
y0 = [20; 10];
h = 0.1;

% 使用欧拉法求解微分方程组
t = t0:h:10;
y = zeros(length(t), 2);
y(1,:) = y0;
for i = 1:length(t)-1
    y(i+1,:) = y(i,:) + h*myODE(t(i), y(i,:)');
end

% 绘制结果
plot(t, y(:,1), 'r-', t, y(:,2), 'b-');
legend('y1', 'y2');
xlabel('t');
ylabel('y');