hat geschrieben:Ja, was für ein törichter Schusselfehler - besten Dank, ich habe es korrigiert.
Genau im selben Moment. Dann waren das eh die selben Fehler die ich auch gefunden habe, dann dürfte das obere Ergebnis passen und Amalthea der im Durchschnitt nächste Mond für alle anderen Monde sein.
hat geschrieben:Kann es sein, dass das Minimum gerade auf die Roche-Grenze fällt ?
Das glaube ich kaum, denn die hängt ja von der Ausdehnung bzw. Dichte des Körpers selber ab, und die darf nicht in unsere Rechnung einfließen.
hat geschrieben:Dateneingabe ist abgeschlossen
- Code: Alles auswählen
(* |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| *)
(* ||| Mathematica Syntax || yukterez.net || n Body Newtonian Mass & Charge Simulator ||| *)
(* |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| *)
ClearAll["Global`*"]; ClearAll["Local`*"];
Amp = 1; kg = 1; m = 1; sek = 1; km = 1000 m; (* SI Einheiten *)
mt1 = {"StiffnessSwitching", Method-> {"ExplicitRungeKutta", Automatic}};
mt2 = {"ImplicitRungeKutta", "DifferenceOrder"-> 20};
mt3 = {"EquationSimplification"-> "Residual"};
mt0 = Automatic;
mta = mt0;
wp = MachinePrecision;
Tmax = -3 yr;
tMax = Min[Tmax, plunge];
(* Konstanten *)
G = 667384/10^16 m^3/kg/sek^2;
ε0 = 8854187817*^-21 Amp^2 sek^4/kg/m^3;
c = 299792458 m/sek;
Au = 149597870700 m;
dy = 24*3600 sek;
yr = 36525*dy/100;
(* Ephemeriden vom 26.03.2019, 0:00:00 TDB *)
n1 = Jupiter;
m1 =+1.89813*^27 kg;
q1 =+0 Amp sek;
x1x=-1.548288165488712*^+00 Au;
y1y=-5.086853499541566*^+00 Au;
z1z=+5.573068942555429*^-02 Au;
v1x=+7.129123597183238*^-03 Au/dy;
v1y=-1.837734375671549*^-03 Au/dy;
v1z=-1.518302302430324*^-04 Au/dy;
n2 = Metis;
m2 =+3.6*^16 kg;
q2 =+0 Amp sek;
x2x=-x1x-1.548043607670377*^+00 Au;
y2y=-y1y-5.087672518468227*^+00 Au;
z2z=-z1z+5.570507524119116*^-02 Au;
v2x=-v1x+2.460878825544765*^-02 Au/dy;
v2y=-v1y+3.362398066415406*^-03 Au/dy;
v2z=-v1z+2.871817613272839*^-04 Au/dy;
n3 = Adrastea;
m3 =+2.0*^15 kg;
q3 =+0 Amp sek;
x3x=-x1x-1.548002955883037*^+00 Au;
y3y=-y1y-5.087666738351646*^+00 Au;
z3z=-z1z+5.570592786396499*^-02 Au;
v3x=-v1x+2.426651595395510*^-02 Au/dy;
v3y=-v1y+4.161688062837289*^-03 Au/dy;
v3z=-v1z+3.127719871486290*^-04 Au/dy;
n4 = Amalthea;
m4 =+2.08*^18 kg;
q4 =+0 Amp sek;
x4x=-x1x-1.549495495563100*^+00 Au;
y4y=-y1y-5.086770056173906*^+00 Au;
z4z=-z1z+5.571260752315928*^-02 Au;
v4x=-v1x+6.038463256460282*^-03 Au/dy;
v4y=-v1y-1.710080909791476*^-02 Au/dy;
v4z=-v1z-8.055082131342978*^-04 Au/dy;
n5 = Thebe;
m5 =+4.3*^17 kg;
q5 =+0 Amp sek;
x5x=-x1x-1.547341116060141*^+00 Au;
y5y=-y1y-5.088008258680005*^+00 Au;
z5z=-z1z+5.570988948392359*^-02 Au;
v5x=-v1x+1.758839294391254*^-02 Au/dy;
v5y=-v1y+7.012149989548776*^-03 Au/dy;
v5z=-v1z+5.666948041750221*^-04 Au/dy;
n6 = Io;
m6 =+8.933*^22 kg;
q6 =+0 Amp sek;
x6x=-x1x-1.545691479501977*^+00 Au;
y6y=-y1y-5.085739318152110*^+00 Au;
z6z=-z1z+5.580826914459231*^-02 Au;
v6x=-v1x+3.215551992239203*^-03 Au/dy;
v6y=-v1y+7.346089625645640*^-03 Au/dy;
v6z=-v1z+1.263026713707490*^-04 Au/dy;
n7 = Europa;
m7 =+4.797*^22 kg;
q7 =+0 Amp sek;
x7x=-x1x-1.551669842589193*^+00 Au;
y7y=-y1y-5.083931466933994*^+00 Au;
z7z=-z1z+5.578340888798599*^-02 Au;
v7x=-v1x+1.870558078587358*^-03 Au/dy;
v7y=-v1y-7.809888623581773*^-03 Au/dy;
v7z=-v1z-5.075203619950251*^-04 Au/dy;
n8 = Ganymede;
m8 =+1.482*^23 kg;
q8 =+0 Amp sek;
x8x=-x1x-1.550245613643152*^+00 Au;
y8y=-y1y-5.079958581591128*^+00 Au;
z8z=-z1z+5.596555163834176*^-02 Au;
v8x=-v1x+1.096081003986504*^-03 Au/dy;
v8y=-v1y-3.539865361898701*^-03 Au/dy;
v8z=-v1z-2.973495686088762*^-04 Au/dy;
n9 = Kallisto;
m9 =+1.076*^23 kg;
q9 =+0 Amp sek;
x9x=-x1x-1.536864188635359*^+00 Au;
y9y=-y1y-5.081792556406656*^+00 Au;
z9z=-z1z+5.604371626803566*^-02 Au;
v9x=-v1x+5.203319956094230*^-03 Au/dy;
v9y=-v1y+2.526469914758841*^-03 Au/dy;
v9z=-v1z-3.995285941360583*^-05 Au/dy;
(* Zentrierung auf Jupiter *)
x1x=0;
y1y=0;
z1z=0;
v1x=0;
v1y=0;
v1z=0;
(* Differentialgleichung *)
nds=NDSolve[{
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],
x4'[t] == vx4[t], y4'[t] == vy4[t], z4'[t] == vz4[t],
x5'[t] == vx5[t], y5'[t] == vy5[t], z5'[t] == vz5[t],
x6'[t] == vx6[t], y6'[t] == vy6[t], z6'[t] == vz6[t],
x7'[t] == vx7[t], y7'[t] == vy7[t], z7'[t] == vz7[t],
x8'[t] == vx8[t], y8'[t] == vy8[t], z8'[t] == vz8[t],
x9'[t] == vx9[t], y9'[t] == vy9[t], z9'[t] == vz9[t],
vx1'[t] ==
(G m2 (x2[t]-x1[t]))/Sqrt[((x2[t]-x1[t])^2+(y2[t]-y1[t])^2+(z2[t]-z1[t])^2)^3]+
(G m3 (x3[t]-x1[t]))/Sqrt[((x3[t]-x1[t])^2+(y3[t]-y1[t])^2+(z3[t]-z1[t])^2)^3]+
(G m4 (x4[t]-x1[t]))/Sqrt[((x4[t]-x1[t])^2+(y4[t]-y1[t])^2+(z4[t]-z1[t])^2)^3]+
(G m5 (x5[t]-x1[t]))/Sqrt[((x5[t]-x1[t])^2+(y5[t]-y1[t])^2+(z5[t]-z1[t])^2)^3]+
(G m6 (x6[t]-x1[t]))/Sqrt[((x6[t]-x1[t])^2+(y6[t]-y1[t])^2+(z6[t]-z1[t])^2)^3]+
(G m7 (x7[t]-x1[t]))/Sqrt[((x7[t]-x1[t])^2+(y7[t]-y1[t])^2+(z7[t]-z1[t])^2)^3]+
(G m8 (x8[t]-x1[t]))/Sqrt[((x8[t]-x1[t])^2+(y8[t]-y1[t])^2+(z8[t]-z1[t])^2)^3]+
(G m9 (x9[t]-x1[t]))/Sqrt[((x9[t]-x1[t])^2+(y9[t]-y1[t])^2+(z9[t]-z1[t])^2)^3]+
If[q1 == 0, 0,
(-q1*q2/(4Pi ε0 )/m1 (x2[t]-x1[t]))/Sqrt[((x2[t]-x1[t])^2+(y2[t]-y1[t])^2+(z2[t]-z1[t])^2)^3]+
(-q1*q3/(4Pi ε0 )/m1 (x3[t]-x1[t]))/Sqrt[((x3[t]-x1[t])^2+(y3[t]-y1[t])^2+(z3[t]-z1[t])^2)^3]+
(-q1*q4/(4Pi ε0 )/m1 (x4[t]-x1[t]))/Sqrt[((x4[t]-x1[t])^2+(y4[t]-y1[t])^2+(z4[t]-z1[t])^2)^3]+
(-q1*q5/(4Pi ε0 )/m1 (x5[t]-x1[t]))/Sqrt[((x5[t]-x1[t])^2+(y5[t]-y1[t])^2+(z5[t]-z1[t])^2)^3]+
(-q1*q6/(4Pi ε0 )/m1 (x6[t]-x1[t]))/Sqrt[((x6[t]-x1[t])^2+(y6[t]-y1[t])^2+(z6[t]-z1[t])^2)^3]+
(-q1*q7/(4Pi ε0 )/m1 (x7[t]-x1[t]))/Sqrt[((x7[t]-x1[t])^2+(y7[t]-y1[t])^2+(z7[t]-z1[t])^2)^3]+
(-q1*q8/(4Pi ε0 )/m1 (x8[t]-x1[t]))/Sqrt[((x8[t]-x1[t])^2+(y8[t]-y1[t])^2+(z8[t]-z1[t])^2)^3]+
(-q1*q9/(4Pi ε0 )/m1 (x9[t]-x1[t]))/Sqrt[((x9[t]-x1[t])^2+(y9[t]-y1[t])^2+(z9[t]-z1[t])^2)^3]],
vy1'[t] ==
(G m2 (y2[t]-y1[t]))/Sqrt[((x2[t]-x1[t])^2+(y2[t]-y1[t])^2+(z2[t]-z1[t])^2)^3]+
(G m3 (y3[t]-y1[t]))/Sqrt[((x3[t]-x1[t])^2+(y3[t]-y1[t])^2+(z3[t]-z1[t])^2)^3]+
(G m4 (y4[t]-y1[t]))/Sqrt[((x4[t]-x1[t])^2+(y4[t]-y1[t])^2+(z4[t]-z1[t])^2)^3]+
(G m5 (y5[t]-y1[t]))/Sqrt[((x5[t]-x1[t])^2+(y5[t]-y1[t])^2+(z5[t]-z1[t])^2)^3]+
(G m6 (y6[t]-y1[t]))/Sqrt[((x6[t]-x1[t])^2+(y6[t]-y1[t])^2+(z6[t]-z1[t])^2)^3]+
(G m7 (y7[t]-y1[t]))/Sqrt[((x7[t]-x1[t])^2+(y7[t]-y1[t])^2+(z7[t]-z1[t])^2)^3]+
(G m8 (y8[t]-y1[t]))/Sqrt[((x8[t]-x1[t])^2+(y8[t]-y1[t])^2+(z8[t]-z1[t])^2)^3]+
(G m9 (y9[t]-y1[t]))/Sqrt[((x9[t]-x1[t])^2+(y9[t]-y1[t])^2+(z9[t]-z1[t])^2)^3]+
If[q1 == 0, 0,
(-q1*q2/(4Pi ε0 )/m1 (y2[t]-y1[t]))/Sqrt[((x2[t]-x1[t])^2+(y2[t]-y1[t])^2+(z2[t]-z1[t])^2)^3]+
(-q1*q3/(4Pi ε0 )/m1 (y3[t]-y1[t]))/Sqrt[((x3[t]-x1[t])^2+(y3[t]-y1[t])^2+(z3[t]-z1[t])^2)^3]+
(-q1*q4/(4Pi ε0 )/m1 (y4[t]-y1[t]))/Sqrt[((x4[t]-x1[t])^2+(y4[t]-y1[t])^2+(z4[t]-z1[t])^2)^3]+
(-q1*q5/(4Pi ε0 )/m1 (y5[t]-y1[t]))/Sqrt[((x5[t]-x1[t])^2+(y5[t]-y1[t])^2+(z5[t]-z1[t])^2)^3]+
(-q1*q6/(4Pi ε0 )/m1 (y6[t]-y1[t]))/Sqrt[((x6[t]-x1[t])^2+(y6[t]-y1[t])^2+(z6[t]-z1[t])^2)^3]+
(-q1*q7/(4Pi ε0 )/m1 (y7[t]-y1[t]))/Sqrt[((x7[t]-x1[t])^2+(y7[t]-y1[t])^2+(z7[t]-z1[t])^2)^3]+
(-q1*q8/(4Pi ε0 )/m1 (y8[t]-y1[t]))/Sqrt[((x8[t]-x1[t])^2+(y8[t]-y1[t])^2+(z8[t]-z1[t])^2)^3]+
(-q1*q9/(4Pi ε0 )/m1 (y9[t]-y1[t]))/Sqrt[((x9[t]-x1[t])^2+(y9[t]-y1[t])^2+(z9[t]-z1[t])^2)^3]],
vz1'[t] ==
(G m2 (z2[t]-z1[t]))/Sqrt[((x2[t]-x1[t])^2+(y2[t]-y1[t])^2+(z2[t]-z1[t])^2)^3]+
(G m3 (z3[t]-z1[t]))/Sqrt[((x3[t]-x1[t])^2+(y3[t]-y1[t])^2+(z3[t]-z1[t])^2)^3]+
(G m4 (z4[t]-z1[t]))/Sqrt[((x4[t]-x1[t])^2+(y4[t]-y1[t])^2+(z4[t]-z1[t])^2)^3]+
(G m5 (z5[t]-z1[t]))/Sqrt[((x5[t]-x1[t])^2+(y5[t]-y1[t])^2+(z5[t]-z1[t])^2)^3]+
(G m6 (z6[t]-z1[t]))/Sqrt[((x6[t]-x1[t])^2+(y6[t]-y1[t])^2+(z6[t]-z1[t])^2)^3]+
(G m7 (z7[t]-z1[t]))/Sqrt[((x7[t]-x1[t])^2+(y7[t]-y1[t])^2+(z7[t]-z1[t])^2)^3]+
(G m8 (z8[t]-z1[t]))/Sqrt[((x8[t]-x1[t])^2+(y8[t]-y1[t])^2+(z8[t]-z1[t])^2)^3]+
(G m9 (z9[t]-z1[t]))/Sqrt[((x9[t]-x1[t])^2+(y9[t]-y1[t])^2+(z9[t]-z1[t])^2)^3]+
If[q1 == 0, 0,
(-q1*q2/(4Pi ε0 )/m1 (z2[t]-z1[t]))/Sqrt[((x2[t]-x1[t])^2+(y2[t]-y1[t])^2+(z2[t]-z1[t])^2)^3]+
(-q1*q3/(4Pi ε0 )/m1 (z3[t]-z1[t]))/Sqrt[((x3[t]-x1[t])^2+(y3[t]-y1[t])^2+(z3[t]-z1[t])^2)^3]+
(-q1*q4/(4Pi ε0 )/m1 (z4[t]-z1[t]))/Sqrt[((x4[t]-x1[t])^2+(y4[t]-y1[t])^2+(z4[t]-z1[t])^2)^3]+
(-q1*q5/(4Pi ε0 )/m1 (z5[t]-z1[t]))/Sqrt[((x5[t]-x1[t])^2+(y5[t]-y1[t])^2+(z5[t]-z1[t])^2)^3]+
(-q1*q6/(4Pi ε0 )/m1 (z6[t]-z1[t]))/Sqrt[((x6[t]-x1[t])^2+(y6[t]-y1[t])^2+(z6[t]-z1[t])^2)^3]+
(-q1*q7/(4Pi ε0 )/m1 (z7[t]-z1[t]))/Sqrt[((x7[t]-x1[t])^2+(y7[t]-y1[t])^2+(z7[t]-z1[t])^2)^3]+
(-q1*q8/(4Pi ε0 )/m1 (z8[t]-z1[t]))/Sqrt[((x8[t]-x1[t])^2+(y8[t]-y1[t])^2+(z8[t]-z1[t])^2)^3]+
(-q1*q9/(4Pi ε0 )/m1 (z9[t]-z1[t]))/Sqrt[((x9[t]-x1[t])^2+(y9[t]-y1[t])^2+(z9[t]-z1[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 (x3[t]-x2[t]))/Sqrt[((x3[t]-x2[t])^2+(y3[t]-y2[t])^2+(z3[t]-z2[t])^2)^3]+
(G m4 (x4[t]-x2[t]))/Sqrt[((x4[t]-x2[t])^2+(y4[t]-y2[t])^2+(z4[t]-z2[t])^2)^3]+
(G m5 (x5[t]-x2[t]))/Sqrt[((x5[t]-x2[t])^2+(y5[t]-y2[t])^2+(z5[t]-z2[t])^2)^3]+
(G m6 (x6[t]-x2[t]))/Sqrt[((x6[t]-x2[t])^2+(y6[t]-y2[t])^2+(z6[t]-z2[t])^2)^3]+
(G m7 (x7[t]-x2[t]))/Sqrt[((x7[t]-x2[t])^2+(y7[t]-y2[t])^2+(z7[t]-z2[t])^2)^3]+
(G m8 (x8[t]-x2[t]))/Sqrt[((x8[t]-x2[t])^2+(y8[t]-y2[t])^2+(z8[t]-z2[t])^2)^3]+
(G m9 (x9[t]-x2[t]))/Sqrt[((x9[t]-x2[t])^2+(y9[t]-y2[t])^2+(z9[t]-z2[t])^2)^3]+
If[q2 == 0, 0,
(-q2*q1/(4Pi ε0 )/m2 (x1[t]-x2[t]))/Sqrt[((x1[t]-x2[t])^2+(y1[t]-y2[t])^2+(z1[t]-z2[t])^2)^3]+
(-q2*q3/(4Pi ε0 )/m2 (x3[t]-x2[t]))/Sqrt[((x3[t]-x2[t])^2+(y3[t]-y2[t])^2+(z3[t]-z2[t])^2)^3]+
(-q2*q4/(4Pi ε0 )/m2 (x4[t]-x2[t]))/Sqrt[((x4[t]-x2[t])^2+(y4[t]-y2[t])^2+(z4[t]-z2[t])^2)^3]+
(-q2*q5/(4Pi ε0 )/m2 (x5[t]-x2[t]))/Sqrt[((x5[t]-x2[t])^2+(y5[t]-y2[t])^2+(z5[t]-z2[t])^2)^3]+
(-q2*q6/(4Pi ε0 )/m2 (x6[t]-x2[t]))/Sqrt[((x6[t]-x2[t])^2+(y6[t]-y2[t])^2+(z6[t]-z2[t])^2)^3]+
(-q2*q7/(4Pi ε0 )/m2 (x7[t]-x2[t]))/Sqrt[((x7[t]-x2[t])^2+(y7[t]-y2[t])^2+(z7[t]-z2[t])^2)^3]+
(-q2*q8/(4Pi ε0 )/m2 (x8[t]-x2[t]))/Sqrt[((x8[t]-x2[t])^2+(y8[t]-y2[t])^2+(z8[t]-z2[t])^2)^3]+
(-q2*q9/(4Pi ε0 )/m2 (x9[t]-x2[t]))/Sqrt[((x9[t]-x2[t])^2+(y9[t]-y2[t])^2+(z9[t]-z2[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 (y3[t]-y2[t]))/Sqrt[((x3[t]-x2[t])^2+(y3[t]-y2[t])^2+(z3[t]-z2[t])^2)^3]+
(G m4 (y4[t]-y2[t]))/Sqrt[((x4[t]-x2[t])^2+(y4[t]-y2[t])^2+(z4[t]-z2[t])^2)^3]+
(G m5 (y5[t]-y2[t]))/Sqrt[((x5[t]-x2[t])^2+(y5[t]-y2[t])^2+(z5[t]-z2[t])^2)^3]+
(G m6 (y6[t]-y2[t]))/Sqrt[((x6[t]-x2[t])^2+(y6[t]-y2[t])^2+(z6[t]-z2[t])^2)^3]+
(G m7 (y7[t]-y2[t]))/Sqrt[((x7[t]-x2[t])^2+(y7[t]-y2[t])^2+(z7[t]-z2[t])^2)^3]+
(G m8 (y8[t]-y2[t]))/Sqrt[((x8[t]-x2[t])^2+(y8[t]-y2[t])^2+(z8[t]-z2[t])^2)^3]+
(G m9 (y9[t]-y2[t]))/Sqrt[((x9[t]-x2[t])^2+(y9[t]-y2[t])^2+(z9[t]-z2[t])^2)^3]+
If[q2 == 0, 0,
(-q2*q1/(4Pi ε0 )/m2 (y1[t]-y2[t]))/Sqrt[((x1[t]-x2[t])^2+(y1[t]-y2[t])^2+(z1[t]-z2[t])^2)^3]+
(-q2*q3/(4Pi ε0 )/m2 (y3[t]-y2[t]))/Sqrt[((x3[t]-x2[t])^2+(y3[t]-y2[t])^2+(z3[t]-z2[t])^2)^3]+
(-q2*q4/(4Pi ε0 )/m2 (y4[t]-y2[t]))/Sqrt[((x4[t]-x2[t])^2+(y4[t]-y2[t])^2+(z4[t]-z2[t])^2)^3]+
(-q2*q5/(4Pi ε0 )/m2 (y5[t]-y2[t]))/Sqrt[((x5[t]-x2[t])^2+(y5[t]-y2[t])^2+(z5[t]-z2[t])^2)^3]+
(-q2*q6/(4Pi ε0 )/m2 (y6[t]-y2[t]))/Sqrt[((x6[t]-x2[t])^2+(y6[t]-y2[t])^2+(z6[t]-z2[t])^2)^3]+
(-q2*q7/(4Pi ε0 )/m2 (y7[t]-y2[t]))/Sqrt[((x7[t]-x2[t])^2+(y7[t]-y2[t])^2+(z7[t]-z2[t])^2)^3]+
(-q2*q8/(4Pi ε0 )/m2 (y8[t]-y2[t]))/Sqrt[((x8[t]-x2[t])^2+(y8[t]-y2[t])^2+(z8[t]-z2[t])^2)^3]+
(-q2*q9/(4Pi ε0 )/m2 (y9[t]-y2[t]))/Sqrt[((x9[t]-x2[t])^2+(y9[t]-y2[t])^2+(z9[t]-z2[t])^2)^3]],
vz2'[t] ==
(G m1 (z1[t]-z2[t]))/Sqrt[((x2[t]-x1[t])^2+(y2[t]-y1[t])^2+(z2[t]-z1[t])^2)^3]+
(G m3 (z3[t]-z2[t]))/Sqrt[((x3[t]-x2[t])^2+(y3[t]-y2[t])^2+(z3[t]-z2[t])^2)^3]+
(G m4 (z4[t]-z2[t]))/Sqrt[((x4[t]-x2[t])^2+(y4[t]-y2[t])^2+(z4[t]-z2[t])^2)^3]+
(G m5 (z5[t]-z2[t]))/Sqrt[((x5[t]-x2[t])^2+(y5[t]-y2[t])^2+(z5[t]-z2[t])^2)^3]+
(G m6 (z6[t]-z2[t]))/Sqrt[((x6[t]-x2[t])^2+(y6[t]-y2[t])^2+(z6[t]-z2[t])^2)^3]+
(G m7 (z7[t]-z2[t]))/Sqrt[((x7[t]-x2[t])^2+(y7[t]-y2[t])^2+(z7[t]-z2[t])^2)^3]+
(G m8 (z8[t]-z2[t]))/Sqrt[((x8[t]-x2[t])^2+(y8[t]-y2[t])^2+(z8[t]-z2[t])^2)^3]+
(G m9 (z9[t]-z2[t]))/Sqrt[((x9[t]-x2[t])^2+(y9[t]-y2[t])^2+(z9[t]-z2[t])^2)^3]+
If[q2 == 0, 0,
(-q2*q1/(4Pi ε0 )/m2 (z1[t]-z2[t]))/Sqrt[((x2[t]-x1[t])^2+(y2[t]-y1[t])^2+(z2[t]-z1[t])^2)^3]+
(-q2*q3/(4Pi ε0 )/m2 (z3[t]-z2[t]))/Sqrt[((x3[t]-x2[t])^2+(y3[t]-y2[t])^2+(z3[t]-z2[t])^2)^3]+
(-q2*q4/(4Pi ε0 )/m2 (z4[t]-z2[t]))/Sqrt[((x4[t]-x2[t])^2+(y4[t]-y2[t])^2+(z4[t]-z2[t])^2)^3]+
(-q2*q5/(4Pi ε0 )/m2 (z5[t]-z2[t]))/Sqrt[((x5[t]-x2[t])^2+(y5[t]-y2[t])^2+(z5[t]-z2[t])^2)^3]+
(-q2*q6/(4Pi ε0 )/m2 (z6[t]-z2[t]))/Sqrt[((x6[t]-x2[t])^2+(y6[t]-y2[t])^2+(z6[t]-z2[t])^2)^3]+
(-q2*q7/(4Pi ε0 )/m2 (z7[t]-z2[t]))/Sqrt[((x7[t]-x2[t])^2+(y7[t]-y2[t])^2+(z7[t]-z2[t])^2)^3]+
(-q2*q8/(4Pi ε0 )/m2 (z8[t]-z2[t]))/Sqrt[((x8[t]-x2[t])^2+(y8[t]-y2[t])^2+(z8[t]-z2[t])^2)^3]+
(-q2*q9/(4Pi ε0 )/m2 (z9[t]-z2[t]))/Sqrt[((x9[t]-x2[t])^2+(y9[t]-y2[t])^2+(z9[t]-z2[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]+
(G m4 (x4[t]-x3[t]))/Sqrt[((x4[t]-x3[t])^2+(y4[t]-y3[t])^2+(z4[t]-z3[t])^2)^3]+
(G m5 (x5[t]-x3[t]))/Sqrt[((x5[t]-x3[t])^2+(y5[t]-y3[t])^2+(z5[t]-z3[t])^2)^3]+
(G m6 (x6[t]-x3[t]))/Sqrt[((x6[t]-x3[t])^2+(y6[t]-y3[t])^2+(z6[t]-z3[t])^2)^3]+
(G m7 (x7[t]-x3[t]))/Sqrt[((x7[t]-x3[t])^2+(y7[t]-y3[t])^2+(z7[t]-z3[t])^2)^3]+
(G m8 (x8[t]-x3[t]))/Sqrt[((x8[t]-x3[t])^2+(y8[t]-y3[t])^2+(z8[t]-z3[t])^2)^3]+
(G m9 (x9[t]-x3[t]))/Sqrt[((x9[t]-x3[t])^2+(y9[t]-y3[t])^2+(z9[t]-z3[t])^2)^3]+
If[q3 == 0, 0,
(-q3*q1/(4Pi ε0 )/m3 (x1[t]-x3[t]))/Sqrt[((x1[t]-x3[t])^2+(y1[t]-y3[t])^2+(z1[t]-z3[t])^2)^3]+
(-q3*q2/(4Pi ε0 )/m3 (x2[t]-x3[t]))/Sqrt[((x2[t]-x3[t])^2+(y2[t]-y3[t])^2+(z2[t]-z3[t])^2)^3]+
(-q3*q4/(4Pi ε0 )/m3 (x4[t]-x3[t]))/Sqrt[((x4[t]-x3[t])^2+(y4[t]-y3[t])^2+(z4[t]-z3[t])^2)^3]+
(-q3*q5/(4Pi ε0 )/m3 (x5[t]-x3[t]))/Sqrt[((x5[t]-x3[t])^2+(y5[t]-y3[t])^2+(z5[t]-z3[t])^2)^3]+
(-q3*q6/(4Pi ε0 )/m3 (x6[t]-x3[t]))/Sqrt[((x6[t]-x3[t])^2+(y6[t]-y3[t])^2+(z6[t]-z3[t])^2)^3]+
(-q3*q7/(4Pi ε0 )/m3 (x7[t]-x3[t]))/Sqrt[((x7[t]-x3[t])^2+(y7[t]-y3[t])^2+(z7[t]-z3[t])^2)^3]+
(-q3*q8/(4Pi ε0 )/m3 (x8[t]-x3[t]))/Sqrt[((x8[t]-x3[t])^2+(y8[t]-y3[t])^2+(z8[t]-z3[t])^2)^3]+
(-q3*q9/(4Pi ε0 )/m3 (x9[t]-x3[t]))/Sqrt[((x9[t]-x3[t])^2+(y9[t]-y3[t])^2+(z9[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]+
(G m4 (y4[t]-y3[t]))/Sqrt[((x4[t]-x3[t])^2+(y4[t]-y3[t])^2+(z4[t]-z3[t])^2)^3]+
(G m5 (y5[t]-y3[t]))/Sqrt[((x5[t]-x3[t])^2+(y5[t]-y3[t])^2+(z5[t]-z3[t])^2)^3]+
(G m6 (y6[t]-y3[t]))/Sqrt[((x6[t]-x3[t])^2+(y6[t]-y3[t])^2+(z6[t]-z3[t])^2)^3]+
(G m7 (y7[t]-y3[t]))/Sqrt[((x7[t]-x3[t])^2+(y7[t]-y3[t])^2+(z7[t]-z3[t])^2)^3]+
(G m8 (y8[t]-y3[t]))/Sqrt[((x8[t]-x3[t])^2+(y8[t]-y3[t])^2+(z8[t]-z3[t])^2)^3]+
(G m9 (y9[t]-y3[t]))/Sqrt[((x9[t]-x3[t])^2+(y9[t]-y3[t])^2+(z9[t]-z3[t])^2)^3]+
If[q3 == 0, 0,
(-q3*q1/(4Pi ε0 )/m3 (y1[t]-y3[t]))/Sqrt[((x1[t]-x3[t])^2+(y1[t]-y3[t])^2+(z1[t]-z3[t])^2)^3]+
(-q3*q2/(4Pi ε0 )/m3 (y2[t]-y3[t]))/Sqrt[((x2[t]-x3[t])^2+(y2[t]-y3[t])^2+(z2[t]-z3[t])^2)^3]+
(-q3*q4/(4Pi ε0 )/m3 (y4[t]-y3[t]))/Sqrt[((x4[t]-x3[t])^2+(y4[t]-y3[t])^2+(z4[t]-z3[t])^2)^3]+
(-q3*q5/(4Pi ε0 )/m3 (y5[t]-y3[t]))/Sqrt[((x5[t]-x3[t])^2+(y5[t]-y3[t])^2+(z5[t]-z3[t])^2)^3]+
(-q3*q6/(4Pi ε0 )/m3 (y6[t]-y3[t]))/Sqrt[((x6[t]-x3[t])^2+(y6[t]-y3[t])^2+(z6[t]-z3[t])^2)^3]+
(-q3*q7/(4Pi ε0 )/m3 (y7[t]-y3[t]))/Sqrt[((x7[t]-x3[t])^2+(y7[t]-y3[t])^2+(z7[t]-z3[t])^2)^3]+
(-q3*q8/(4Pi ε0 )/m3 (y8[t]-y3[t]))/Sqrt[((x8[t]-x3[t])^2+(y8[t]-y3[t])^2+(z8[t]-z3[t])^2)^3]+
(-q3*q9/(4Pi ε0 )/m3 (y9[t]-y3[t]))/Sqrt[((x9[t]-x3[t])^2+(y9[t]-y3[t])^2+(z9[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]+
(G m4 (z4[t]-z3[t]))/Sqrt[((x4[t]-x3[t])^2+(y4[t]-y3[t])^2+(z4[t]-z3[t])^2)^3]+
(G m5 (z5[t]-z3[t]))/Sqrt[((x5[t]-x3[t])^2+(y5[t]-y3[t])^2+(z5[t]-z3[t])^2)^3]+
(G m6 (z6[t]-z3[t]))/Sqrt[((x6[t]-x3[t])^2+(y6[t]-y3[t])^2+(z6[t]-z3[t])^2)^3]+
(G m7 (z7[t]-z3[t]))/Sqrt[((x7[t]-x3[t])^2+(y7[t]-y3[t])^2+(z7[t]-z3[t])^2)^3]+
(G m8 (z8[t]-z3[t]))/Sqrt[((x8[t]-x3[t])^2+(y8[t]-y3[t])^2+(z8[t]-z3[t])^2)^3]+
(G m9 (z9[t]-z3[t]))/Sqrt[((x9[t]-x3[t])^2+(y9[t]-y3[t])^2+(z9[t]-z3[t])^2)^3]+
If[q3 == 0, 0,
(-q3*q1/(4Pi ε0 )/m3 (z1[t]-z3[t]))/Sqrt[((x1[t]-x3[t])^2+(y1[t]-y3[t])^2+(z1[t]-z3[t])^2)^3]+
(-q3*q2/(4Pi ε0 )/m3 (z2[t]-z3[t]))/Sqrt[((x2[t]-x3[t])^2+(y2[t]-y3[t])^2+(z2[t]-z3[t])^2)^3]+
(-q3*q4/(4Pi ε0 )/m3 (z4[t]-z3[t]))/Sqrt[((x4[t]-x3[t])^2+(y4[t]-y3[t])^2+(z4[t]-z3[t])^2)^3]+
(-q3*q5/(4Pi ε0 )/m3 (z5[t]-z3[t]))/Sqrt[((x5[t]-x3[t])^2+(y5[t]-y3[t])^2+(z5[t]-z3[t])^2)^3]+
(-q3*q6/(4Pi ε0 )/m3 (z6[t]-z3[t]))/Sqrt[((x6[t]-x3[t])^2+(y6[t]-y3[t])^2+(z6[t]-z3[t])^2)^3]+
(-q3*q7/(4Pi ε0 )/m3 (z7[t]-z3[t]))/Sqrt[((x7[t]-x3[t])^2+(y7[t]-y3[t])^2+(z7[t]-z3[t])^2)^3]+
(-q3*q8/(4Pi ε0 )/m3 (z8[t]-z3[t]))/Sqrt[((x8[t]-x3[t])^2+(y8[t]-y3[t])^2+(z8[t]-z3[t])^2)^3]+
(-q3*q9/(4Pi ε0 )/m3 (z9[t]-z3[t]))/Sqrt[((x9[t]-x3[t])^2+(y9[t]-y3[t])^2+(z9[t]-z3[t])^2)^3]],
vx4'[t] ==
(G m1 (x1[t]-x4[t]))/Sqrt[((x1[t]-x4[t])^2+(y1[t]-y4[t])^2+(z1[t]-z4[t])^2)^3]+
(G m2 (x2[t]-x4[t]))/Sqrt[((x2[t]-x4[t])^2+(y2[t]-y4[t])^2+(z2[t]-z4[t])^2)^3]+
(G m3 (x3[t]-x4[t]))/Sqrt[((x3[t]-x4[t])^2+(y3[t]-y4[t])^2+(z3[t]-z4[t])^2)^3]+
(G m5 (x5[t]-x4[t]))/Sqrt[((x5[t]-x4[t])^2+(y5[t]-y4[t])^2+(z5[t]-z4[t])^2)^3]+
(G m6 (x6[t]-x4[t]))/Sqrt[((x6[t]-x4[t])^2+(y6[t]-y4[t])^2+(z6[t]-z4[t])^2)^3]+
(G m7 (x7[t]-x4[t]))/Sqrt[((x7[t]-x4[t])^2+(y7[t]-y4[t])^2+(z7[t]-z4[t])^2)^3]+
(G m8 (x8[t]-x4[t]))/Sqrt[((x8[t]-x4[t])^2+(y8[t]-y4[t])^2+(z8[t]-z4[t])^2)^3]+
(G m9 (x9[t]-x4[t]))/Sqrt[((x9[t]-x4[t])^2+(y9[t]-y4[t])^2+(z9[t]-z4[t])^2)^3]+
If[q4 == 0, 0,
(-q4*q1/(4Pi ε0 )/m4 (x1[t]-x4[t]))/Sqrt[((x1[t]-x4[t])^2+(y1[t]-y4[t])^2+(z1[t]-z4[t])^2)^3]+
(-q4*q2/(4Pi ε0 )/m4 (x2[t]-x4[t]))/Sqrt[((x2[t]-x4[t])^2+(y2[t]-y4[t])^2+(z2[t]-z4[t])^2)^3]+
(-q4*q3/(4Pi ε0 )/m4 (x3[t]-x4[t]))/Sqrt[((x3[t]-x4[t])^2+(y3[t]-y4[t])^2+(z3[t]-z4[t])^2)^3]+
(-q4*q5/(4Pi ε0 )/m4 (x5[t]-x4[t]))/Sqrt[((x5[t]-x4[t])^2+(y5[t]-y4[t])^2+(z5[t]-z4[t])^2)^3]+
(-q4*q6/(4Pi ε0 )/m4 (x6[t]-x4[t]))/Sqrt[((x6[t]-x4[t])^2+(y6[t]-y4[t])^2+(z6[t]-z4[t])^2)^3]+
(-q4*q7/(4Pi ε0 )/m4 (x7[t]-x4[t]))/Sqrt[((x7[t]-x4[t])^2+(y7[t]-y4[t])^2+(z7[t]-z4[t])^2)^3]+
(-q4*q8/(4Pi ε0 )/m4 (x8[t]-x4[t]))/Sqrt[((x8[t]-x4[t])^2+(y8[t]-y4[t])^2+(z8[t]-z4[t])^2)^3]+
(-q4*q9/(4Pi ε0 )/m4 (x9[t]-x4[t]))/Sqrt[((x9[t]-x4[t])^2+(y9[t]-y4[t])^2+(z9[t]-z4[t])^2)^3]],
vy4'[t] ==
(G m1 (y1[t]-y4[t]))/Sqrt[((x1[t]-x4[t])^2+(y1[t]-y4[t])^2+(z1[t]-z4[t])^2)^3]+
(G m2 (y2[t]-y4[t]))/Sqrt[((x2[t]-x4[t])^2+(y2[t]-y4[t])^2+(z2[t]-z4[t])^2)^3]+
(G m3 (y3[t]-y4[t]))/Sqrt[((x3[t]-x4[t])^2+(y3[t]-y4[t])^2+(z3[t]-z4[t])^2)^3]+
(G m5 (y5[t]-y4[t]))/Sqrt[((x5[t]-x4[t])^2+(y5[t]-y4[t])^2+(z5[t]-z4[t])^2)^3]+
(G m6 (y6[t]-y4[t]))/Sqrt[((x6[t]-x4[t])^2+(y6[t]-y4[t])^2+(z6[t]-z4[t])^2)^3]+
(G m7 (y7[t]-y4[t]))/Sqrt[((x7[t]-x4[t])^2+(y7[t]-y4[t])^2+(z7[t]-z4[t])^2)^3]+
(G m8 (y8[t]-y4[t]))/Sqrt[((x8[t]-x4[t])^2+(y8[t]-y4[t])^2+(z8[t]-z4[t])^2)^3]+
(G m9 (y9[t]-y4[t]))/Sqrt[((x9[t]-x4[t])^2+(y9[t]-y4[t])^2+(z9[t]-z4[t])^2)^3]+
If[q4 == 0, 0,
(-q4*q1/(4Pi ε0 )/m4 (y1[t]-y4[t]))/Sqrt[((x1[t]-x4[t])^2+(y1[t]-y4[t])^2+(z1[t]-z4[t])^2)^3]+
(-q4*q2/(4Pi ε0 )/m4 (y2[t]-y4[t]))/Sqrt[((x2[t]-x4[t])^2+(y2[t]-y4[t])^2+(z2[t]-z4[t])^2)^3]+
(-q4*q3/(4Pi ε0 )/m4 (y3[t]-y4[t]))/Sqrt[((x3[t]-x4[t])^2+(y3[t]-y4[t])^2+(z3[t]-z4[t])^2)^3]+
(-q4*q5/(4Pi ε0 )/m4 (y5[t]-y4[t]))/Sqrt[((x5[t]-x4[t])^2+(y5[t]-y4[t])^2+(z5[t]-z4[t])^2)^3]+
(-q4*q6/(4Pi ε0 )/m4 (y6[t]-y4[t]))/Sqrt[((x6[t]-x4[t])^2+(y6[t]-y4[t])^2+(z6[t]-z4[t])^2)^3]+
(-q4*q7/(4Pi ε0 )/m4 (y7[t]-y4[t]))/Sqrt[((x7[t]-x4[t])^2+(y7[t]-y4[t])^2+(z7[t]-z4[t])^2)^3]+
(-q4*q8/(4Pi ε0 )/m4 (y8[t]-y4[t]))/Sqrt[((x8[t]-x4[t])^2+(y8[t]-y4[t])^2+(z8[t]-z4[t])^2)^3]+
(-q4*q9/(4Pi ε0 )/m4 (y9[t]-y4[t]))/Sqrt[((x9[t]-x4[t])^2+(y9[t]-y4[t])^2+(z9[t]-z4[t])^2)^3]],
vz4'[t] ==
(G m1 (z1[t]-z4[t]))/Sqrt[((x1[t]-x4[t])^2+(y1[t]-y4[t])^2+(z1[t]-z4[t])^2)^3]+
(G m2 (z2[t]-z4[t]))/Sqrt[((x2[t]-x4[t])^2+(y2[t]-y4[t])^2+(z2[t]-z4[t])^2)^3]+
(G m3 (z3[t]-z4[t]))/Sqrt[((x3[t]-x4[t])^2+(y3[t]-y4[t])^2+(z3[t]-z4[t])^2)^3]+
(G m5 (z5[t]-z4[t]))/Sqrt[((x5[t]-x4[t])^2+(y5[t]-y4[t])^2+(z5[t]-z4[t])^2)^3]+
(G m6 (z6[t]-z4[t]))/Sqrt[((x6[t]-x4[t])^2+(y6[t]-y4[t])^2+(z6[t]-z4[t])^2)^3]+
(G m7 (z7[t]-z4[t]))/Sqrt[((x7[t]-x4[t])^2+(y7[t]-y4[t])^2+(z7[t]-z4[t])^2)^3]+
(G m8 (z8[t]-z4[t]))/Sqrt[((x8[t]-x4[t])^2+(y8[t]-y4[t])^2+(z8[t]-z4[t])^2)^3]+
(G m9 (z9[t]-z4[t]))/Sqrt[((x9[t]-x4[t])^2+(y9[t]-y4[t])^2+(z9[t]-z4[t])^2)^3]+
If[q4 == 0, 0,
(-q4*q1/(4Pi ε0 )/m4 (z1[t]-z4[t]))/Sqrt[((x1[t]-x4[t])^2+(y1[t]-y4[t])^2+(z1[t]-z4[t])^2)^3]+
(-q4*q2/(4Pi ε0 )/m4 (z2[t]-z4[t]))/Sqrt[((x2[t]-x4[t])^2+(y2[t]-y4[t])^2+(z2[t]-z4[t])^2)^3]+
(-q4*q3/(4Pi ε0 )/m4 (z3[t]-z4[t]))/Sqrt[((x3[t]-x4[t])^2+(y3[t]-y4[t])^2+(z3[t]-z4[t])^2)^3]+
(-q4*q5/(4Pi ε0 )/m4 (z5[t]-z4[t]))/Sqrt[((x5[t]-x4[t])^2+(y5[t]-y4[t])^2+(z5[t]-z4[t])^2)^3]+
(-q4*q6/(4Pi ε0 )/m4 (z6[t]-z4[t]))/Sqrt[((x6[t]-x4[t])^2+(y6[t]-y4[t])^2+(z6[t]-z4[t])^2)^3]+
(-q4*q7/(4Pi ε0 )/m4 (z7[t]-z4[t]))/Sqrt[((x7[t]-x4[t])^2+(y7[t]-y4[t])^2+(z7[t]-z4[t])^2)^3]+
(-q4*q8/(4Pi ε0 )/m4 (z8[t]-z4[t]))/Sqrt[((x8[t]-x4[t])^2+(y8[t]-y4[t])^2+(z8[t]-z4[t])^2)^3]+
(-q4*q9/(4Pi ε0 )/m4 (z9[t]-z4[t]))/Sqrt[((x9[t]-x4[t])^2+(y9[t]-y4[t])^2+(z9[t]-z4[t])^2)^3]],
vx5'[t] ==
(G m1 (x1[t]-x5[t]))/Sqrt[((x1[t]-x5[t])^2+(y1[t]-y5[t])^2+(z1[t]-z5[t])^2)^3]+
(G m2 (x2[t]-x5[t]))/Sqrt[((x2[t]-x5[t])^2+(y2[t]-y5[t])^2+(z2[t]-z5[t])^2)^3]+
(G m3 (x3[t]-x5[t]))/Sqrt[((x3[t]-x5[t])^2+(y3[t]-y5[t])^2+(z3[t]-z5[t])^2)^3]+
(G m4 (x4[t]-x5[t]))/Sqrt[((x4[t]-x5[t])^2+(y4[t]-y5[t])^2+(z4[t]-z5[t])^2)^3]+
(G m6 (x6[t]-x5[t]))/Sqrt[((x6[t]-x5[t])^2+(y6[t]-y5[t])^2+(z6[t]-z5[t])^2)^3]+
(G m7 (x7[t]-x5[t]))/Sqrt[((x7[t]-x5[t])^2+(y7[t]-y5[t])^2+(z7[t]-z5[t])^2)^3]+
(G m8 (x8[t]-x5[t]))/Sqrt[((x8[t]-x5[t])^2+(y8[t]-y5[t])^2+(z8[t]-z5[t])^2)^3]+
(G m9 (x9[t]-x5[t]))/Sqrt[((x9[t]-x5[t])^2+(y9[t]-y5[t])^2+(z9[t]-z5[t])^2)^3]+
If[q5 == 0, 0,
(-q5*q1/(4Pi ε0 )/m5 (x1[t]-x5[t]))/Sqrt[((x1[t]-x5[t])^2+(y1[t]-y5[t])^2+(z1[t]-z5[t])^2)^3]+
(-q5*q2/(4Pi ε0 )/m5 (x2[t]-x5[t]))/Sqrt[((x2[t]-x5[t])^2+(y2[t]-y5[t])^2+(z2[t]-z5[t])^2)^3]+
(-q5*q3/(4Pi ε0 )/m5 (x3[t]-x5[t]))/Sqrt[((x3[t]-x5[t])^2+(y3[t]-y5[t])^2+(z3[t]-z5[t])^2)^3]+
(-q5*q4/(4Pi ε0 )/m5 (x4[t]-x5[t]))/Sqrt[((x4[t]-x5[t])^2+(y4[t]-y5[t])^2+(z4[t]-z5[t])^2)^3]+
(-q5*q6/(4Pi ε0 )/m5 (x6[t]-x5[t]))/Sqrt[((x6[t]-x5[t])^2+(y6[t]-y5[t])^2+(z6[t]-z5[t])^2)^3]+
(-q5*q7/(4Pi ε0 )/m5 (x7[t]-x5[t]))/Sqrt[((x7[t]-x5[t])^2+(y7[t]-y5[t])^2+(z7[t]-z5[t])^2)^3]+
(-q5*q8/(4Pi ε0 )/m5 (x8[t]-x5[t]))/Sqrt[((x8[t]-x5[t])^2+(y8[t]-y5[t])^2+(z8[t]-z5[t])^2)^3]+
(-q5*q9/(4Pi ε0 )/m5 (x9[t]-x5[t]))/Sqrt[((x9[t]-x5[t])^2+(y9[t]-y5[t])^2+(z9[t]-z5[t])^2)^3]],
vy5'[t] ==
(G m1 (y1[t]-y5[t]))/Sqrt[((x1[t]-x5[t])^2+(y1[t]-y5[t])^2+(z1[t]-z5[t])^2)^3]+
(G m2 (y2[t]-y5[t]))/Sqrt[((x2[t]-x5[t])^2+(y2[t]-y5[t])^2+(z2[t]-z5[t])^2)^3]+
(G m3 (y3[t]-y5[t]))/Sqrt[((x3[t]-x5[t])^2+(y3[t]-y5[t])^2+(z3[t]-z5[t])^2)^3]+
(G m4 (y4[t]-y5[t]))/Sqrt[((x4[t]-x5[t])^2+(y4[t]-y5[t])^2+(z4[t]-z5[t])^2)^3]+
(G m6 (y6[t]-y5[t]))/Sqrt[((x6[t]-x5[t])^2+(y6[t]-y5[t])^2+(z6[t]-z5[t])^2)^3]+
(G m7 (y7[t]-y5[t]))/Sqrt[((x7[t]-x5[t])^2+(y7[t]-y5[t])^2+(z7[t]-z5[t])^2)^3]+
(G m8 (y8[t]-y5[t]))/Sqrt[((x8[t]-x5[t])^2+(y8[t]-y5[t])^2+(z8[t]-z5[t])^2)^3]+
(G m9 (y9[t]-y5[t]))/Sqrt[((x9[t]-x5[t])^2+(y9[t]-y5[t])^2+(z9[t]-z5[t])^2)^3]+
If[q5 == 0, 0,
(-q5*q1/(4Pi ε0 )/m5 (y1[t]-y5[t]))/Sqrt[((x1[t]-x5[t])^2+(y1[t]-y5[t])^2+(z1[t]-z5[t])^2)^3]+
(-q5*q2/(4Pi ε0 )/m5 (y2[t]-y5[t]))/Sqrt[((x2[t]-x5[t])^2+(y2[t]-y5[t])^2+(z2[t]-z5[t])^2)^3]+
(-q5*q3/(4Pi ε0 )/m5 (y3[t]-y5[t]))/Sqrt[((x3[t]-x5[t])^2+(y3[t]-y5[t])^2+(z3[t]-z5[t])^2)^3]+
(-q5*q4/(4Pi ε0 )/m5 (y4[t]-y5[t]))/Sqrt[((x4[t]-x5[t])^2+(y4[t]-y5[t])^2+(z4[t]-z5[t])^2)^3]+
(-q5*q6/(4Pi ε0 )/m5 (y6[t]-y5[t]))/Sqrt[((x6[t]-x5[t])^2+(y6[t]-y5[t])^2+(z6[t]-z5[t])^2)^3]+
(-q5*q7/(4Pi ε0 )/m5 (y7[t]-y5[t]))/Sqrt[((x7[t]-x5[t])^2+(y7[t]-y5[t])^2+(z7[t]-z5[t])^2)^3]+
(-q5*q8/(4Pi ε0 )/m5 (y8[t]-y5[t]))/Sqrt[((x8[t]-x5[t])^2+(y8[t]-y5[t])^2+(z8[t]-z5[t])^2)^3]+
(-q5*q9/(4Pi ε0 )/m5 (y9[t]-y5[t]))/Sqrt[((x9[t]-x5[t])^2+(y9[t]-y5[t])^2+(z9[t]-z5[t])^2)^3]],
vz5'[t] ==
(G m1 (z1[t]-z5[t]))/Sqrt[((x1[t]-x5[t])^2+(y1[t]-y5[t])^2+(z1[t]-z5[t])^2)^3]+
(G m2 (z2[t]-z5[t]))/Sqrt[((x2[t]-x5[t])^2+(y2[t]-y5[t])^2+(z2[t]-z5[t])^2)^3]+
(G m3 (z3[t]-z5[t]))/Sqrt[((x3[t]-x5[t])^2+(y3[t]-y5[t])^2+(z3[t]-z5[t])^2)^3]+
(G m4 (z4[t]-z5[t]))/Sqrt[((x4[t]-x5[t])^2+(y4[t]-y5[t])^2+(z4[t]-z5[t])^2)^3]+
(G m6 (z6[t]-z5[t]))/Sqrt[((x6[t]-x5[t])^2+(y6[t]-y5[t])^2+(z6[t]-z5[t])^2)^3]+
(G m7 (z7[t]-z5[t]))/Sqrt[((x7[t]-x5[t])^2+(y7[t]-y5[t])^2+(z7[t]-z5[t])^2)^3]+
(G m8 (z8[t]-z5[t]))/Sqrt[((x8[t]-x5[t])^2+(y8[t]-y5[t])^2+(z8[t]-z5[t])^2)^3]+
(G m9 (z9[t]-z5[t]))/Sqrt[((x9[t]-x5[t])^2+(y9[t]-y5[t])^2+(z9[t]-z5[t])^2)^3]+
If[q5 == 0, 0,
(-q5*q1/(4Pi ε0 )/m5 (z1[t]-z5[t]))/Sqrt[((x1[t]-x5[t])^2+(y1[t]-y5[t])^2+(z1[t]-z5[t])^2)^3]+
(-q5*q2/(4Pi ε0 )/m5 (z2[t]-z5[t]))/Sqrt[((x2[t]-x5[t])^2+(y2[t]-y5[t])^2+(z2[t]-z5[t])^2)^3]+
(-q5*q3/(4Pi ε0 )/m5 (z3[t]-z5[t]))/Sqrt[((x3[t]-x5[t])^2+(y3[t]-y5[t])^2+(z3[t]-z5[t])^2)^3]+
(-q5*q4/(4Pi ε0 )/m5 (z4[t]-z5[t]))/Sqrt[((x4[t]-x5[t])^2+(y4[t]-y5[t])^2+(z4[t]-z5[t])^2)^3]+
(-q5*q6/(4Pi ε0 )/m5 (z6[t]-z5[t]))/Sqrt[((x6[t]-x5[t])^2+(y6[t]-y5[t])^2+(z6[t]-z5[t])^2)^3]+
(-q5*q7/(4Pi ε0 )/m5 (z7[t]-z5[t]))/Sqrt[((x7[t]-x5[t])^2+(y7[t]-y5[t])^2+(z7[t]-z5[t])^2)^3]+
(-q5*q8/(4Pi ε0 )/m5 (z8[t]-z5[t]))/Sqrt[((x8[t]-x5[t])^2+(y8[t]-y5[t])^2+(z8[t]-z5[t])^2)^3]+
(-q5*q9/(4Pi ε0 )/m5 (z9[t]-z5[t]))/Sqrt[((x9[t]-x5[t])^2+(y9[t]-y5[t])^2+(z9[t]-z5[t])^2)^3]],
vx6'[t] ==
(G m1 (x1[t]-x6[t]))/Sqrt[((x1[t]-x6[t])^2+(y1[t]-y6[t])^2+(z1[t]-z6[t])^2)^3]+
(G m2 (x2[t]-x6[t]))/Sqrt[((x2[t]-x6[t])^2+(y2[t]-y6[t])^2+(z2[t]-z6[t])^2)^3]+
(G m3 (x3[t]-x6[t]))/Sqrt[((x3[t]-x6[t])^2+(y3[t]-y6[t])^2+(z3[t]-z6[t])^2)^3]+
(G m4 (x4[t]-x6[t]))/Sqrt[((x4[t]-x6[t])^2+(y4[t]-y6[t])^2+(z4[t]-z6[t])^2)^3]+
(G m5 (x5[t]-x6[t]))/Sqrt[((x5[t]-x6[t])^2+(y5[t]-y6[t])^2+(z5[t]-z6[t])^2)^3]+
(G m7 (x7[t]-x6[t]))/Sqrt[((x7[t]-x6[t])^2+(y7[t]-y6[t])^2+(z7[t]-z6[t])^2)^3]+
(G m8 (x8[t]-x6[t]))/Sqrt[((x8[t]-x6[t])^2+(y8[t]-y6[t])^2+(z8[t]-z6[t])^2)^3]+
(G m9 (x9[t]-x6[t]))/Sqrt[((x9[t]-x6[t])^2+(y9[t]-y6[t])^2+(z9[t]-z6[t])^2)^3]+
If[q6 == 0, 0,
(-q6*q1/(4Pi ε0 )/m6 (x1[t]-x6[t]))/Sqrt[((x1[t]-x6[t])^2+(y1[t]-y6[t])^2+(z1[t]-z6[t])^2)^3]+
(-q6*q2/(4Pi ε0 )/m6 (x2[t]-x6[t]))/Sqrt[((x2[t]-x6[t])^2+(y2[t]-y6[t])^2+(z2[t]-z6[t])^2)^3]+
(-q6*q3/(4Pi ε0 )/m6 (x3[t]-x6[t]))/Sqrt[((x3[t]-x6[t])^2+(y3[t]-y6[t])^2+(z3[t]-z6[t])^2)^3]+
(-q6*q4/(4Pi ε0 )/m6 (x4[t]-x6[t]))/Sqrt[((x4[t]-x6[t])^2+(y4[t]-y6[t])^2+(z4[t]-z6[t])^2)^3]+
(-q6*q5/(4Pi ε0 )/m6 (x5[t]-x6[t]))/Sqrt[((x5[t]-x6[t])^2+(y5[t]-y6[t])^2+(z5[t]-z6[t])^2)^3]+
(-q6*q7/(4Pi ε0 )/m6 (x7[t]-x6[t]))/Sqrt[((x7[t]-x6[t])^2+(y7[t]-y6[t])^2+(z7[t]-z6[t])^2)^3]+
(-q6*q8/(4Pi ε0 )/m6 (x8[t]-x6[t]))/Sqrt[((x8[t]-x6[t])^2+(y8[t]-y6[t])^2+(z8[t]-z6[t])^2)^3]+
(-q6*q9/(4Pi ε0 )/m6 (x9[t]-x6[t]))/Sqrt[((x9[t]-x6[t])^2+(y9[t]-y6[t])^2+(z9[t]-z6[t])^2)^3]],
vy6'[t] ==
(G m1 (y1[t]-y6[t]))/Sqrt[((x1[t]-x6[t])^2+(y1[t]-y6[t])^2+(z1[t]-z6[t])^2)^3]+
(G m2 (y2[t]-y6[t]))/Sqrt[((x2[t]-x6[t])^2+(y2[t]-y6[t])^2+(z2[t]-z6[t])^2)^3]+
(G m3 (y3[t]-y6[t]))/Sqrt[((x3[t]-x6[t])^2+(y3[t]-y6[t])^2+(z3[t]-z6[t])^2)^3]+
(G m4 (y4[t]-y6[t]))/Sqrt[((x4[t]-x6[t])^2+(y4[t]-y6[t])^2+(z4[t]-z6[t])^2)^3]+
(G m5 (y5[t]-y6[t]))/Sqrt[((x5[t]-x6[t])^2+(y5[t]-y6[t])^2+(z5[t]-z6[t])^2)^3]+
(G m7 (y7[t]-y6[t]))/Sqrt[((x7[t]-x6[t])^2+(y7[t]-y6[t])^2+(z7[t]-z6[t])^2)^3]+
(G m8 (y8[t]-y6[t]))/Sqrt[((x8[t]-x6[t])^2+(y8[t]-y6[t])^2+(z8[t]-z6[t])^2)^3]+
(G m9 (y9[t]-y6[t]))/Sqrt[((x9[t]-x6[t])^2+(y9[t]-y6[t])^2+(z9[t]-z6[t])^2)^3]+
If[q6 == 0, 0,
(-q6*q1/(4Pi ε0 )/m6 (y1[t]-y6[t]))/Sqrt[((x1[t]-x6[t])^2+(y1[t]-y6[t])^2+(z1[t]-z6[t])^2)^3]+
(-q6*q2/(4Pi ε0 )/m6 (y2[t]-y6[t]))/Sqrt[((x2[t]-x6[t])^2+(y2[t]-y6[t])^2+(z2[t]-z6[t])^2)^3]+
(-q6*q3/(4Pi ε0 )/m6 (y3[t]-y6[t]))/Sqrt[((x3[t]-x6[t])^2+(y3[t]-y6[t])^2+(z3[t]-z6[t])^2)^3]+
(-q6*q4/(4Pi ε0 )/m6 (y4[t]-y6[t]))/Sqrt[((x4[t]-x6[t])^2+(y4[t]-y6[t])^2+(z4[t]-z6[t])^2)^3]+
(-q6*q5/(4Pi ε0 )/m6 (y5[t]-y6[t]))/Sqrt[((x5[t]-x6[t])^2+(y5[t]-y6[t])^2+(z5[t]-z6[t])^2)^3]+
(-q6*q7/(4Pi ε0 )/m6 (y7[t]-y6[t]))/Sqrt[((x7[t]-x6[t])^2+(y7[t]-y6[t])^2+(z7[t]-z6[t])^2)^3]+
(-q6*q8/(4Pi ε0 )/m6 (y8[t]-y6[t]))/Sqrt[((x8[t]-x6[t])^2+(y8[t]-y6[t])^2+(z8[t]-z6[t])^2)^3]+
(-q6*q9/(4Pi ε0 )/m6 (y9[t]-y6[t]))/Sqrt[((x9[t]-x6[t])^2+(y9[t]-y6[t])^2+(z9[t]-z6[t])^2)^3]],
vz6'[t] ==
(G m1 (z1[t]-z6[t]))/Sqrt[((x1[t]-x6[t])^2+(y1[t]-y6[t])^2+(z1[t]-z6[t])^2)^3]+
(G m2 (z2[t]-z6[t]))/Sqrt[((x2[t]-x6[t])^2+(y2[t]-y6[t])^2+(z2[t]-z6[t])^2)^3]+
(G m3 (z3[t]-z6[t]))/Sqrt[((x3[t]-x6[t])^2+(y3[t]-y6[t])^2+(z3[t]-z6[t])^2)^3]+
(G m4 (z4[t]-z6[t]))/Sqrt[((x4[t]-x6[t])^2+(y4[t]-y6[t])^2+(z4[t]-z6[t])^2)^3]+
(G m5 (z5[t]-z6[t]))/Sqrt[((x5[t]-x6[t])^2+(y5[t]-y6[t])^2+(z5[t]-z6[t])^2)^3]+
(G m7 (z7[t]-z6[t]))/Sqrt[((x7[t]-x6[t])^2+(y7[t]-y6[t])^2+(z7[t]-z6[t])^2)^3]+
(G m8 (z8[t]-z6[t]))/Sqrt[((x8[t]-x6[t])^2+(y8[t]-y6[t])^2+(z8[t]-z6[t])^2)^3]+
(G m9 (z9[t]-z6[t]))/Sqrt[((x9[t]-x6[t])^2+(y9[t]-y6[t])^2+(z9[t]-z6[t])^2)^3]+
If[q6 == 0, 0,
(-q6*q1/(4Pi ε0 )/m6 (z1[t]-z6[t]))/Sqrt[((x1[t]-x6[t])^2+(y1[t]-y6[t])^2+(z1[t]-z6[t])^2)^3]+
(-q6*q2/(4Pi ε0 )/m6 (z2[t]-z6[t]))/Sqrt[((x2[t]-x6[t])^2+(y2[t]-y6[t])^2+(z2[t]-z6[t])^2)^3]+
(-q6*q3/(4Pi ε0 )/m6 (z3[t]-z6[t]))/Sqrt[((x3[t]-x6[t])^2+(y3[t]-y6[t])^2+(z3[t]-z6[t])^2)^3]+
(-q6*q4/(4Pi ε0 )/m6 (z4[t]-z6[t]))/Sqrt[((x4[t]-x6[t])^2+(y4[t]-y6[t])^2+(z4[t]-z6[t])^2)^3]+
(-q6*q5/(4Pi ε0 )/m6 (z5[t]-z6[t]))/Sqrt[((x5[t]-x6[t])^2+(y5[t]-y6[t])^2+(z5[t]-z6[t])^2)^3]+
(-q6*q7/(4Pi ε0 )/m6 (z7[t]-z6[t]))/Sqrt[((x7[t]-x6[t])^2+(y7[t]-y6[t])^2+(z7[t]-z6[t])^2)^3]+
(-q6*q8/(4Pi ε0 )/m6 (z8[t]-z6[t]))/Sqrt[((x8[t]-x6[t])^2+(y8[t]-y6[t])^2+(z8[t]-z6[t])^2)^3]+
(-q6*q9/(4Pi ε0 )/m6 (z9[t]-z6[t]))/Sqrt[((x9[t]-x6[t])^2+(y9[t]-y6[t])^2+(z9[t]-z6[t])^2)^3]],
vx7'[t] ==
(G m1 (x1[t]-x7[t]))/Sqrt[((x1[t]-x7[t])^2+(y1[t]-y7[t])^2+(z1[t]-z7[t])^2)^3]+
(G m2 (x2[t]-x7[t]))/Sqrt[((x2[t]-x7[t])^2+(y2[t]-y7[t])^2+(z2[t]-z7[t])^2)^3]+
(G m3 (x3[t]-x7[t]))/Sqrt[((x3[t]-x7[t])^2+(y3[t]-y7[t])^2+(z3[t]-z7[t])^2)^3]+
(G m4 (x4[t]-x7[t]))/Sqrt[((x4[t]-x7[t])^2+(y4[t]-y7[t])^2+(z4[t]-z7[t])^2)^3]+
(G m5 (x5[t]-x7[t]))/Sqrt[((x5[t]-x7[t])^2+(y5[t]-y7[t])^2+(z5[t]-z7[t])^2)^3]+
(G m6 (x6[t]-x7[t]))/Sqrt[((x6[t]-x7[t])^2+(y6[t]-y7[t])^2+(z6[t]-z7[t])^2)^3]+
(G m8 (x8[t]-x7[t]))/Sqrt[((x8[t]-x7[t])^2+(y8[t]-y7[t])^2+(z8[t]-z7[t])^2)^3]+
(G m9 (x9[t]-x7[t]))/Sqrt[((x9[t]-x7[t])^2+(y9[t]-y7[t])^2+(z9[t]-z7[t])^2)^3]+
If[q7 == 0, 0,
(-q7*q1/(4Pi ε0 )/m7 (x1[t]-x7[t]))/Sqrt[((x1[t]-x7[t])^2+(y1[t]-y7[t])^2+(z1[t]-z7[t])^2)^3]+
(-q7*q2/(4Pi ε0 )/m7 (x2[t]-x7[t]))/Sqrt[((x2[t]-x7[t])^2+(y2[t]-y7[t])^2+(z2[t]-z7[t])^2)^3]+
(-q7*q3/(4Pi ε0 )/m7 (x3[t]-x7[t]))/Sqrt[((x3[t]-x7[t])^2+(y3[t]-y7[t])^2+(z3[t]-z7[t])^2)^3]+
(-q7*q4/(4Pi ε0 )/m7 (x4[t]-x7[t]))/Sqrt[((x4[t]-x7[t])^2+(y4[t]-y7[t])^2+(z4[t]-z7[t])^2)^3]+
(-q7*q5/(4Pi ε0 )/m7 (x5[t]-x7[t]))/Sqrt[((x5[t]-x7[t])^2+(y5[t]-y7[t])^2+(z5[t]-z7[t])^2)^3]+
(-q7*q6/(4Pi ε0 )/m7 (x6[t]-x7[t]))/Sqrt[((x6[t]-x7[t])^2+(y6[t]-y7[t])^2+(z6[t]-z7[t])^2)^3]+
(-q7*q8/(4Pi ε0 )/m7 (x8[t]-x7[t]))/Sqrt[((x8[t]-x7[t])^2+(y8[t]-y7[t])^2+(z8[t]-z7[t])^2)^3]+
(-q7*q9/(4Pi ε0 )/m7 (x9[t]-x7[t]))/Sqrt[((x9[t]-x7[t])^2+(y9[t]-y7[t])^2+(z9[t]-z7[t])^2)^3]],
vy7'[t] ==
(G m1 (y1[t]-y7[t]))/Sqrt[((x1[t]-x7[t])^2+(y1[t]-y7[t])^2+(z1[t]-z7[t])^2)^3]+
(G m2 (y2[t]-y7[t]))/Sqrt[((x2[t]-x7[t])^2+(y2[t]-y7[t])^2+(z2[t]-z7[t])^2)^3]+
(G m3 (y3[t]-y7[t]))/Sqrt[((x3[t]-x7[t])^2+(y3[t]-y7[t])^2+(z3[t]-z7[t])^2)^3]+
(G m4 (y4[t]-y7[t]))/Sqrt[((x4[t]-x7[t])^2+(y4[t]-y7[t])^2+(z4[t]-z7[t])^2)^3]+
(G m5 (y5[t]-y7[t]))/Sqrt[((x5[t]-x7[t])^2+(y5[t]-y7[t])^2+(z5[t]-z7[t])^2)^3]+
(G m6 (y6[t]-y7[t]))/Sqrt[((x6[t]-x7[t])^2+(y6[t]-y7[t])^2+(z6[t]-z7[t])^2)^3]+
(G m8 (y8[t]-y7[t]))/Sqrt[((x8[t]-x7[t])^2+(y8[t]-y7[t])^2+(z8[t]-z7[t])^2)^3]+
(G m9 (y9[t]-y7[t]))/Sqrt[((x9[t]-x7[t])^2+(y9[t]-y7[t])^2+(z9[t]-z7[t])^2)^3]+
If[q7 == 0, 0,
(-q7*q1/(4Pi ε0 )/m7 (y1[t]-y7[t]))/Sqrt[((x1[t]-x7[t])^2+(y1[t]-y7[t])^2+(z1[t]-z7[t])^2)^3]+
(-q7*q2/(4Pi ε0 )/m7 (y2[t]-y7[t]))/Sqrt[((x2[t]-x7[t])^2+(y2[t]-y7[t])^2+(z2[t]-z7[t])^2)^3]+
(-q7*q3/(4Pi ε0 )/m7 (y3[t]-y7[t]))/Sqrt[((x3[t]-x7[t])^2+(y3[t]-y7[t])^2+(z3[t]-z7[t])^2)^3]+
(-q7*q4/(4Pi ε0 )/m7 (y4[t]-y7[t]))/Sqrt[((x4[t]-x7[t])^2+(y4[t]-y7[t])^2+(z4[t]-z7[t])^2)^3]+
(-q7*q5/(4Pi ε0 )/m7 (y5[t]-y7[t]))/Sqrt[((x5[t]-x7[t])^2+(y5[t]-y7[t])^2+(z5[t]-z7[t])^2)^3]+
(-q7*q6/(4Pi ε0 )/m7 (y6[t]-y7[t]))/Sqrt[((x6[t]-x7[t])^2+(y6[t]-y7[t])^2+(z6[t]-z7[t])^2)^3]+
(-q7*q8/(4Pi ε0 )/m7 (y8[t]-y7[t]))/Sqrt[((x8[t]-x7[t])^2+(y8[t]-y7[t])^2+(z8[t]-z7[t])^2)^3]+
(-q7*q9/(4Pi ε0 )/m7 (y9[t]-y7[t]))/Sqrt[((x9[t]-x7[t])^2+(y9[t]-y7[t])^2+(z9[t]-z7[t])^2)^3]],
vz7'[t] ==
(G m1 (z1[t]-z7[t]))/Sqrt[((x1[t]-x7[t])^2+(y1[t]-y7[t])^2+(z1[t]-z7[t])^2)^3]+
(G m2 (z2[t]-z7[t]))/Sqrt[((x2[t]-x7[t])^2+(y2[t]-y7[t])^2+(z2[t]-z7[t])^2)^3]+
(G m3 (z3[t]-z7[t]))/Sqrt[((x3[t]-x7[t])^2+(y3[t]-y7[t])^2+(z3[t]-z7[t])^2)^3]+
(G m4 (z4[t]-z7[t]))/Sqrt[((x4[t]-x7[t])^2+(y4[t]-y7[t])^2+(z4[t]-z7[t])^2)^3]+
(G m5 (z5[t]-z7[t]))/Sqrt[((x5[t]-x7[t])^2+(y5[t]-y7[t])^2+(z5[t]-z7[t])^2)^3]+
(G m6 (z6[t]-z7[t]))/Sqrt[((x6[t]-x7[t])^2+(y6[t]-y7[t])^2+(z6[t]-z7[t])^2)^3]+
(G m8 (z8[t]-z7[t]))/Sqrt[((x8[t]-x7[t])^2+(y8[t]-y7[t])^2+(z8[t]-z7[t])^2)^3]+
(G m9 (z9[t]-z7[t]))/Sqrt[((x9[t]-x7[t])^2+(y9[t]-y7[t])^2+(z9[t]-z7[t])^2)^3]+
If[q7 == 0, 0,
(-q7*q1/(4Pi ε0 )/m7 (z1[t]-z7[t]))/Sqrt[((x1[t]-x7[t])^2+(y1[t]-y7[t])^2+(z1[t]-z7[t])^2)^3]+
(-q7*q2/(4Pi ε0 )/m7 (z2[t]-z7[t]))/Sqrt[((x2[t]-x7[t])^2+(y2[t]-y7[t])^2+(z2[t]-z7[t])^2)^3]+
(-q7*q3/(4Pi ε0 )/m7 (z3[t]-z7[t]))/Sqrt[((x3[t]-x7[t])^2+(y3[t]-y7[t])^2+(z3[t]-z7[t])^2)^3]+
(-q7*q4/(4Pi ε0 )/m7 (z4[t]-z7[t]))/Sqrt[((x4[t]-x7[t])^2+(y4[t]-y7[t])^2+(z4[t]-z7[t])^2)^3]+
(-q7*q5/(4Pi ε0 )/m7 (z5[t]-z7[t]))/Sqrt[((x5[t]-x7[t])^2+(y5[t]-y7[t])^2+(z5[t]-z7[t])^2)^3]+
(-q7*q6/(4Pi ε0 )/m7 (z6[t]-z7[t]))/Sqrt[((x6[t]-x7[t])^2+(y6[t]-y7[t])^2+(z6[t]-z7[t])^2)^3]+
(-q7*q8/(4Pi ε0 )/m7 (z8[t]-z7[t]))/Sqrt[((x8[t]-x7[t])^2+(y8[t]-y7[t])^2+(z8[t]-z7[t])^2)^3]+
(-q7*q9/(4Pi ε0 )/m7 (z9[t]-z7[t]))/Sqrt[((x9[t]-x7[t])^2+(y9[t]-y7[t])^2+(z9[t]-z7[t])^2)^3]],
vx8'[t] ==
(G m1 (x1[t]-x8[t]))/Sqrt[((x1[t]-x8[t])^2+(y1[t]-y8[t])^2+(z1[t]-z8[t])^2)^3]+
(G m2 (x2[t]-x8[t]))/Sqrt[((x2[t]-x8[t])^2+(y2[t]-y8[t])^2+(z2[t]-z8[t])^2)^3]+
(G m3 (x3[t]-x8[t]))/Sqrt[((x3[t]-x8[t])^2+(y3[t]-y8[t])^2+(z3[t]-z8[t])^2)^3]+
(G m4 (x4[t]-x8[t]))/Sqrt[((x4[t]-x8[t])^2+(y4[t]-y8[t])^2+(z4[t]-z8[t])^2)^3]+
(G m5 (x5[t]-x8[t]))/Sqrt[((x5[t]-x8[t])^2+(y5[t]-y8[t])^2+(z5[t]-z8[t])^2)^3]+
(G m6 (x6[t]-x8[t]))/Sqrt[((x6[t]-x8[t])^2+(y6[t]-y8[t])^2+(z6[t]-z8[t])^2)^3]+
(G m7 (x7[t]-x8[t]))/Sqrt[((x7[t]-x8[t])^2+(y7[t]-y8[t])^2+(z7[t]-z8[t])^2)^3]+
(G m9 (x9[t]-x8[t]))/Sqrt[((x9[t]-x8[t])^2+(y9[t]-y8[t])^2+(z9[t]-z8[t])^2)^3]+
If[q8 == 0, 0,
(-q8*q1/(4Pi ε0 )/m8 (x1[t]-x8[t]))/Sqrt[((x1[t]-x8[t])^2+(y1[t]-y8[t])^2+(z1[t]-z8[t])^2)^3]+
(-q8*q2/(4Pi ε0 )/m8 (x2[t]-x8[t]))/Sqrt[((x2[t]-x8[t])^2+(y2[t]-y8[t])^2+(z2[t]-z8[t])^2)^3]+
(-q8*q3/(4Pi ε0 )/m8 (x3[t]-x8[t]))/Sqrt[((x3[t]-x8[t])^2+(y3[t]-y8[t])^2+(z3[t]-z8[t])^2)^3]+
(-q8*q4/(4Pi ε0 )/m8 (x4[t]-x8[t]))/Sqrt[((x4[t]-x8[t])^2+(y4[t]-y8[t])^2+(z4[t]-z8[t])^2)^3]+
(-q8*q5/(4Pi ε0 )/m8 (x5[t]-x8[t]))/Sqrt[((x5[t]-x8[t])^2+(y5[t]-y8[t])^2+(z5[t]-z8[t])^2)^3]+
(-q8*q6/(4Pi ε0 )/m8 (x6[t]-x8[t]))/Sqrt[((x6[t]-x8[t])^2+(y6[t]-y8[t])^2+(z6[t]-z8[t])^2)^3]+
(-q8*q7/(4Pi ε0 )/m8 (x7[t]-x8[t]))/Sqrt[((x7[t]-x8[t])^2+(y7[t]-y8[t])^2+(z7[t]-z8[t])^2)^3]+
(-q8*q9/(4Pi ε0 )/m8 (x9[t]-x8[t]))/Sqrt[((x9[t]-x8[t])^2+(y9[t]-y8[t])^2+(z9[t]-z8[t])^2)^3]],
vy8'[t] ==
(G m1 (y1[t]-y8[t]))/Sqrt[((x1[t]-x8[t])^2+(y1[t]-y8[t])^2+(z1[t]-z8[t])^2)^3]+
(G m2 (y2[t]-y8[t]))/Sqrt[((x2[t]-x8[t])^2+(y2[t]-y8[t])^2+(z2[t]-z8[t])^2)^3]+
(G m3 (y3[t]-y8[t]))/Sqrt[((x3[t]-x8[t])^2+(y3[t]-y8[t])^2+(z3[t]-z8[t])^2)^3]+
(G m4 (y4[t]-y8[t]))/Sqrt[((x4[t]-x8[t])^2+(y4[t]-y8[t])^2+(z4[t]-z8[t])^2)^3]+
(G m5 (y5[t]-y8[t]))/Sqrt[((x5[t]-x8[t])^2+(y5[t]-y8[t])^2+(z5[t]-z8[t])^2)^3]+
(G m6 (y6[t]-y8[t]))/Sqrt[((x6[t]-x8[t])^2+(y6[t]-y8[t])^2+(z6[t]-z8[t])^2)^3]+
(G m7 (y7[t]-y8[t]))/Sqrt[((x7[t]-x8[t])^2+(y7[t]-y8[t])^2+(z7[t]-z8[t])^2)^3]+
(G m9 (y9[t]-y8[t]))/Sqrt[((x9[t]-x8[t])^2+(y9[t]-y8[t])^2+(z9[t]-z8[t])^2)^3]+
If[q8 == 0, 0,
(-q8*q1/(4Pi ε0 )/m8 (y1[t]-y8[t]))/Sqrt[((x1[t]-x8[t])^2+(y1[t]-y8[t])^2+(z1[t]-z8[t])^2)^3]+
(-q8*q2/(4Pi ε0 )/m8 (y2[t]-y8[t]))/Sqrt[((x2[t]-x8[t])^2+(y2[t]-y8[t])^2+(z2[t]-z8[t])^2)^3]+
(-q8*q3/(4Pi ε0 )/m8 (y3[t]-y8[t]))/Sqrt[((x3[t]-x8[t])^2+(y3[t]-y8[t])^2+(z3[t]-z8[t])^2)^3]+
(-q8*q4/(4Pi ε0 )/m8 (y4[t]-y8[t]))/Sqrt[((x4[t]-x8[t])^2+(y4[t]-y8[t])^2+(z4[t]-z8[t])^2)^3]+
(-q8*q5/(4Pi ε0 )/m8 (y5[t]-y8[t]))/Sqrt[((x5[t]-x8[t])^2+(y5[t]-y8[t])^2+(z5[t]-z8[t])^2)^3]+
(-q8*q6/(4Pi ε0 )/m8 (y6[t]-y8[t]))/Sqrt[((x6[t]-x8[t])^2+(y6[t]-y8[t])^2+(z6[t]-z8[t])^2)^3]+
(-q8*q7/(4Pi ε0 )/m8 (y7[t]-y8[t]))/Sqrt[((x7[t]-x8[t])^2+(y7[t]-y8[t])^2+(z7[t]-z8[t])^2)^3]+
(-q8*q9/(4Pi ε0 )/m8 (y9[t]-y8[t]))/Sqrt[((x9[t]-x8[t])^2+(y9[t]-y8[t])^2+(z9[t]-z8[t])^2)^3]],
vz8'[t] ==
(G m1 (z1[t]-z8[t]))/Sqrt[((x1[t]-x8[t])^2+(y1[t]-y8[t])^2+(z1[t]-z8[t])^2)^3]+
(G m2 (z2[t]-z8[t]))/Sqrt[((x2[t]-x8[t])^2+(y2[t]-y8[t])^2+(z2[t]-z8[t])^2)^3]+
(G m3 (z3[t]-z8[t]))/Sqrt[((x3[t]-x8[t])^2+(y3[t]-y8[t])^2+(z3[t]-z8[t])^2)^3]+
(G m4 (z4[t]-z8[t]))/Sqrt[((x4[t]-x8[t])^2+(y4[t]-y8[t])^2+(z4[t]-z8[t])^2)^3]+
(G m5 (z5[t]-z8[t]))/Sqrt[((x5[t]-x8[t])^2+(y5[t]-y8[t])^2+(z5[t]-z8[t])^2)^3]+
(G m6 (z6[t]-z8[t]))/Sqrt[((x6[t]-x8[t])^2+(y6[t]-y8[t])^2+(z6[t]-z8[t])^2)^3]+
(G m7 (z7[t]-z8[t]))/Sqrt[((x7[t]-x8[t])^2+(y7[t]-y8[t])^2+(z7[t]-z8[t])^2)^3]+
(G m9 (z9[t]-z8[t]))/Sqrt[((x9[t]-x8[t])^2+(y9[t]-y8[t])^2+(z9[t]-z8[t])^2)^3]+
If[q8 == 0, 0,
(-q8*q1/(4Pi ε0 )/m8 (z1[t]-z8[t]))/Sqrt[((x1[t]-x8[t])^2+(y1[t]-y8[t])^2+(z1[t]-z8[t])^2)^3]+
(-q8*q2/(4Pi ε0 )/m8 (z2[t]-z8[t]))/Sqrt[((x2[t]-x8[t])^2+(y2[t]-y8[t])^2+(z2[t]-z8[t])^2)^3]+
(-q8*q3/(4Pi ε0 )/m8 (z3[t]-z8[t]))/Sqrt[((x3[t]-x8[t])^2+(y3[t]-y8[t])^2+(z3[t]-z8[t])^2)^3]+
(-q8*q4/(4Pi ε0 )/m8 (z4[t]-z8[t]))/Sqrt[((x4[t]-x8[t])^2+(y4[t]-y8[t])^2+(z4[t]-z8[t])^2)^3]+
(-q8*q5/(4Pi ε0 )/m8 (z5[t]-z8[t]))/Sqrt[((x5[t]-x8[t])^2+(y5[t]-y8[t])^2+(z5[t]-z8[t])^2)^3]+
(-q8*q6/(4Pi ε0 )/m8 (z6[t]-z8[t]))/Sqrt[((x6[t]-x8[t])^2+(y6[t]-y8[t])^2+(z6[t]-z8[t])^2)^3]+
(-q8*q7/(4Pi ε0 )/m8 (z7[t]-z8[t]))/Sqrt[((x7[t]-x8[t])^2+(y7[t]-y8[t])^2+(z7[t]-z8[t])^2)^3]+
(-q8*q9/(4Pi ε0 )/m8 (z9[t]-z8[t]))/Sqrt[((x9[t]-x8[t])^2+(y9[t]-y8[t])^2+(z9[t]-z8[t])^2)^3]],
vx9'[t] ==
(G m1 (x1[t]-x9[t]))/Sqrt[((x1[t]-x9[t])^2+(y1[t]-y9[t])^2+(z1[t]-z9[t])^2)^3]+
(G m2 (x2[t]-x9[t]))/Sqrt[((x2[t]-x9[t])^2+(y2[t]-y9[t])^2+(z2[t]-z9[t])^2)^3]+
(G m3 (x3[t]-x9[t]))/Sqrt[((x3[t]-x9[t])^2+(y3[t]-y9[t])^2+(z3[t]-z9[t])^2)^3]+
(G m4 (x4[t]-x9[t]))/Sqrt[((x4[t]-x9[t])^2+(y4[t]-y9[t])^2+(z4[t]-z9[t])^2)^3]+
(G m5 (x5[t]-x9[t]))/Sqrt[((x5[t]-x9[t])^2+(y5[t]-y9[t])^2+(z5[t]-z9[t])^2)^3]+
(G m6 (x6[t]-x9[t]))/Sqrt[((x6[t]-x9[t])^2+(y6[t]-y9[t])^2+(z6[t]-z9[t])^2)^3]+
(G m7 (x7[t]-x9[t]))/Sqrt[((x7[t]-x9[t])^2+(y7[t]-y9[t])^2+(z7[t]-z9[t])^2)^3]+
(G m8 (x8[t]-x9[t]))/Sqrt[((x8[t]-x9[t])^2+(y8[t]-y9[t])^2+(z8[t]-z9[t])^2)^3]+
If[q9 == 0, 0,
(-q9*q1/(4Pi ε0 )/m9 (x1[t]-x9[t]))/Sqrt[((x1[t]-x9[t])^2+(y1[t]-y9[t])^2+(z1[t]-z9[t])^2)^3]+
(-q9*q2/(4Pi ε0 )/m9 (x2[t]-x9[t]))/Sqrt[((x2[t]-x9[t])^2+(y2[t]-y9[t])^2+(z2[t]-z9[t])^2)^3]+
(-q9*q3/(4Pi ε0 )/m9 (x3[t]-x9[t]))/Sqrt[((x3[t]-x9[t])^2+(y3[t]-y9[t])^2+(z3[t]-z9[t])^2)^3]+
(-q9*q4/(4Pi ε0 )/m9 (x4[t]-x9[t]))/Sqrt[((x4[t]-x9[t])^2+(y4[t]-y9[t])^2+(z4[t]-z9[t])^2)^3]+
(-q9*q5/(4Pi ε0 )/m9 (x5[t]-x9[t]))/Sqrt[((x5[t]-x9[t])^2+(y5[t]-y9[t])^2+(z5[t]-z9[t])^2)^3]+
(-q9*q6/(4Pi ε0 )/m9 (x6[t]-x9[t]))/Sqrt[((x6[t]-x9[t])^2+(y6[t]-y9[t])^2+(z6[t]-z9[t])^2)^3]+
(-q9*q7/(4Pi ε0 )/m9 (x7[t]-x9[t]))/Sqrt[((x7[t]-x9[t])^2+(y7[t]-y9[t])^2+(z7[t]-z9[t])^2)^3]+
(-q9*q8/(4Pi ε0 )/m9 (x8[t]-x9[t]))/Sqrt[((x8[t]-x9[t])^2+(y8[t]-y9[t])^2+(z8[t]-z9[t])^2)^3]],
vy9'[t] ==
(G m1 (y1[t]-y9[t]))/Sqrt[((x1[t]-x9[t])^2+(y1[t]-y9[t])^2+(z1[t]-z9[t])^2)^3]+
(G m2 (y2[t]-y9[t]))/Sqrt[((x2[t]-x9[t])^2+(y2[t]-y9[t])^2+(z2[t]-z9[t])^2)^3]+
(G m3 (y3[t]-y9[t]))/Sqrt[((x3[t]-x9[t])^2+(y3[t]-y9[t])^2+(z3[t]-z9[t])^2)^3]+
(G m4 (y4[t]-y9[t]))/Sqrt[((x4[t]-x9[t])^2+(y4[t]-y9[t])^2+(z4[t]-z9[t])^2)^3]+
(G m5 (y5[t]-y9[t]))/Sqrt[((x5[t]-x9[t])^2+(y5[t]-y9[t])^2+(z5[t]-z9[t])^2)^3]+
(G m6 (y6[t]-y9[t]))/Sqrt[((x6[t]-x9[t])^2+(y6[t]-y9[t])^2+(z6[t]-z9[t])^2)^3]+
(G m7 (y7[t]-y9[t]))/Sqrt[((x7[t]-x9[t])^2+(y7[t]-y9[t])^2+(z7[t]-z9[t])^2)^3]+
(G m8 (y8[t]-y9[t]))/Sqrt[((x8[t]-x9[t])^2+(y8[t]-y9[t])^2+(z8[t]-z9[t])^2)^3]+
If[q9 == 0, 0,
(-q9*q1/(4Pi ε0 )/m9 (y1[t]-y9[t]))/Sqrt[((x1[t]-x9[t])^2+(y1[t]-y9[t])^2+(z1[t]-z9[t])^2)^3]+
(-q9*q2/(4Pi ε0 )/m9 (y2[t]-y9[t]))/Sqrt[((x2[t]-x9[t])^2+(y2[t]-y9[t])^2+(z2[t]-z9[t])^2)^3]+
(-q9*q3/(4Pi ε0 )/m9 (y3[t]-y9[t]))/Sqrt[((x3[t]-x9[t])^2+(y3[t]-y9[t])^2+(z3[t]-z9[t])^2)^3]+
(-q9*q4/(4Pi ε0 )/m9 (y4[t]-y9[t]))/Sqrt[((x4[t]-x9[t])^2+(y4[t]-y9[t])^2+(z4[t]-z9[t])^2)^3]+
(-q9*q5/(4Pi ε0 )/m9 (y5[t]-y9[t]))/Sqrt[((x5[t]-x9[t])^2+(y5[t]-y9[t])^2+(z5[t]-z9[t])^2)^3]+
(-q9*q6/(4Pi ε0 )/m9 (y6[t]-y9[t]))/Sqrt[((x6[t]-x9[t])^2+(y6[t]-y9[t])^2+(z6[t]-z9[t])^2)^3]+
(-q9*q7/(4Pi ε0 )/m9 (y7[t]-y9[t]))/Sqrt[((x7[t]-x9[t])^2+(y7[t]-y9[t])^2+(z7[t]-z9[t])^2)^3]+
(-q9*q8/(4Pi ε0 )/m9 (y8[t]-y9[t]))/Sqrt[((x8[t]-x9[t])^2+(y8[t]-y9[t])^2+(z8[t]-z9[t])^2)^3]],
vz9'[t] ==
(G m1 (z1[t]-z9[t]))/Sqrt[((x1[t]-x9[t])^2+(y1[t]-y9[t])^2+(z1[t]-z9[t])^2)^3]+
(G m2 (z2[t]-z9[t]))/Sqrt[((x2[t]-x9[t])^2+(y2[t]-y9[t])^2+(z2[t]-z9[t])^2)^3]+
(G m3 (z3[t]-z9[t]))/Sqrt[((x3[t]-x9[t])^2+(y3[t]-y9[t])^2+(z3[t]-z9[t])^2)^3]+
(G m4 (z4[t]-z9[t]))/Sqrt[((x4[t]-x9[t])^2+(y4[t]-y9[t])^2+(z4[t]-z9[t])^2)^3]+
(G m5 (z5[t]-z9[t]))/Sqrt[((x5[t]-x9[t])^2+(y5[t]-y9[t])^2+(z5[t]-z9[t])^2)^3]+
(G m6 (z6[t]-z9[t]))/Sqrt[((x6[t]-x9[t])^2+(y6[t]-y9[t])^2+(z6[t]-z9[t])^2)^3]+
(G m7 (z7[t]-z9[t]))/Sqrt[((x7[t]-x9[t])^2+(y7[t]-y9[t])^2+(z7[t]-z9[t])^2)^3]+
(G m8 (z8[t]-z9[t]))/Sqrt[((x8[t]-x9[t])^2+(y8[t]-y9[t])^2+(z8[t]-z9[t])^2)^3]+
If[q9 == 0, 0,
(-q9*q1/(4Pi ε0 )/m9 (z1[t]-z9[t]))/Sqrt[((x1[t]-x9[t])^2+(y1[t]-y9[t])^2+(z1[t]-z9[t])^2)^3]+
(-q9*q2/(4Pi ε0 )/m9 (z2[t]-z9[t]))/Sqrt[((x2[t]-x9[t])^2+(y2[t]-y9[t])^2+(z2[t]-z9[t])^2)^3]+
(-q9*q3/(4Pi ε0 )/m9 (z3[t]-z9[t]))/Sqrt[((x3[t]-x9[t])^2+(y3[t]-y9[t])^2+(z3[t]-z9[t])^2)^3]+
(-q9*q4/(4Pi ε0 )/m9 (z4[t]-z9[t]))/Sqrt[((x4[t]-x9[t])^2+(y4[t]-y9[t])^2+(z4[t]-z9[t])^2)^3]+
(-q9*q5/(4Pi ε0 )/m9 (z5[t]-z9[t]))/Sqrt[((x5[t]-x9[t])^2+(y5[t]-y9[t])^2+(z5[t]-z9[t])^2)^3]+
(-q9*q6/(4Pi ε0 )/m9 (z6[t]-z9[t]))/Sqrt[((x6[t]-x9[t])^2+(y6[t]-y9[t])^2+(z6[t]-z9[t])^2)^3]+
(-q9*q7/(4Pi ε0 )/m9 (z7[t]-z9[t]))/Sqrt[((x7[t]-x9[t])^2+(y7[t]-y9[t])^2+(z7[t]-z9[t])^2)^3]+
(-q9*q8/(4Pi ε0 )/m9 (z8[t]-z9[t]))/Sqrt[((x8[t]-x9[t])^2+(y8[t]-y9[t])^2+(z8[t]-z9[t])^2)^3]],
x1[0] == x1x, y1[0] == y1y, z1[0] == z1z,
x2[0] == x2x, y2[0] == y2y, z2[0] == z2z,
x3[0] == x3x, y3[0] == y3y, z3[0] == z3z,
x4[0] == x4x, y4[0] == y4y, z4[0] == z4z,
x5[0] == x5x, y5[0] == y5y, z5[0] == z5z,
x6[0] == x6x, y6[0] == y6y, z6[0] == z6z,
x7[0] == x7x, y7[0] == y7y, z7[0] == z7z,
x8[0] == x8x, y8[0] == y8y, z8[0] == z8z,
x9[0] == x9x, y9[0] == y9y, z9[0] == z9z,
vx1[0] == v1x, vy1[0] == v1y, vz1[0] == v1z,
vx2[0] == v2x, vy2[0] == v2y, vz2[0] == v2z,
vx3[0] == v3x, vy3[0] == v3y, vz3[0] == v3z,
vx4[0] == v4x, vy4[0] == v4y, vz4[0] == v4z,
vx5[0] == v5x, vy5[0] == v5y, vz5[0] == v5z,
vx6[0] == v6x, vy6[0] == v6y, vz6[0] == v6z,
vx7[0] == v7x, vy7[0] == v7y, vz7[0] == v7z,
vx8[0] == v8x, vy8[0] == v8y, vz8[0] == v8z,
vx9[0] == v9x, vy9[0] == v9y, vz9[0] == v9z},
{x1, x2, x3, x4, x5, x6, x7, x8, x9, y1, y2, y3, y4, y5, y6, y7, y8, y9, z1, z2, z3, z4, z5, z6, z7, z8, z9,
vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vx9, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vy9, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, vz9},
{t, 0, Tmax},
WorkingPrecision-> wp,
MaxSteps-> Infinity,
Method-> mta,
InterpolationOrder-> All,
StepMonitor :> (laststep=plunge; plunge=t;
stepsize=plunge-laststep;), Method->{"EventLocator",
"Event" :> (If[stepsize<1*^-4, 0, 1])}];
(* Position, Geschwindigkeit *)
f2p[t_]={{x1[t], y1[t], z1[t]}, {x2[t], y2[t], z2[t]}, {x3[t], y3[t], z3[t]}, {x4[t], y4[t], z4[t]}, {x5[t], y5[t], z5[t]}, {x6[t], y6[t], z6[t]}, {x7[t], y7[t], z7[t]}, {x8[t], y8[t], z8[t]}, {x9[t], y9[t], z9[t]}}/.nds[[1]];
f2v[t_]={{vx1[t], vy1[t], vz1[t]}, {vx2[t], vy2[t], vz2[t]}, {vx3[t], vy3[t], vz3[t]}, {vx4[t], vy4[t], vz4[t]}, {vx5[t], vy5[t], vz5[t]}, {vx6[t], vy6[t], vz6[t]}, {vx7[t], vy7[t], vz7[t]}, {vx8[t], vy8[t], vz8[t]}, {vx9[t], vy9[t], vz9[t]}}/.nds[[1]];
swp[t_]=(m1 Evaluate[f2p[t][[1]]]+m2 Evaluate[f2p[t][[2]]]+m3 Evaluate[f2p[t][[3]]]+m4 Evaluate[f2p[t][[4]]]+m5 Evaluate[f2p[t][[5]]]+m6 Evaluate[f2p[t][[6]]]+m7 Evaluate[f2p[t][[7]]]+m8 Evaluate[f2p[t][[8]]]+m9 Evaluate[f2p[t][[9]]])/(m1+m2+m3+m4+m5+m6+m7+m8+m9);
(* mittlere Abstände *)
table=Table[Subscript[d,{n1,n2}]=Quiet[NIntegrate[Norm[f2p[t][[n1]]-f2p[t][[n2]]]/km,{t,0,Tmax}]]/Tmax "km",{n1,1,9},{n2,1,9}];
MatrixForm[table]
Grid[{
{n2,n1,Subscript[d, {2,1}], 1/2},
{n2,n3,Subscript[d, {2,3}], 1/2},
{n2,n4,Subscript[d, {2,4}], 1/2},
{n2,n5,Subscript[d, {2,5}], 1/2},
{n2,n6,Subscript[d, {2,6}], 1/2},
{n2,n7,Subscript[d, {2,7}], 1/2},
{n2,n8,Subscript[d, {2,8}], 1/2},
{n2,n9,Subscript[d, {2,9}], 1/2},
{},
{n3,n1,Subscript[d, {3,1}], 1/2},
{n3,n2,Subscript[d, {3,2}], 1/2},
{n3,n4,Subscript[d, {3,4}], 1/2},
{n3,n5,Subscript[d, {3,5}], 1/2},
{n3,n6,Subscript[d, {3,6}], 1/2},
{n3,n7,Subscript[d, {3,7}], 1/2},
{n3,n8,Subscript[d, {3,8}], 1/2},
{n3,n9,Subscript[d, {3,9}], 1/2},
{},
{n4,n1,Subscript[d, {4,1}], 1/2},
{n4,n2,Subscript[d, {4,2}], 1/2},
{n4,n3,Subscript[d, {4,3}], 1/2},
{n4,n5,Subscript[d, {4,5}], 1/2},
{n4,n6,Subscript[d, {4,6}], 1/2},
{n4,n7,Subscript[d, {4,7}], 1/2},
{n4,n8,Subscript[d, {4,8}], 1/2},
{n4,n9,Subscript[d, {4,9}], 1/2},
{},
{n5,n1,Subscript[d, {5,1}], 1/2},
{n5,n2,Subscript[d, {5,2}], 1/2},
{n5,n3,Subscript[d, {5,3}], 1/2},
{n5,n4,Subscript[d, {5,4}], 1/2},
{n5,n6,Subscript[d, {5,6}], 1/2},
{n5,n7,Subscript[d, {5,7}], 1/2},
{n5,n8,Subscript[d, {5,8}], 1/2},
{n5,n9,Subscript[d, {5,9}], 1/2},
{},
{n6,n1,Subscript[d, {6,1}], 1/2},
{n6,n2,Subscript[d, {6,2}], 1/2},
{n6,n3,Subscript[d, {6,3}], 1/2},
{n6,n4,Subscript[d, {6,4}], 1/2},
{n6,n5,Subscript[d, {6,5}], 1/2},
{n6,n7,Subscript[d, {6,7}], 1/2},
{n6,n8,Subscript[d, {6,8}], 1/2},
{n6,n9,Subscript[d, {6,9}], 1/2},
{},
{n7,n1,Subscript[d, {7,1}], 1/2},
{n7,n2,Subscript[d, {7,2}], 1/2},
{n7,n3,Subscript[d, {7,3}], 1/2},
{n7,n4,Subscript[d, {7,4}], 1/2},
{n7,n5,Subscript[d, {7,5}], 1/2},
{n7,n6,Subscript[d, {7,6}], 1/2},
{n7,n8,Subscript[d, {7,8}], 1/2},
{n7,n9,Subscript[d, {7,9}], 1/2},
{},
{n8,n1,Subscript[d, {8,1}], 1/2},
{n8,n2,Subscript[d, {8,2}], 1/2},
{n8,n3,Subscript[d, {8,3}], 1/2},
{n8,n4,Subscript[d, {8,4}], 1/2},
{n8,n5,Subscript[d, {8,5}], 1/2},
{n8,n6,Subscript[d, {8,6}], 1/2},
{n8,n7,Subscript[d, {8,7}], 1/2},
{n8,n9,Subscript[d, {8,9}], 1/2},
{},
{n9,n1,Subscript[d, {9,1}], 1/2},
{n9,n2,Subscript[d, {9,2}], 1/2},
{n9,n3,Subscript[d, {9,3}], 1/2},
{n9,n4,Subscript[d, {9,4}], 1/2},
{n9,n5,Subscript[d, {9,5}], 1/2},
{n9,n6,Subscript[d, {9,6}], 1/2},
{n9,n7,Subscript[d, {9,7}], 1/2},
{n9,n8,Subscript[d, {9,8}], 1/2},
{},
{" ", " "}},
Alignment->Left]
Integrationszeit 3 julianische Jahre in die Zukunft:
Da es zu viel RAM frisst wenn ich mehr als 3 Jahre im Stück integriere (einige der Monde machen ja mehrere Umrundungen pro Tag) noch mal 3 Jahre in die Vergangenheit:
Metis und Adrastea haben halt fast den gleichen Bahnradius, und für Kallisto müsste man da er die längste Periode hat wohl noch ein paar zusätzliche Runden integrieren bis sich in seinem System ein eindeutiger Sieger herauskristallisiert, aber man sieht schon dass es ein ziemlich knappes Rennen ist. Ich werde vielleicht später meinen Code noch insofern erweitern dass er automatisch mehrere Runden hintereinander ausführt und dazwischen den Speicher freigibt ohne dass ich daneben sitzen muss, nach mehreren Durchgängen wird wahrscheinlich eh nur Metis übrigbleiben.
Updatend,
