> Mein Problem ist das die Schleife (ich habe schon alle moeglichen Konstrukte mit
> Do, For, While durchgespielt) nur bis zu einem bestimmten Punkt durchgearbeitet
> wird und dann aaufhoert. Mathematica bleibt einfach stehen. Es wird nur ein Teil vom
> Hauptspeicher benutzt (also keine Auslastung) und die CPU wird stark belastet.
Das scheint entweder doch auf Speicherplatzprobleme oder nichtnumerische
große Zwischenresultate hinzuweisen. Sie können ja einmal die Schleife
"von Hand" ein paar Schritte weit laufen lassen und sich jeweils die
Zwischenergebnisse anschauen:
M = 4000;
v[0] = Table[....,{i,1,M+1}];
tau[0] = Table[...,{i,1,M+1}];
tet[0] = Table[...,{i,1,M+1}];
n=0
v[n+1] = v[n] + lam (tauy[n]+ k)
tau[n+1] = tau[n] + ......
tet[n+1] = v[n] + .... + tauyy[n] ......
n=1
etc...
Sie haben ins nicht gesagt, was tauy[n_], tauyy[n_],... sind. Diese sollten
wohl zu einer Zahl oder Liste von M (Gleitpunkt-) Zahlen evaluieren,
sonst sieht es nicht gut aus. Fehler bei der Definition der tauy
(oder der Anfangsbedingungen) können zum Aufbau von symbolischen Werten
führen, die sich dann in jedem Schleifendurchgang weiter auftürmen;
die geschilderten Symptome weisen auf so etwas hin.
Ein weiteres Problem, das aber offenbar hier nicht entscheidend ist,
liegt darin, daß Sie in jedem Schleifendurchgang neue Werte speichern;
dies ist nur dann sinnvoll, wenn Sie alle Zwischenergebnisse am Schluss
auch wirklich benötigen, sonst geht es effizienter mit
v = Table[....,{i,1,M+1}];
tau = Table[...,{i,1,M+1}];
tet = Table[...,{i,1,M+1}];
Do[
v1 = v + lam (tauy[n]+ k);
tau1 = tau + ......;
tet1 = v + .... + tauyy[n] ......;
v=v1; tau=tau1; tet=tet1;
, {n, 0, nstep}];
Womit Sie die alten Zwischenergebnisse überschreiben.
mit freundlichen Grüßen,
Roman Mäder
-----------------------------------------------------------------------
MathConsult Dr. R. Mäder Samstagernstrasse 58a
Mathematik- und Informatik-Beratung CH-8832 Wollerau
T: +41-1-687 4050 mailto:maeder@XXXXXXX.ch
F: +41-1-687 4054 http://www.mathconsult.ch/
-----------------------------------------------------------------------
|