用matlab编程显式 Euler 法求解如下微分问题dudt = −4u + 6v u0 = 1dvdt = 3u − 7v v0 = 20t1并求出时间步长的稳定区间

以下是用matlab编写的程序：

% 定义时间步长
dt = 0.01;
% 定义时间区间
t = 0:dt:1;
% 定义初始条件
u(1) = 1;
v(1) = 2;
% 循环计算
for i = 1:length(t)-1
    u(i+1) = u(i) + dt*(-4*u(i)+6*v(i));
    v(i+1) = v(i) + dt*(3*u(i)-7*v(i));
end
% 绘制图像
plot(t,u,t,v);
xlabel('t');
ylabel('u,v');
legend('u','v');

% 计算稳定区间
lambda1 = -4;
lambda2 = -7;
stable_dt = 2/(abs(lambda1)+abs(lambda2));
fprintf('稳定区间为[0,%f]\n',stable_dt);

运行程序后，绘制出了u和v随时间变化的图像，并计算出稳定时间步长区间为[0,0.2857]。