Déplacer M pour modifier la latitude

α = 0.12049944426401576 rad (slider from 0 to 6.283185307179586, step 0.98630137°)
A = (0, 0, 0)
d = Sphere(A, 1)
inclinaison = 0.40898882138400455 rad
M = PointIn(d)
u: (0, sin(inclinaison), cos(inclinaison))
H = Intersect(OrthogonalPlane[M, u], Line[A, u])
a = OrthogonalPlane(M, u)
b = Line(A, u)
c = Circle(b, M)
HiverNord = (0, -2, 0)
e = Circle(zAxis, HiverNord)
Soleil: ((2 * sin(α)), (-2 * cos(α)), 0)
f = IntersectPath(OrthogonalPlane[A, Vector[A, Soleil]], d)
Lever = Intersect(f, c, 1)
Coucher = Intersect(f, c, 2)
β = Angle(Lever, H, Coucher)
βi = Angle(Lever, H, Coucher)
gd: α > 90° ∧ α < 270°
numrect: (α / ° / 360 * 365) - 10
nummois = Sequence(Sum[First[{31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31}, k]], k, 1, 12)
mois = If(numrect ≤ 0, 12, CountIf[x < numrect, nummois] + 1)
jour = If(numrect ≤ 0, numrect + 31, numrect - Sum[First[{31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31}, mois - 1]])
latitude: "à la latitude " + (If[(Vector[A, H] * u) ≥ 0, asind((Vector[A, H] * u)) + " Nord", asind((Vector[H, A] * u)) + " Sud"])
dj = If((Vector[A, H] * u) ≥ 0, If[gd, (β / ° * 24 / 360), (βi / ° * 24 / 360)], If[¬gd, (β / ° * 24 / 360), (βi / ° * 24 / 360)])
B = (-29, -62)
dureejour: "la durée du jour est de : " + (floor(dj)) + " h " + (floor(((dj - floor(dj)) * 60))) + " min"
text3: "On est le : " + jour + " " + (Element[{"janvier", "février", "mars", "avrl", "mai", "juin", "juillet", "août", "septembre", "octobre", "novembre", "décembre"}, mois]) + " d'une année à 365 jours," + latitude + " " + dureejour + ""
C = (-29, -22)
D = (0, -62)
BD = Distance(B, D)
BC = Distance(B, C)
F: B + ((α / ° / 360 * BD), (dj / 24 * BC))
michel: (βi / ° * 24 / 360)
E = Midpoint(C, D)
B' = Mirror(B, E)
poly1, c_1, b_1, d_1, b' = Polygon(C, B, D, B')