Die Frage ist nur, was wollt ihr Großmeister der Kraftausdrücke denn schon groß mit einer solchen Gleichung anfangen?
Code:
settings = {AccuracyGoal -> 24, PrecisionGoal -> 24,
WorkingPrecision -> 84, MaxSteps -> \[Infinity]};
kg = 1; m = 1; sek = 1;
G = 667384*^-16 m^3 kg^-1 sek^-2; (*Gravitation*)
c = 299792458 m/sek; (*Licht*)
yr = 36525*24*36 sek; (*Jahr*)
Au = 149597870700 m; (*Astronomische Einheit*)
r = 696000000 m; (*Radius Sonne*)
rs = 2 G M/c^2; (*Schwarzschildradius*)
ps = 3 G M/c^2; (*Photonensphäre*)
r0 = r; (*Distanz der kleinsten Annäherung*)
x1 = 12 ps; (*radiale Distanz der Quelle zur Masse*)
x2 = x1; (*radiale Distanz Beobachter zur Masse*)
M = 19891*^26 kg; (*Masse der Linse*)
b = y[x]/Sqrt[y[x]^2 + (x)^2]; (*transversale Komponente*)
v = c/Sqrt[1^2 + Tan[\[Delta][x]]^2]; (*radiale Geschwindigkeitskomponente*)
ds = Sqrt[y[x]^2 + x^2]; (*Abstandsbetrag*)
sq = Sqrt[1 - rs/ds]; (*Schwarzschildfaktor*)
n = 1 + (2*G*M)/(ds c^2)/sq;(*Brechungsindex*)
w = Abs[(2*G*M*b)/ds^2/sq/v^2];(*Winkeldifferential*)
X = x1; (* Plot Range *)
DGL1 = {\[Delta]'[x] == w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == 1 ps};
DGL2 = {\[Delta]'[x] == w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == 1.5 ps};
DGL3 = {\[Delta]'[x] == w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == 2 ps};
DGL4 = {\[Delta]'[x] == w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == 2.5 ps};
DGL5 = {\[Delta]'[x] == w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == 3 ps};
DGL6 = {\[Delta]'[x] == w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == 3.5 ps};
DGL7 = {\[Delta]'[x] == w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == 4 ps};
DGL1 = {\[Delta]'[x] == -w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == - 1 ps};
DGL12 = {\[Delta]'[x] == -w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == -1.5 ps};
DGL13 = {\[Delta]'[x] == -w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == -2 ps};
DGL14 = {\[Delta]'[x] == -w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == -2.5 ps};
DGL15 = {\[Delta]'[x] == -w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == -3 ps};
DGL16 = {\[Delta]'[x] == -w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == -3.5 ps};
DGL17 = {\[Delta]'[x] == -w, \[Delta][0] == 0,
y'[x] == -Tan[\[Delta][x]], y[0] == -4 ps};
(*Differentialgleichung*)
s1 = NDSolve[DGL1, {\[Delta], y}, {x, -x1, x1}, settings];
s2 = NDSolve[DGL2, {\[Delta], y}, {x, -x1, x1}, settings];
s3 = NDSolve[DGL3, {\[Delta], y}, {x, -x1, x1}, settings];
s4 = NDSolve[DGL4, {\[Delta], y}, {x, -x1, x1}, settings];
s5 = NDSolve[DGL5, {\[Delta], y}, {x, -x1, x1}, settings];
s6 = NDSolve[DGL6, {\[Delta], y}, {x, -x1, x1}, settings];
s7 = NDSolve[DGL7, {\[Delta], y}, {x, -x1, x1}, settings];
s11 = NDSolve[DGL11, {\[Delta], y}, {x, -x1, x1}, settings];
s12 = NDSolve[DGL12, {\[Delta], y}, {x, -x1, x1}, settings];
s13 = NDSolve[DGL13, {\[Delta], y}, {x, -x1, x1}, settings];
s14 = NDSolve[DGL14, {\[Delta], y}, {x, -x1, x1}, settings];
s15 = NDSolve[DGL15, {\[Delta], y}, {x, -x1, x1}, settings];
s16 = NDSolve[DGL16, {\[Delta], y}, {x, -x1, x1}, settings];
s17 = NDSolve[DGL17, {\[Delta], y}, {x, -x1, x1}, settings];
Quiet[Plot[{Flatten[({y[x]} /. s1)][[1]],
Flatten[({y[x]} /. s2)][[1]], Flatten[({y[x]} /. s3)][[1]],
Flatten[({y[x]} /. s4)][[1]], Flatten[({y[x]} /. s5)][[1]],
Flatten[({y[x]} /. s6)][[1]], Flatten[({y[x]} /. s7)][[1]],
Flatten[({y[x]} /. s11)][[1]], Flatten[({y[x]} /. s12)][[1]],
Flatten[({y[x]} /. s13)][[1]], Flatten[({y[x]} /. s14)][[1]],
Flatten[({y[x]} /. s15)][[1]], Flatten[({y[x]} /. s16)][[1]],
Flatten[({y[x]} /. s17)][[1]]}, {x, -x1, x1},
GridLines -> {Table[k rs, {k, -15, 15}], Table[k rs, {k, -15, 15}]},
Frame -> True, AspectRatio -> 1/2, ImageSize -> 600,
PlotRange -> {{-x1/2, x1/2}, {-x1/4, x1/4}},
PlotStyle -> {{Red, Thick}, {Orange, Thick}, {Green, Thick}, {Cyan,
Thick}, {Blue, Thick}, {Purple, Thick}, {Pink, Thick}}]]
Quit[] (* yukterez.net *)

Gridabstand: 1 rs; Photonensphäre bei 1.5 rs如果你是一个白痴每个发票听起来中国,
