Das Prinzip des Seins

ClearAll["Global`*"]
tmax = 1/Sqrt[2];
f[t_] := 1/Sin[t π + π/4]/2 + t; F[t_] := -(f[t] - t) + f[0];
plot2[f_, ta_, tb_, c1_, c2_, c3_, a_, b_, c_] := Do[Print[
P1 = {t Sin[n π/180 + π/4], t Cos[n π/180 + π/4], 0};
L1 = {P1, {0, 0, 0}};
TP1 = {c3, Table[{PointSize[0.01], Point[P1]}, {n, 0, 45}]};
TL1 = Table[Line[L1, VertexColors -> {c2, c1}], {n, 0, 45}];
P2[t_] := {(tmax - t) Sin[n π/180 + π/4], (tmax - t) Cos[n π/180 + π/4], 0};
L2 = {P2[t], p}; p = P2[0];
TP2 = {c3, Table[{PointSize[0.01], Point[P2[t]]}, {n, 0, 45}]};
TL2 = Table[Line[L2, VertexColors -> {c1, c2}], {n, 0, 45}];
P3 = {(t - F[n/180]) Sin[n π/180 + π/4], (t - F[n/180]) Cos[n π/180 + π/4], 0};
L3 = {P3, {0, 0, 0}};
TP3 = {c3, Table[{PointSize[0.01], Point[P3]}, {n, 0, 45}]};
TL3 = Table[Line[L3, VertexColors -> {c2, c1}], {n, 0, 45}];
Rasterize[Grid[{{
Grid[{{Rotate[
Graphics3D[{TL1, TP1, {Dashed, Line[{{1/2, 1/2, 0}, {1/2, -1/2, 0}}]}},
PlotRange -> {{1/2, 0}, {0, 1/2}, {-1/16, 1/16}},
ImagePadding -> 1, ImageSize -> 240,
ViewPoint -> {a, b, c}], -π/2]}, {N[2 t]}}, Alignment -> Left],
Grid[{{Rotate[
Graphics3D[{TL3, TP3, {Dashed, Line[{{1/2, 1/2, 0}, {1/2, -1/2, 0}}]}},
PlotRange -> {{1/2, 0}, {0, 1/2}, {-1/16, 1/16}},
ImageSize -> 240, ImagePadding -> 1,
ViewPoint -> {a, b, c}], -π/2]}, {"Projektionseffekt"}}, Alignment -> Center],
Grid[{{Rotate[
Graphics3D[{TL2, TP2, {Dashed, Line[{{1/2, 1/2, 0}, {1/2, -1/2, 0}}]}},
PlotRange -> {{1/2, 0}, {0, 1/2}, {-1/16, 1/16}},
ImageSize -> 240, ImagePadding -> 1,
ViewPoint -> {a, b, c}], -π/2]}, {"yukterez.net"}}, Alignment -> Right]}}]]],
{t, ta, tb, 1/Sqrt[2]/10}];
p3 = plot2[1, 0, 1/Sqrt[2], Red, Magenta, Purple, 0, 0, Infinity]

ClearAll["Global`*"]
G = 1; M = 1; c = 1; rs = 2 G M/c^2;
wp = MachinePrecision;
j[v_] := Sqrt[1 - v^2/c^2];
J = j[v0];
k[r_] := Sqrt[1 - rs/r];
к = k[r0];
r0 = 20;
θ0 = -Pi/4;
тmax = 200;
Ф = β Pi/18000;
vr0 = v0 Sin[Ф] к/J;
vθ0 = v0/r0 Cos[Ф]/J;
v0 = 999999/1000000;
Table[Subscript[sol, β] =
NDSolve[{
r''[t] == -((G M)/r[t]^2) + r[t] θ'[t]^2 - (3 G M)/c^2 θ'[t]^2,
r'[0] == vr0,
r[0] == r0,
θ''[t] == -((2 r'[t] θ'[t])/r[t]),
θ'[0] == vθ0, θ[0] == θ0,
τ'[t] == Sqrt[c^2 r[t] + r[t] r'[t]^2 - c^2 rs + r[t]^3 θ'[t]^2 - r[t]^2 rs θ'[t]^2]/(c Sqrt[r[t] - rs] Sqrt[1 - rs/r[t]]),
τ[0] == 0,
cl'[t] == ((r'[t]/k[r[t]])^2 + (θ'[t] r[t])^2)/c^2,
cl[0] == 0}, {r, θ, τ, cl}, {t, 0, тmax},
MaxSteps -> Infinity, Method -> Automatic, WorkingPrecision -> wp, InterpolationOrder -> All],
{β, 25565, 25575, 0.1}];
t[Χ_, β_] := Quiet[ξ /. FindRoot[
Evaluate[τ[ξ] /. Subscript[sol, β]][[1]] - Χ, {ξ, 0},
WorkingPrecision -> wp, Method -> Automatic]];
Τ[β_] := Quiet[t[ι, β]];
x[t_, β_] := (Sin[Evaluate[θ[t] /. Subscript[sol, β]]] Evaluate[r[t] /. Subscript[sol, β]])[[1]]
y[t_, β_] := (Cos[Evaluate[θ[t] /. Subscript[sol, β]]] Evaluate[r[t] /. Subscript[sol, β]])[[1]]
R[t_, β_] := Evaluate[r[t] /. Subscript[sol, β]][[1]];
γ[t_, β_] := Evaluate[τ'[t] /. Subscript[sol, β]][[1]];
и[t_, β_] := Evaluate[τ[t] /. Subscript[sol, β]][[1]];
crθ[t_, β_] := Evaluate[cl'[t] /. Subscript[sol, β]][[1]];
vrθ[t_, β_] := crθ[t]/Sqrt[1 + crθ[t]^2];
clr[t_, β_] := Evaluate[r'[t] /. Subscript[sol, β]][[1]];
clθ[t_, β_] := R[t] Evaluate[θ'[t] /. Subscript[sol, β]][[1]];
s[text_] := Style[text, FontSize -> font]; font = 11;
PR = Sqrt[2] 10;
Do[Print[
Rasterize[Grid[{{Show[
Graphics[{{LightGray, Disk[{0, 0}, rs]}, {Lighter[Gray], Dashed, Circle[{0, 0}, r0]}},
Frame -> True, ImageSize -> 400, PlotRange -> PR, ImagePadding -> 1],
Table[
{Graphics[{RGBColor[1 - (25573 - β)/10, (25575 - β)/10, 0],
Point[{x[Τ[β], β], y[Τ[β], β]}]},
PlotRange -> {{-10, 10}, {-10, 10}}],
ParametricPlot[{x[η, β], y[η, β]}, {η, 0, Τ[β]},
PlotStyle -> {Thickness[0.001],
RGBColor[1 - (25575 - β)/10, (25575 - β)/10, 0]}]},
{β, 25565, 25573, 0.1}]]},
{Grid[{
{s["t"], "=", s[N[ι]], s[" GM/c³"]}}, Alignment -> Left,
Spacings -> {0, 1/2}]}}, Alignment -> Left]]],
{ι,10, 120, 10}]