Vers pôle nord céleste

Vue dans le plan horizontal de l'observateur

Latitude  \; du \; lieu

Angle \;  du \;  cône

Hauteur \; du \; gnomon

latitude = 0.9250245035569946 rad (slider from 0 to 90*°, step 0.017453292519943295)
hauteur = 2 (slider from 0 to 10, step 0.05)
O: (0, 0, hauteur)
α = 0.5672320068981571 rad (slider from 0*° to 90*°, step 0.5*°)
A: O + 7 * (0, cos(latitude), sin(latitude))
cone = InfiniteCone(O, A, α)
K = (0, 0, 10)
β = Angle(A, O, K)
Axe = Line(O, A)
o = (0, 0, 0)
oxy = PerpendicularPlane(o, Line[O, K])
inter = IntersectionPaths(oxy, cone)
d = Segment(o, O)
u = Vector(O, A)
B = Intersect(Axe, oxy)
a = Segment(B, o)
lat = Angle(O, B, o)