1) Newton
Code: Alles auswählen
(* Gegeben *)
G = 667384*^-16 m^3/kg/sek^2; (* Newtons Konstante *)
M1 = 5972 10^21 kg; (* Erdmasse *)
M2 = 1 kg; (* Testmasse *)
M = M1 + M2; (* Gesamtmasse *)
R1 = 6371000 m; (* Erdradius *)
R2 = 1 m; (* Objektradius *)
R = R1 + R2; (* minimaler Schwerpunktabstand *)
h = 10 m; (* Fallhöhe *)
v0 = 0 m/sek; (* Initialgeschwindigkeit *)
(* Gesucht *)
t = f[\!\(
\*SubsuperscriptBox[\(\[Integral]\), \(R\), \(R + h\)]\(
\*FractionBox[\(1\),
SqrtBox[\(2\ G\ M\ \((
\*FractionBox[\(1\), \(r\)] -
\*FractionBox[\(1\), \(R + h\)])\) +
v0\ Abs[v0]\)]] \[DifferentialD]r\)\)] (* Fallzeit, integriert *)
t = f[1/(-2 G M + (h + R) v0 Abs[v0])^(3/2) Sqrt[
h + R] (-Sqrt[2 G h M + R (h + R) v0 Abs[v0]] Sqrt[-2 G M R +
R (h + R) v0 Abs[v0]] + (h +
R) (Sqrt[(h + R) v0] Sqrt[Abs[v0]]
Sqrt[-2 G M + (h + R) v0 Abs[v0]] +
2 G M (ArcSinh[Sqrt[-(R/(h + R)) + (R v0 Abs[v0])/(2 G M)]] -
ArcSinh[
Sqrt[-1 + ((h + R) v0 Abs[v0])/(
2 G M)]])))] (* Fallzeit, analytisch *)
v = f[Sqrt[\!\(
\*SubsuperscriptBox[\(\[Integral]\), \(R\), \(R + h\)]\(\((
\*FractionBox[\(2\ G\ M\),
SuperscriptBox[\(r\), \(2\)]] +
\*FractionBox[\(v0\ Abs[
v0]\), \(h\)])\) \[DifferentialD]r\)\)]] (* insgesamte Geschwindigkeit, integriert *)
v = f[Sqrt[(2 G M h)/(R (R + h)) +
v0 Abs[v0]]] (* insgesamte Geschwindigkeit, analytisch *)
v1 = v M2/(M1 + M2) (* Geschwindigkeit Erde *)
v2 = v M1/(M1 + M2) (* Geschwindigkeit Objekt *)
Zeit t, Integral:
Zeit t, analytisch:
Zeit t, numerisch:
Gesamtgeschwindigkeit v, Integral:
Gesamtgeschwindigkeit v, analytisch:
Gesamtgeschwindigkeit v, numerisch:
Einzelgeschwindigkeiten v1 & v2:
Strecke statt Ziel:
Hier ist h die initiale Höhe und H die bereits gefallene Strecke.
.nb