Spacerat hat geschrieben:Zu dem Masseproblem mit Hammer, Feder und Mond... Auch Dinge, die sich rauskürzen, haben meistens Seiteneffekte, die in der weiteren Rechnung ins Gewicht fallen.
Auch du, mein Sohn Brutus? Sieh dir lieber den Code an bevor du den Trollen noch Wind in ihre zerissenen Segel bläst.
- Code: Alles auswählen
(* DGL *)
x1'[t] == vx1[t], y1'[t] == vy1[t], z1'[t] == vz1[t],
x2'[t] == vx2[t], y2'[t] == vy2[t], z2'[t] == vz2[t],
x3'[t] == vx3[t], y3'[t] == vy3[t], z3'[t] == vz3[t],
vx1'[t] == -((G m2 (x1[t] - x2[t]))/
Sqrt[((x1[t] - x2[t])^2 + (y1[t] - y2[t])^2 + (z1[t] -
z2[t])^2)^3]) - (G m3 (x1[t] - x3[t]))/
Sqrt[((x1[t] - x3[t])^2 + (y1[t] - y3[t])^2 + (z1[t] -
z3[t])^2)^3],
vy1'[t] == -((G m2 (y1[t] - y2[t]))/
Sqrt[((x1[t] - x2[t])^2 + (y1[t] - y2[t])^2 + (z1[t] -
z2[t])^2)^3]) - (G m3 (y1[t] - y3[t]))/
Sqrt[((x1[t] - x3[t])^2 + (y1[t] - y3[t])^2 + (z1[t] -
z3[t])^2)^3],
vz1'[t] == -((G m2 (z1[t] - z2[t]))/
Sqrt[((x1[t] - x2[t])^2 + (y1[t] - y2[t])^2 + (z1[t] -
z2[t])^2)^3]) - (G m3 (z1[t] - z3[t]))/
Sqrt[((x1[t] - x3[t])^2 + (y1[t] - y3[t])^2 + (z1[t] -
z3[t])^2)^3],
vx2'[t] == (G m1 (x1[t] - x2[t]))/
Sqrt[((x1[t] - x2[t])^2 + (y1[t] - y2[t])^2 + (z1[t] -
z2[t])^2)^3] - (G m3 (x2[t] - x3[t]))/
Sqrt[((x2[t] - x3[t])^2 + (y2[t] - y3[t])^2 + (z2[t] -
z3[t])^2)^3],
vy2'[t] == (G m1 (y1[t] - y2[t]))/
Sqrt[((x1[t] - x2[t])^2 + (y1[t] - y2[t])^2 + (z1[t] -
z2[t])^2)^3] - (G m3 (y2[t] - y3[t]))/
Sqrt[((x2[t] - x3[t])^2 + (y2[t] - y3[t])^2 + (z2[t] -
z3[t])^2)^3],
vz2'[t] == (G m1 (z1[t] - z2[t]))/
Sqrt[((x1[t] - x2[t])^2 + (y1[t] - y2[t])^2 + (z1[t] -
z2[t])^2)^3] - (G m3 (z2[t] - z3[t]))/
Sqrt[((x2[t] - x3[t])^2 + (y2[t] - y3[t])^2 + (z2[t] -
z3[t])^2)^3],
vx3'[t] == (G m1 (x1[t] - x3[t]))/
Sqrt[((x1[t] - x3[t])^2 + (y1[t] - y3[t])^2 + (z1[t] -
z3[t])^2)^3] + (G m2 (x2[t] - x3[t]))/
Sqrt[((x2[t] - x3[t])^2 + (y2[t] - y3[t])^2 + (z2[t] -
z3[t])^2)^3],
vy3'[t] == (G m1 (y1[t] - y3[t]))/
Sqrt[((x1[t] - x3[t])^2 + (y1[t] - y3[t])^2 + (z1[t] -
z3[t])^2)^3] + (G m2 (y2[t] - y3[t]))/
Sqrt[((x2[t] - x3[t])^2 + (y2[t] - y3[t])^2 + (z2[t] -
z3[t])^2)^3],
vz3'[t] == (G m1 (z1[t] - z3[t]))/
Sqrt[((x1[t] - x3[t])^2 + (y1[t] - y3[t])^2 + (z1[t] -
z3[t])^2)^3] + (G m2 (z2[t] - z3[t]))/
Sqrt[((x2[t] - x3[t])^2 + (y2[t] - y3[t])^2 + (z2[t] -
z3[t])^2)^3]
(* Startpositionen *)
x1[0] == Sin[0], y1[0] == Cos[0], z1[0] == 0,
x2[0] == Sin[2/3*p], y2[0] == Cos[2/3*p], z2[0] == 0,
x3[0] == Sin[4/3*p], y3[0] == Cos[4/3*p], z3[0] == 0
(* Anfangsgeschwindigkeiten *)
vx1[0] == 0, vy1[0] == 0, vz1[0] == 0,
vx2[0] == 0, vy2[0] == 0, vz2[0] == 0,
vx3[0] == 0, vy3[0] == 0, vz3[0] == 0
(* Dreikörperproblem Simulation, Yukterez - http://bit.ly/1aMbgGH *)
Wie du siehst wird Masse 1 von Masse 2 und Masse 3 beschleunigt, nicht aber von sich selber. Masse 2 wird dafür von Masse 1 und Masse 3 angezogen, und Masse 3 von Masse 1 und Masse 2. So einfach ist das!
, 

hat geschrieben:Ja nee, is klar.
hat geschrieben:Du hast bis heute noch nicht einmal die Fragestellung verstanden
hat geschrieben:Du selber kommst immer mit Dingen die vernachlässigbar klein sind, damit Deine Idiotenrechnungen aufgehen
,