Hallo Alexander,
Ich habe gerade mit meinem Physik Studium begonnen, daher nochnicht viel Erfahrung mit Mathematica.
"daher" ... man kann seit 1988 mit Mma rechnen, dazu muss man nicht vorgängig ein Physikstudium aufnehmen ;-)
Die Fehler sind jetzt weg
Clear[m1, m2, g, l0, d, \[Alpha], \[Phi]]
m1 = 2; m2 = 2; g = 9.81; l0 = 4; d = 10^4; \[Alpha] = 30.21;
Clear[\[Phi], Fsx1, Fsx2, Fsy1, Fsy2]
\[Phi] = d*
E^(-2*\[Alpha]*(Sqrt[(x2[t] - x1[t])^2 + (y2[t] - y1[t])^2] -
l0)) - 2*d*
E^(-\[Alpha]*(Sqrt[(x2[t] - x1[t])^2 + (y2[t] - y1[t])^2] - l0));
Fsx1 = D[\[Phi], x1[t]];
Fsx2 = D[\[Phi], x2[t]];
Fsy1 = D[\[Phi], y1[t]];
Fsy2 = D[\[Phi], y2[t]];
s = NDSolve[{x1''[t] == Fsx1/m1, x1[0] == 0, x1'[0] == 30,
x2''[t] == Fsx2/m2, x2[0] == 50/Sqrt[2], x2'[0] == 0,
y1''[t] == -g + Fsy1/m1, y1[0] == 0, y1'[0] == 70,
y2''[t] == -g + Fsy2/m2, y2[0] == 50/Sqrt[2], y2'[0] == 0},
{x1[t], x2[t], y1[t], y2[t]}, {t, 0, 10}];
die Werte sind gestrichen, sie werden direkt gesetzt, all die
überflüssigen Zweien in den Funktionennamen sind auch weg.
Das Bildchen sieht sehr komisch aus, die Masse 2 stürzt direkt herab, Masse 1 fliegt auf einer Parabel: das Morsepotential ist allenfalls nicht korrekt, beiden Massen erscheinen völlig entkoppelt in der Lösung.
Gruss Udo.
knoettnerHantelsystem.nb
Description: application/vnd.wolfram.mathematica
