如何用 maple 12 绘制和求解这个微分方程？答案

【问题标题】：How to plot and solve this differential equation with maple 12?如何用 maple 12 绘制和求解这个微分方程？
【发布时间】：2013-04-07 10:16:20
【问题描述】：

考虑初始条件为 $y\left(0\right)=1$ 的微分方程 $y^{\prime}=y-2$ 。

a) 使用欧拉法以 4 步大小 0.2 估计 $y\left(0.8\right)$

我知道如何手工操作；但是，我安装了 Maple 12，并试图弄清楚如何使用 Maple 执行此操作，然后制作一个图表来显示该功能的每个步骤。有什么建议么。我曾尝试查看 mapleprimes，但它一直指向我没有的新版本 maplesoft 的功能。

我发布了这个问题以用作模型，因为我已经手动解决了这个问题，它将帮助我为其他微分方程编辑它。

ps。我希望这是问这个问题的合适地方，如果不是，请告诉我哪里会更好。

【问题讨论】：

标签： maple differential-equations

【解决方案1】：

restart:

sys := diff(y(x),x) = y(x) - 2;
IC := y(0) = 1;

sol := dsolve({sys,IC}, numeric,
          output = listprocedure,
          method = classical[foreuler], stepsize = 0.2):
yest := eval(y(x),sol):

for i from 1 to 4 do
  yest(0.2*i);
end do;

Peul:=plots:-odeplot(sol,x=0..1):
Peul;

exactsol:=dsolve({sys,IC});
yexact:=eval(y(x),exactsol);

for i from 1 to 4 do
  eval(yexact,x=0.2*i);
end do;

Pexact:=DEtools[DEplot](sys,y(x),x=0..1,[[y(0)=1]]):

plots:-display([Pexact,Peul],color=[green,red]);

【讨论】：