назначения в системе дифференциальных уравнений

Question

, когда я численно решить систему двух дифференциальных уравнений:

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

по отношению к «Т», то решение (например) для «Т = 4» из форма, золь (4):

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

сейчас, как можно использовать значения n1 (Т) и n2 (Т) для всех «т» 'ы в другом уравнении, а именно „р“, в котором участвуют n1 (t) или n2 (t) (например: {p = a + n1 (t) * n2 (t) + f (t)}, где определены «a» и «f (t)»), и сюжет «p» для интервала «t»?

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, которая специализируется на построении решений ode.