无detools的maple plot方向场

无detools的maple plot方向场,plot,maple,Plot,Maple,这是一个家庭作业问题。我将使用for循环和or过程绘制ode y'=x*y/2的方向场。我有以下代码: X := [seq(.1*x, x = -20 .. 20)]; Y := [seq(y, y = -2 .. 2, .1)]; P := Array([seq(0, x = 0 .. 41^2-1)]); with(plots); for k to 41 do for j to 41 do slope := (1/2)*X[k]*Y[j]; Xp := [X[k], X[k]+.1]

这是一个家庭作业问题。我将使用for循环和or过程绘制ode y'=x*y/2的方向场。我有以下代码:

X := [seq(.1*x, x = -20 .. 20)]; 
Y := [seq(y, y = -2 .. 2, .1)]; 
P := Array([seq(0, x = 0 .. 41^2-1)]);
with(plots); 
for k to 41 do for j to 41 do
slope := (1/2)*X[k]*Y[j]; 
Xp := [X[k], X[k]+.1]; Yp := [Y[j], Y[j]+0.1*slope]; 
P[41*(k-1)+j] := plot(Xp, Yp) end do end do;
display(P);
我了解很多情节,使用detools会更容易,但正如前面所说的,它是用来做家庭作业的。我得到这样的东西:

如果执行
display(p[1..2])
,您可以看到为什么
显示(p)
看起来是错误的-它是一个“打印窗口窗口窗口”

您需要告诉Maple,p不是一个对象,但显示应为1681个绘图,例如序列功能

display(seq(P[i], i = 1..numelems(P)));

如果执行
display(p[1..2])
,您可以看到为什么
显示(p)
看起来是错误的-它是一个“打印窗口窗口窗口”

您需要告诉Maple,p不是一个对象,但显示应为1681个绘图,例如序列功能

display(seq(P[i], i = 1..numelems(P)));

建议:注意不要与分号(
)一致。此外,缩进代码使其更具可读性。:)建议:注意不要与分号(
)一致。此外,缩进代码使其更具可读性。:)