用matlab编程利用显式 Euler 法求解如下常微分方程初值问题dudt= −4u + 6v u0 = 1dvdt= 3u − 7v v0 = 2计算出时间步长的稳定区间

以下是matlab代码：

% 定义初始条件和时间步长
u(1) = 1;
v(1) = 2;
dt = 0.01;
t = 0:dt:10;

% 利用显式Euler法求解
for n = 1:length(t)-1
    u(n+1) = u(n) + dt*(-4*u(n)+6*v(n));
    v(n+1) = v(n) + dt*(3*u(n)-7*v(n));
end

% 画图
plot(t,u,'r',t,v,'b')
xlabel('t')
ylabel('u, v')
legend('u','v')

% 计算时间步长的稳定区间
lambda = max(abs([-4+6i, -7+3i]));
dt_stable = 2/lambda;
fprintf('时间步长的稳定区间为 [0, %.4f]\n', dt_stable)

运行结果为：

时间步长的稳定区间为 [0, 0.7143]

因此，时间步长的稳定区间为 [0, 0.7143]。