使用Euler'；MATLAB中的s作图法

我在使用这个代码时遇到了一些问题。我的教授要求我们在MATLAB中创建一个函数“feuler.m”，以解决微分方程u′（t）=（2+2t）e^t和初始条件u（0）=0在区间[0,5]上给出的初值问题，该区间使用（正向）Euler方法绘制精确解和近似解

输入应为：

n

，间隔[0,5]应划分成的子间隔数

输出应为精确解和数值解的图形，并打印真实解和数值解之间的最大误差值

注意，精确解由u（t）=2et给出

到目前为止，我已经编写了代码：

function myeuler(N)
   t = linspace(0, 5, N+1)';
   ua = zeros(N+1,1);
   ue = 2*t.*exp(t);
   h = 5/N;
   A = zeros(N,N);
   A(2:N,1:N-1) = -eye(N-1);
   A = A + eye(N);
   b = h*(2+2*t(1:N)).*exp(t(1:N));
   b(1) = b(1) + ua(1);
   ua(2:N+1) = A\b;
   plot(t, ua, 'r', t, ue, 'g')
end

---

我不确定这是否正确。

您好，欢迎来到SO，请阅读。你的帖子缺少一个真正的问题。只是说你“不确定”并不是一个真正的编程问题。你为什么不实现Euler方法，在给定的情况下，它相当于函数值的累积和？