
hat geschrieben:
Vorschlag: Wien
Soll mir auch recht sein. Das Ereignis wird dann ungefähr in 161 Tagen, also am 11.10.2015 und in 304 Tagen, also am 02.03.2016 eintreten (heute war Tag 1 am Plot). Die Achsen am Bild sind in Grad gehalten:

Code:
(* Mathematica Syntax. Yukterez 2015 *)
grad = π/180; (* Winkelmaß *)
φ = 482081743/1*^7 grad; (* Breitengrad Wien *)
λ = 163738189/1*^7 grad; (* Längengrad Wien *)
n = 5601 + d; JD = n + 2451545; (* Tage seit 1.1.2000 *)
l = 28046/100 grad + 9856474/1*^7 grad n; (* Mittlere ekliptikale Länge *)
L = Mod[l, 2 Pi]; (* Skalierung *)
g = 357528/1000 grad + 9856003/1*^7 grad n; (* mittlere Anomalie *)
г = Mod[g, 2 Pi]; (* Skalierung *)
Λ = L + 1915/1000 grad Sin[г] + 1997/1*^5 grad Sin[2 г]; (* Ekliptikale Länge *)
ε = 23439/1000 grad - 4*^-7 grad n; (* Ekliptik *)
α = ArcTan[Cos[ε] Sin[Λ]/
Cos[Λ]] + π; (* Rektaszension *)
δ = ArcSin[Sin[ε] Sin[Λ]]; (* Deklination *)
t0 = n/36525; (* Start *)
θhG = 6697376*^-6 + 240005134*^-5 t0 + 1002738*^-6 T;
θG = (15 Mod[θhG, 24]) π/180; θ = θG + λ; τ = θ - α; (* Stundenwinkel *)
a = ArcTan[Sin[τ]/(Cos[τ] Sin[φ] - Tan[δ] Cos[φ])]; (* Azimut *)
h = ArcSin[Cos[δ] Cos[τ] Cos[φ] + Sin[δ] Sin[φ]]; (* Höhe *)
Do[Print[Show[
Table[ParametricPlot[{a 180/Pi + z, h 180/Pi}, {T, 0, 24},
Frame -> True, PlotRange -> {{-180, 180}, {-90, 90}}], {z, -180,
180, 180}],
Graphics[Point[{176.2, 34.5}]]]], {d, 0, 365}]
Gerne helfend,
