fez
Author: PhM

Rotation P

Ecliptique

α_{Psol} = 286.11
δ_{Psol} = 63.85
β: 90 - δ_{Psol}
L_{Ter} = 180 (slider from 0 to 360, step 1)
ε = 0.4092797095926703 rad
sinδ_T: sin(L_{Ter}°) * sin(ε)
δ_T: asin(sinδ_T) / π * 180
cosα_T: cos(L_{Ter}°) / cos(δ_T°)
sinα_T: cos(ε) * sin(L_{Ter}°)
α_T: x((Angle[(1, 0), (0, 0), (cosα_T, sinα_T)], 0)) * 180 / π / 15
aa: α_T * 15 - α_{Psol} + If[α_T * 15 - α_{Psol} < 0, 360, 0]
cosP_ST: cos(β°) * sin(δ_T°) + cos(δ_T°) * sin(β°) * cos(aa°)
arcP_ST: acos(cosP_ST) * 180 / π
Pcos: (cos(β°) - sin(δ_T°) * cos(arcP_ST°)) / (cos(δ_T°) * sin(arcP_ST°))
Psin: sin(β°) * sin(aa°) / sin(arcP_ST°)
PP: x((Angle[(1, 0), (0, 0), (Pcos, Psin)], 0)) * 180 / π
A = (0.33867264305355704, 0.32057482978877355)
h_{ima} = 2.05
R_{sol} = 1
z_P: y(A) / R_{sol}
y_P: x(A) / R_{sol}
P = If(PP > 180, PP - 360, PP)
c_{sol} = Circle((0, 0), R_{sol})
B_0: 90 - arcP_ST
aP: α_T * 15 - α_{Psol}
x_P: sqrt(1 - y_P^(2) - z_P^(2))
l_P: atan(z_P / sqrt(x_P^(2) + y_P^(2))) * 180 / π
L_P: atan(y_P / x_P) * 180 / π
par_P = Segment((sqrt(1 - z_P^(2)) * R_{sol}, z_P * R_{sol}), (-sqrt(1 - z_P^(2)) * R_{sol}, z_P * R_{sol}))
mer_P = Ellipse((0, sqrt(1 - y_P^(2)) * R_{sol}), (0, -sqrt(1 - y_P^(2)) * R_{sol}), R_{sol})
b: z_P / sqrt(1 - y_P^(2))
ec: y_P / sqrt(1 - z_P^(2))
cel = Rotate(Function[sqrt(R_{sol}^(2) - x^(2)) * ec, -R_{sol}, R_{sol}], -1.5707963267948966, (0, 0))
cima: h_{ima} * sqrt(2) / 2
f_{rot} = false
θ: If[f_{rot}, -P, 0] * π / 180
mrdctr = Segment((0, -R_{sol}), (0, R_{sol}))
equator = Segment((-R_{sol}, 0), (R_{sol}, 0))
e_{equat}: sqrt(1 - sin(B_0°)^(2))
axeP = Rotate(Segment[(0, -h_{ima}), (0, h_{ima})], P°)
A' = Rotate(A, (-P)°)
y'_P: x(A') / R_{sol}
z'_P: y(A') / R_{sol}
ec': y'_P / sqrt(1 - z'_P^(2))
cel2 = Rotate(Function[sqrt(R_{sol}^(2) - x^(2)) * ec', -R_{sol}, R_{sol}], -1.5707963267948966)
par'_P = Segment((sqrt(1 - z'_P^(2)) * R_{sol}, z'_P * R_{sol}), (-sqrt(1 - z'_P^(2)) * R_{sol}, z'_P * R_{sol}))
equatB = Rotate(Function[-sqrt(R_{sol}^(2) - x^(2)) * sin(B_0°), -R_{sol}, R_{sol}], P°, (0, 0))
x'_P: sqrt(1 - y'_P^(2) - z'_P^(2))
L'_P: atan(y'_P / x'_P) * 180 / π
l'_P: atan(z'_P / sqrt(x'_P^(2) + y'_P^(2))) * 180 / π
texte2: "Coordonnées de A' : " + L'_P + "°   " + l'_P + "°"
texte1: "Coordonnées de A :  " + L_P + "°    " + l_P + "°"
texte3: "P = " + P + "°   B0 = " + B_0 + "°"
x_{PS}: R_{sol} * sin(B_0°)
z_{PS}: R_{sol} * cos(B_0°) * cos((-P)°)
y_{PS}: R_{sol} * cos(B_0°) * sin((-P)°)
P_S = If(x_{PS} > 0, (y_{PS}, z_{PS}), (-y_{PS}, -z_{PS}))
qq: sqrt(y_{PS}^(2) + z_{PS}^(2))
ε': asin(sin(ε) * sin((90 + L_{Ter})°)) * 180 / π
sgecl = Rotate(Segment[(-1, 0), (1, 0)], ε'° + θ)
f_{ecl]} = false