DgL's in Maple lösen < Maple < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) beantwortet  | | Datum: | 22:58 Di 24.10.2006 | | Autor: | Thanas |
Hallo
ich möchte mit maple zwei voneinader abängige DGL's lösen.
ich habe dazu folgenden Ansatz:
restart:
with(plots);
dgl:=D(N[1])(t)=-r[1]*N[1](t)-beta[1]*N[1](t), D(N[2])(t)=r[2]*N[1](t)-lambda[2]*N[2](t)-beta[2]*N[2](t);
anf:=N[1](0)=N[1,0],N[2](0)=N[2,0];
sol:=dsolve({dgl,anf},{N[1](t),N[2](t)});
simplify(sol);
assign(sol);
N[1](t);
N[2](t);
N[1,plot]:=subs(N[1,0]=10,r[1]=0.000001,beta[1]=1,N[1](t));
plot(N[1,plot](t), t=0..20);
N[2,plot]:=subs(N[1,0]=10,N[2,0]=10,r[1]=0.000001,r[2]=0.0001,lambda[2]=0.00047,beta[1]=1,beta[2]=1,N[2](t));
plot(N[2,plot](t), t=0..20);
Das ganze funktioniert auch soweit. Nun mein Problem. Die Gleichung N[2](t) (vorletzter Schritt) wird nicht berechnet und somit wird auch kein Plot durchgeführt. Warum ? Wenn ich N[2](t) vor N1[t] berechen lasse wird N[1](t) nicht berechet. Hat jemand da ne Idee?
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
|
|
| |
|
Hallo,
es sind wohl die plot-Variablen, die Ärger machen. Am einfachsten wäre es, sie umzubenennen, das Maple anscheinend bei der Zuweisung
N[1,plot]:=...
die alten N-Belegungen löscht. Wie das genau aussieht, weiß ich nicht, aber mit Nplot[1] und Nplot[2] funktioniert es.
Gruß
Martin
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 18:00 Mi 25.10.2006 | | Autor: | Thanas |
Danke. So klappt es. Kleine Ursache Große Wirkung.
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig | | Datum: | 20:46 Do 26.10.2006 | | Autor: | Thanas |
Wenn ich obiges 'DGL-System' um 2 DGL's erweitere :
D(N[3])(t)=r[2]*N[2](t)-lambda[3]*N[3](t)-beta[3]*N[3](t) und D(N[4])(t)=r[3]*N[3](t)-lambda[4]*N[4](t)-beta[4]*N[4](t)
dann berechnet er mir auch noch die 3. DGL, die vierte aber wird nicht vollständig gelöst.
Woran könnte den das liegen
|
|
|
| |
|
| Status: |
(Frage) überfällig  | | Datum: | 20:49 Do 26.10.2006 | | Autor: | Thanas |
Wenn ich obiges 'DGL-System' um 2 DGL's erweitere :
D(N[3])(t)=r[2]*N[2](t)-lambda[3]*N[3](t)-beta[3]*N[3](t) und D(N[4])(t)=r[3]*N[3](t)-lambda[4]*N[4](t)-beta[4]*N[4](t)
dann berechnet er mir auch noch die 3. DGL, die vierte aber wird nicht vollständig gelöst.
Woran könnte den das liegen ?
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig | | Datum: | 21:20 Sa 28.10.2006 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
