|
Hallo, also das sieht mir nach einer singularit"at in Deinen Gleichungen aus. Wenn Du mal die explixite Form f"ur y[x]<=1 berechnest, siehst Du, dass ein term ... /y[x] auftaucht. In Kurzform ergibt sich etwas wie y[x]^2*( _ (_ + _/y[x])^10) und der y[x]^2 Term kann die 1/y[x]^10 nicht kompensieren. Deine Parameter sind nun so gew"ahlt das y[x] tats"achlich gegen 0 geht. Ob das sinnvoll ist weiss ich nat"urlich nicht. Da ich auch keinen Schimmer von den Parametern habe kann ich Dir auch da keine Tip geben. Da die meisten numerischen Verfahren f"ur Anfangswertprobleme mit Polynom-Approximationen arbeiten also lokal y[x+h]=y[x]+y'[x]*h+y''[x]*h^2 ... ann"ahern, kann man damit keine Polstelle "uberspringen. Der Anfangswertl"oser macht immer kleinere Schritte wenn er auf einen Pol st"osst und gibt irgendwann auf. Deshalb geht Deine Interpolation[] nicht "uber das volle Zeitinterval. Gruss Jens cz wrote: > > Lieber Jens-Peer Kuska, > > anbei noch einmal eine Kombination von NDSolve und Fallunterscheidung. > Ich muß da noch irgendwo einen wahrscheinlich wieder trivialen Fehler > haben, den ich aber nicht finde und der zu der angegebenen Fehlermeldung > führt. Könntest Du vielleicht so nett sein, da einfach mal drüber > schauen? Ich habe das File als Attachment beigefügt. > > Vielen Dank. > > Christian Zupanc > Technische Chemie > Ruhruniversität Bochum |