這是一個家庭作業問題。我將使用for循環和/或特效來繪製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
繪製在楓的函數(偏導數),我有以下功能 u := proc (x, t) options operator, arrow;
50+sum((100*n*Pi*sin(n*Pi)+100*cos(n*Pi)-100)*cos((1/10)*n*Pi*x)*exp(-(1/100)*n^2*Pi^2*t)/(n^2*Pi^2), n = 1 .. 20)
end proc;
p