long geoc

long geoc

long geoc

long geoc

long helioc

du soleil

de Mars

du soleil

de Mars

de la Terre

Date

Jour Julien

L_T(S)

L_T(M)

L_T(S)

L_T(M)

L_H(T)

339°23'

135°12'

359°41'

131°48'

295°21'

182°08'

316°06'

184°42'

 16°50'

168°12'

333°42'

218°48'

185°47'

284°18'

143°26'

346°56'

265°53'

  3°04'

221°42'

 49°42'

Périodes

Terre

Terre

Mars

G12 = 265.88
I12 = Mod(G12 - 180, 360)
e_T = 0.01671
l_T = 1.7966419320029627 rad
x_T: e_T * cos(l_T + 180°)
y_T: e_T * sin(l_T + 180°)
E: (x_T, y_T)
r = 1.5224799999998198 (slider from 1 to 2, step 1.0E-5)
y_C = 0.07163999999999943 (slider from -0.5 to 0.5, step 1.0E-5)
x_C = -0.12045999999999983 (slider from -0.5 to 0.5, step 1.0E-5)
yc_T = -0.016687099525748927
A4 = 1585
B4 = 2
C4 = 17
D4 = 2300016
G4 = 339.38
H4 = 135.2
K4 = 21
L4 = 687
M4 = 708
A5 = 1585
B5 = 3
C5 = 10
D5 = 2300037
G5 = 359.68
H5 = 131.8
K5 = 666
L5 = 687
M5 = 748
A6 = 1587
B6 = 1
C6 = 5
D6 = 2300703
G6 = 295.35
H6 = 182.13
K6 = 21
L6 = 82
M6 = 769
A7 = 1587
B7 = 1
C7 = 26
D7 = 2300724
G7 = 316.1
H7 = 184.7
K7 = 61
L7 = 748
M7 = 1697
A8 = 1587
B8 = 3
C8 = 28
D8 = 2300785
G8 = 16.83
H8 = 168.2
K8 = 687
L8 = 1636
M8 = 2323
A9 = 1589
B9 = 2
C9 = 12
D9 = 2301472
G9 = 333.7
H9 = 218.8
K9 = 949
L9 = 1636
M9 = 1759
A10 = 1591
B10 = 9
C10 = 19
D10 = 2302421
G10 = 185.78
H10 = 284.3
K10 = 687
L10 = 810
M10 = 1497
A11 = 1593
B11 = 8
C11 = 6
D11 = 2303108
G11 = 143.43
H11 = 346.93
K11 = 123
L11 = 810
A12 = 1593
B12 = 12
C12 = 7
D12 = 2303231
H12 = 3.07
K12 = 687
A13 = 1595
B13 = 10
C13 = 25
D13 = 2303918
G13 = 221.7
H13 = 49.7000000000002
D16 = 365.25
E16 = 686.980000000001
I4 = Mod(G4 - 180, 360)
I5 = Mod(G5 - 180, 360)
I6 = Mod(G6 - 180, 360)
I7 = Mod(G7 - 180, 360)
I8 = Mod(G8 - 180, 360)
I9 = Mod(G9 - 180, 360)
I10 = Mod(G10 - 180, 360)
I11 = Mod(G11 - 180, 360)
I13 = Mod(G13 - 180, 360)
S = (0, 0)
γ = (2, 0)
d_1 = Ray(S, Vector[S, (cos(I4 * π / 180), sin(I4 * π / 180))])
C_T = Circle((x_T, y_T), 1)
T_1 = Intersect(C_T, Ray[S, Vector[S, (cos(I4 * π / 180), sin(I4 * π / 180))]])
T'_1 = Intersect(C_T, Ray[S, Vector[S, (cos(I6 * π / 180), sin(I6 * π / 180))]])
T_2 = Intersect(C_T, Ray[S, Vector[S, (cos(I5 * π / 180), sin(I5 * π / 180))]])
T'_2 = Intersect(C_T, Ray[S, Vector[S, (cos(I7 * π / 180), sin(I7 * π / 180))]])
T_3 = Intersect(C_T, Ray[S, Vector[S, (cos(I8 * π / 180), sin(I8 * π / 180))]])
T'_3 = Intersect(C_T, Ray[S, Vector[S, (cos(I9 * π / 180), sin(I9 * π / 180))]])
T_4 = Intersect(C_T, Ray[S, Vector[S, (cos(I10 * π / 180), sin(I10 * π / 180))]])
T'_4 = Intersect(C_T, Ray[S, Vector[S, (cos(I11 * π / 180), sin(I11 * π / 180))]])
T'_5 = Intersect(C_T, Ray[S, Vector[S, (cos(I13 * π / 180), sin(I13 * π / 180))]])
m_1 = Ray(T_1, Vector[S, (cos(H4 * π / 180), sin(H4 * π / 180))])
m'_1 = Ray(T'_1, Vector[S, (cos(H6 * π / 180), sin(H6 * π / 180))])
M_1 = Intersect(m_1, m'_1)
m'_2 = Ray(T'_2, Vector[S, (cos(H7 * π / 180), sin(H7 * π / 180))])
m_2 = Ray(T_2, Vector[S, (cos(H5 * π / 180), sin(H5 * π / 180))])
M_2 = Intersect(m_2, m'_2)
m_3 = Ray(T_3, Vector[S, (cos(H8 * π / 180), sin(H8 * π / 180))])
m'_3 = Ray(T'_3, Vector[S, (cos(H9 * π / 180), sin(H9 * π / 180))])
M_3 = Intersect(m'_3, m_3)
m'_4 = Ray(T'_4, Vector[S, (cos(H11 * π / 180), sin(H11 * π / 180))])
m_4 = Ray(T_4, Vector[S, (cos(H10 * π / 180), sin(H10 * π / 180))])
M_4 = Intersect(m_4, m'_4)
T_5 = Intersect(C_T, Ray[S, Vector[S, (cos(I12 * π / 180), sin(I12 * π / 180))]])
m_5 = Ray(T_5, Vector[S, (cos(H12 * π / 180), sin(H12 * π / 180))])
m'_5 = Ray(T'_5, Vector[S, (cos(H13 * π / 180), sin(H13 * π / 180))])
M_5 = Intersect(m_5, m'_5)
C: (x_C, y_C)
C_M = Circle(C, r)
I_1 = Intersect(C_M, Ray[C, M_1])
I_2 = Intersect(C_M, Ray[C, M_2])
I_3 = Intersect(C_M, Ray[C, M_3])
I_4 = Intersect(C_M, Ray[C, M_4])
I_5 = Intersect(C_M, Ray[C, M_5])
l_1 = Segment(M_1, I_1)
l_2 = Segment(M_2, I_2)
l_3 = Segment(M_3, I_3)
l_4 = Segment(M_4, I_4)
l_5 = Segment(M_5, I_5)
Exces: l_1^(2) + l_2^(2) + l_3^(2) + l_4^(2) + l_5^(2)
gaxe = Line(C, S)
J_1, J_2 = Intersect(C_M, gaxe)
a = Segment(C, J_1)
e: Segment[C, S] / a
texte1: "E=" + Exces
Lper = Angle(Vector[(0, 0), (1, 0)], Vector[C, S])
texte3: "demi-grand axe a=" + a
texte4: "Excentricité e=" + e
texte5: "Longitude du périhélie Lper=" + Lper
texte2: "Excès=" + Exces + ""

Tool: Terre
Input: S, I4
Output: T_1

C_T = Circle(S, 1)
d_1 = Ray(S, Vector[S, (cos(I4 * π / 180), sin(I4 * π / 180))])
T_1 = Intersect(C_T, d_1)

Tool: ddte
Input: S, G4
Output: d_1

I4 = Mod(G4 - 180, 360)
d_1 = Ray(S, Vector[S, (cos(I4 * π / 180), sin(I4 * π / 180))])

Tool: ddtem
Input: S, G4, H4
Output: m_1

I4 = Mod(G4 - 180, 360)
C_T = Circle(S, 1)
d_1 = Ray(S, Vector[S, (cos(I4 * π / 180), sin(I4 * π / 180))])
T_1 = Intersect(C_T, d_1)
m_1 = Ray(T_1, Vector[S, (cos(H4 * π / 180), sin(H4 * π / 180))])