jour de l'année

Période de l'année

Est

Ouest

Levant

Couchant

angles de vue figure

Solstice \;  d'hiver

solstice d'été

équinoxe de printemps

équinoxe d'automne

   \;  \;  \; Soleil

Trace de l'ombre

Déclinaison

gnomon vertical

gnomon polaire

latitude = 54 (slider from -90 to 90, step 0.1)
latitude_1: latitude°
α = 1.8325954499999997 rad (slider from 0 to 3.141592653589793, step 0.9999998556178845*°)
O = (0, 0, 0)
β = 0 rad (slider from 0*° to 20*°, step 0.017453292519943295)
I': (cos(α) * cos(β), -sin(α), -cos(α) * sin(β))
u = Vector(O, I')
K': (sin(β), 0, cos(β))
w = Vector(O, K')
verspolaire: -u * cos(latitude_1) + sin(latitude_1) * w
dirbaton: w
taillebaton = 3 (slider from 0 to 10, step 1)
heure = 10.033333534 (slider from 0 to 24, step 0.016666667)
azimut: (heure - 12) * 15 * π / 180
echelle = 1 (slider from 0 to 5, step 0.1)
ptverspolaire = Translate(O, verspolaire)
axerotationterre = Line(O, ptverspolaire)
lattrop = 23.5
jour = 173 (slider from 0 to 365, step 1)
declinaison: lattrop * sin((jour - 81) * π / 182.5)
dirculminsol: echelle * (cos(π / 2 - latitude_1 + declinaison * π / 180) * u + sin(π / 2 - latitude_1 + declinaison * π / 180) * w)
ptculminsol = Translate(O, dirculminsol)
f = Circle(axerotationterre, ptculminsol)
Soleil_1 = Point(f)
ecliptique = PerpendicularPlane(Soleil_1, axerotationterre)
O'' = Intersect(ecliptique, axerotationterre)
Soleil: O'' + cos(-azimut) * Vector[O'', ptculminsol] + sin(-azimut) * verspolaire ⊗ Vector[O'', ptculminsol] / Length[verspolaire]
I = (1, 0, 0)
J = (0, 1, 0)
K = (0, 0, 1)
J': (sin(α) * cos(β), cos(α), -sin(α) * sin(β))
a = Line(O, I')
b = Line(O, K')
c = Line(O, J')
v = Vector(O, J')
d = Circle(O, echelle, b)
Sud, Nord = Intersect(a, d)
W, E = Intersect(c, d)
Intersect()
C, L = Intersect(f, d)
O' = Translate(O, Vector[taillebaton / 10 * dirbaton])
g = Segment(O, O')
i = PerpendicularPlane(O, b)
k = Ray(Soleil, O')
A = Intersect(k, i)
h = Segment(O, A)
l = Segment(Sud, Nord)
m = Segment(W, E)
n = PerpendicularPlane(O, u)
o = PerpendicularPlane(O, w)
TracedeA: (x(A), y(A))
affichageombre: Vector[O, Soleil] * Vector[O, A]
Aff = (1.3227497902601857, 1.1487285092444464, 5)
j = Segment(Soleil, A)
B = (-4.54879134092208, 1.9062609242822368)
e = false
texte3: "" + (jour * " \; ème")
p = true
texte4: "Déclinaison \;  " + (LaTeX[declinaison]) + "°"
γ = Angle(Nord, O, O'')