Feldlinien mit Mathematica < Mathematica < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) beantwortet  | | Datum: | 17:09 So 14.08.2005 | | Autor: | Martini |
Hallo!
Auf der Website www.wolfram.com kann man sich zu dem Buch von Tom Wickham-Jones "Mathematica Graphics: Techniques & Applications" (welches dort vorgestellt wird, daß ich aber leider nicht besitze) die zusätzliche Software "ExtendGraphics" downloaden. U.a. findet man dort das Package "FieldLines.m" mit dem Feldlinien von Vektorfeldern dargestellt werden können. Obwohl das Package einwandfrei geladen wird, funktioniert aber die entsprechende Funktion nicht richtig:
FieldLine[ {x, ex, x0}, {y, ey, y0}, {t, t1}] gibt bei mir leider mit jeder Eingabe eine Fehlermeldung.
Folgender Text zu dieser Funktion kann abgerufen werden: "FieldLine[ {x, ex, x0}, {y, ey, y0}, {t, t1}] will calculate the field line from the field {ex, ey}, starting at {x0,y0}, of length t1. x and y are the variables of the field and t is the variable of length down the trajectory."
Ich habe es z.B. mit einem ganz einfachen Vektorfeld ex=x, ey=y ausprobiert.
Was mache ich falsch???
(Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt..)
|
|
| |
|
Hallo Martini!
Ich habe ExtendGraphics nicht bei wolfram.com heruntergeladen - deshalb zunächst:
Die Package Version, über die ich verfüge, hat die Nummer 1.00 und stammt aus dem Sommer 1993. Im Quelltext FieldLines.m befindet sich dort in der Funktion FieldLine eine Zeile (bei mir Nr. 64): t2 = Part[ sol, 1, 0, 1, 2] ;. Das pfriemelt aus dem Vektor der beiden InterpolatingFunction-Objekte die obere Grenze, bis zu der NDSolve die Gleichungen lösen konnte. In neueren Mathematica-Versionen funktioniert dies aber mit t2 = Part[ sol, 1, 0, 1, 1, 2] ;.
Nach dieser Änderung ergab sich bei mir nach
| 1: | << ExtendGraphics'FieldLines'
| | 2: | fl=FieldLine[{x,1-Sin[Pi*t]-Sqrt[x^2+y^2],1},{y,Cos[x+y-t],-1/2},{t,1.8}];
| | 3: | Show[Graphics[fl],Axes->True]; |
[Dateianhang nicht öffentlich]
ein blödsinniges aber funktionierendes Beispiel ![[happy]](/images/smileys/happy.gif "[happy]")
Viel Erfolg,
Peter
Dateianhänge:
Anhang Nr. 1 (Typ: png) [nicht öffentlich]
|
|
|
|
