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:
(* 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!

,
