maple - 微分方程式のシステムでの割り当て

Question

2つの微分方程式のシステムを数値的に解くとき：

s1:=diff(n[Di](t), t)=...;
s2:=diff(n[T](t), t)=...;
ics:={...};   #initial condition.
sys := {s1, s2, ics}: 
sol:=dsolve(sys,numeric);

「t」に関しては、「t = 4」の解は（たとえば）sol（4）の形式になります。

[t=4, n1(t)=const1, n2(t)=const2].

ここで、n1（t）またはn2（t）を含む別の方程式「p」のすべての「t」にn1（t）およびn2（t）の値を使用するにはどうすればよいですか（例：{p = a + n1（t）* n2（t）+ f（t）}、ここで「a」と「f（t）」が定義されています）、「t」の間隔で「p」をプロットするには？

Answer 1

おそらく例が役立つでしょう。

s1:=diff(n1(t), t)=sin(t);
s2:=diff(n2(t), t)=cos(t);
ics:={n1(0)=1,n2(0)=0}:
sys := {s1, s2} union ics:
sol:=dsolve(sys,numeric,output=listprocedure);
N1:=eval(n1(t),sol);
N2:=eval(n2(t),sol);
N1(0);
N1(1.1);
N2(0);
N2(1.1);
p := 11.3 + N1*N2 + tan:
plot(p, 0..1.1);

常微分方程式のプロットに特化したDEplotルーチンをご覧ください。