variables orbites
Author: aaa

lg périh

année périh

mois péri

jour péri

heure péri

min péri

excentricité

a en ua

Terre

Mars

jj périhélie

Année

Mois

Jour

Heures

Minutes

n°jour

jour

Période

Mars

Terre

Jour Julien 

Mercure

Vénus

Jupiter

Saturne

Uranus

Neptune

a en km

Mercure

Vénus

u0

263°52'

89°6'

139°32'

352°43'

324°8'

u (anomalie exentrique)

v (anomalie vraie)

val IMCCE au 21/12/2009

Date origine

Date en cours

324°5'

184°0'

θ Longitude héliocentrique au jour en cours

15°21'

Mercure, ..., Mars

C20 = 4
D7 = 365.256
D17 = 2009
D18 = 1
D19 = 4
D20 = 15
D21 = 20
D22: 367.0 * D17 - floor(7.0 * (D17 + floor((D18 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(D17 + floor((D18 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * D18 / 9.0) + D19 + 1721029.0 - 0.5 + D20 / 24.0 + D21 / 24.0 / 60.0
B33 = 12
B34 = 21
B35 = 0
B36 = 0
B32 = 2009
B37: 367.0 * B32 - floor(7.0 * (B32 + floor((B33 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(B32 + floor((B33 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * B33 / 9.0) + B34 + 1721029.0 - 0.5 + B35 / 24.0 + B36 / 24.0 / 60.0
tps = 769 (slider from 0 to 6000, step 0.1)
C37: B37 + tps
D26: 2.0 * 3.141592653589793 * ($C$37 - D22) / D7
D6 = 0.017
D4 = 1.4959807E8
D8 = 1.796589572125403 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
E7 = 686.98
E21 = 46
E20 = 9
E18 = 4
E19 = 21
E17 = 2009
E22: 367.0 * E17 - floor(7.0 * (E17 + floor((E18 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(E17 + floor((E18 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * E18 / 9.0) + E19 + 1721029.0 - 0.5 + E20 / 24.0 + E21 / 24.0 / 60.0
E26: 2.0 * 3.141592653589793 * ($C$37 - E22) / E7
E6 = 0.09341233
E8 = 5.865353484252144 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
B6 = 0.20563069
B7 = 87.969
B17 = 2009
B18 = 10
B19 = 4
B21 = 0
B20 = 13
B22: 367.0 * B17 - floor(7.0 * (B17 + floor((B18 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(B17 + floor((B18 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * B18 / 9.0) + B19 + 1721029.0 - 0.5 + B20 / 24.0 + B21 / 24.0 / 60.0
B26: 2.0 * 3.141592653589793 * ($C$37 - B22) / B7
fT(x) = x + (D26 + D6 * sin(x) - x) / (1.0 - D6 * cos(x))
alp: floor((floor(C37 + 0.5) - 1867216.25) / 36524.25)
a: floor(C37 + 0.5) + 1.0 + alp - floor(alp / 4.0)
b: a + 1524.0
c: floor((b - 122.1) / 365.25)
d: floor(365.25 * c)
e: floor((b - d) / 30.6001)
C32: c - If[B33 > 2.5, 4716.0, 4715.0]
C6 = 0.0067
C21 = 0
C18 = 10
C19 = 4
C17 = 2009
C22: 367.0 * C17 - floor(7.0 * (C17 + floor((C18 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(C17 + floor((C18 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * C18 / 9.0) + C19 + 1721029.0 - 0.5 + C20 / 24.0 + C21 / 24.0 / 60.0
C7 = 224.701
C26: 2.0 * 3.141592653589793 * ($C$37 - C22) / C7
B8 = 1.351862225424728 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
C8 = 2.295509392100502 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
hr: (C37 + 0.5 - floor(C37 + 0.5)) * 24.0
u_T = Iteration(fT(x), D26, 4.0)
v_T: 2.0 * atan(tan(u_T / 2.0) * sqrt((1.0 + D6) / (1.0 - D6))) * 180.0° / 3.141592653589793
D5 = 1
ρ_T: D5 * (1.0 - D6^(2.0)) / (1.0 + D6 * cos(v_T))
θ_T: v_T + D8
T: (ρ_T; θ_T)
H = (0, 0)
fM(x) = x + (E26 + E6 * sin(x) - x) / (1.0 - E6 * cos(x))
F2_T: (2.0 * D5 * D6; D8 + 180.0°)
u_M = Iteration(fM(x), E26, 4.0)
v_M: 2.0 * atan(tan(u_M / 2.0) * sqrt((1.0 + E6) / (1.0 - E6))) * 180.0° / 3.141592653589793
θ_M: v_M + E8
E5 = 1.52368
ρ_M: E5 * (1.0 - E6^(2.0)) / (1.0 + E6 * cos(v_M))
M: (ρ_M; θ_M)
F2_M: (2.0 * E5 * E6; E8 + 180.0°)
C35: floor((C37 + 0.5 - floor(C37 + 0.5)) * 24.0)
C36: floor((hr - C35) * 60.0)
C34: floor(b - d - floor(30.6001 * e) + Mod[C37 + 0.5, 1.0])
C33: e - If[e < 13.5, 1.0, 13.0]
texte1: "Date : " + C34 + "/" + C33 + "/" + C32 + "    " + C35 + "h" + C36 + "min"
D34 = Mod(floor(C37 - 0.5) + 0.5 + 1.5, 7.0)
D36 = Element({"Dimanche", "Lundi", "Mardi", "Mercredi", "Jeudi", "Vendredi", "Samedi"}, D34 + 1.0)
texte3: "" + D36
el_M = Ellipse(H, F2_M, E5)
el_T = Ellipse(H, F2_T, D5)
B5 = 0.38709893
C5 = 0.723332
B4: B5 * $D$4
C4: C5 * $D$4
E4: E5 * $D$4
F2_{Me}: (2.0 * B5 * B6; B8 + 180.0°)
fMe(x) = x + (B26 + B6 * sin(x) - x) / (1.0 - B6 * cos(x))
u_{Me} = Iteration(fMe(x), B26, 4.0)
v_{Me}: 2.0 * atan(tan(u_{Me} / 2.0) * sqrt((1.0 + B6) / (1.0 - B6))) * 180.0° / 3.141592653589793
θ_{Me}: v_{Me} + B8
ρ_{Me}: B5 * (1.0 - B6^(2.0)) / (1.0 + B6 * cos(v_{Me}))
Me: (ρ_{Me}; θ_{Me})
el_{Me} = Ellipse(H, F2_{Me}, B5)
F2_V: (2.0 * C5 * C6; C8 + 180.0°)
fV(x) = x + (C26 + C6 * sin(x) - x) / (1.0 - C6 * cos(x))
u_V = Iteration(fV(x), C26, 4.0)
v_{V}: 2.0 * atan(tan(u_V / 2.0) * sqrt((1.0 + C6) / (1.0 - C6))) * 180.0° / 3.141592653589793
el_{V} = Ellipse(H, F2_V, C5)
θ_{V}: v_{V} + C8
ρ_V: C5 * (1.0 - C6^(2.0)) / (1.0 + C6 * cos(v_{V}))
V: (ρ_V; θ_{V})
F5 = 5.20336301
F6 = 0.04839266
F7 = 4335.3545
G7 = 10757.7365
G5 = 9.53707032
G6 = 0.0541506
H5 = 19.22941195
H6 = 0.044405586
H7 = 30799.095
F4: F5 * $D$4
G4: G5 * $D$4
H4: H5 * $D$4
I5 = 30.06896
I4: I5 * $D$4
I6 = 0.00858587
I7 = 60224.9036
G8 = 1.6182151772423305 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
H8 = 2.9763274267184503 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
I8 = 0.7848920179143699 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
F17 = 2011
F18 = 3
F19 = 17
F20 = 17
F21 = 0
F22: 367.0 * F17 - floor(7.0 * (F17 + floor((F18 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(F17 + floor((F18 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * F18 / 9.0) + F19 + 1721029.0 - 0.5 + F20 / 24.0 + F21 / 24.0 / 60.0
F26: 2.0 * 3.141592653589793 * ($C$37 - F22) / F7
G17 = 2003
G18 = 7
G19 = 26
G20 = 0
G21 = 0
G22: 367.0 * G17 - floor(7.0 * (G17 + floor((G18 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(G17 + floor((G18 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * G18 / 9.0) + G19 + 1721029.0 - 0.5 + G20 / 24.0 + G21 / 24.0 / 60.0
G26: 2.0 * 3.141592653589793 * ($C$37 - G22) / G7
F8 = 2.0125391604746614 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
F2_J: (2.0 * F5 * F6; F8 + 180.0°)
fJ(x) = x + (F26 + F6 * sin(x) - x) / (1.0 - F6 * cos(x))
u_J = Iteration(fJ(x), F26, 10.0)
v_{J}: 2.0 * atan(tan(u_J / 2.0) * sqrt((1.0 + F6) / (1.0 - F6))) * 180.0° / 3.141592653589793
el_{J} = Ellipse(H, F2_J, F5)
ρ_J: F5 * (1.0 - F6^(2.0)) / (1.0 + F6 * cos(v_{J}))
F2_S: (2.0 * G5 * G6; G8 + 180.0°)
fS(x) = x + (G26 + G6 * sin(x) - x) / (1.0 - G6 * cos(x))
u_S = Iteration(fS(x), G26, 4.0)
v_{S}: 2.0 * atan(tan(u_S / 2.0) * sqrt((1.0 + G6) / (1.0 - G6))) * 180.0° / 3.141592653589793
el_{S} = Ellipse(H, F2_S, G5)
θ_{S}: v_{S} + G8
ρ_S: G5 * (1.0 - G6^(2.0)) / (1.0 + G6 * cos(v_{S}))
S: (ρ_S; θ_{S})
H17 = 1965
H18 = 7
H19 = 26
H20 = 0
H21 = 0
H22: 367.0 * H17 - floor(7.0 * (H17 + floor((H18 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(H17 + floor((H18 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * H18 / 9.0) + H19 + 1721029.0 - 0.5 + H20 / 24.0 + H21 / 24.0 / 60.0
I17 = 2040
I18 = 7
I19 = 26
I20 = 0
I21 = 0
I22: 367.0 * I17 - floor(7.0 * (I17 + floor((I18 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(I17 + floor((I18 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * I18 / 9.0) + I19 + 1721029.0 - 0.5 + I20 / 24.0 + I21 / 24.0 / 60.0
H26: 2.0 * 3.141592653589793 * ($C$37 - H22) / H7
I26: 2.0 * 3.141592653589793 * ($C$37 - I22) / I7
F2_U: (2.0 * H5 * H6; H8 + 180.0°)
fU(x) = x + (H26 + H6 * sin(x) - x) / (1.0 - H6 * cos(x))
u_U = Iteration(fU(x), H26, 4.0)
v_{U}: 2.0 * atan(tan(u_U / 2.0) * sqrt((1.0 + H6) / (1.0 - H6))) * 180.0° / 3.141592653589793
el_{U} = Ellipse(H, F2_U, H5)
θ_{U}: v_{U} + H8
ρ_U: H5 * (1.0 - H6^(2.0)) / (1.0 + H6 * cos(v_{U}))
U: (ρ_U; θ_{U})
F2_N: (2.0 * I5 * I6; I8 + 180.0°)
fN(x) = x + (I26 + I6 * sin(x) - x) / (1.0 - I6 * cos(x))
u_N = Iteration(fN(x), I26, 4.0)
v_{N}: 2.0 * atan(tan(u_N / 2.0) * sqrt((1.0 + I6) / (1.0 - I6))) * 180.0° / 3.141592653589793
el_{N} = Ellipse(H, F2_N, I5)
θ_{N}: v_{N} + I8
ρ_N: I5 * (1.0 - I6^(2.0)) / (1.0 + I6 * cos(v_{N}))
N: (ρ_N; θ_{N})
θ_{J}: v_{J} + F8
J: (ρ_J; θ_{J})
C29 = 4.482882721603491 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
D29 = 1.5315684810599972 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
E29 = 1.9477092544751824 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
F29 = 1.2960272227912522 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
G29 = 3.0922413745945514 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
H29 = 6.271083892277959 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
I29 = 5.886354682075542 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
B27 = Iteration(fMe(x), B26, 4.0)
C27 = Iteration(fV(x), C26, 4.0)
D27 = Iteration(fT(x), D26, 4.0)
E27 = Iteration(fM(x), E26, 4.0)
F27 = Iteration(fJ(x), F26, 4.0)
G27 = Iteration(fS(x), G26, 4.0)
H27 = Iteration(fU(x), H26, 4.0)
I27 = Iteration(fN(x), I26, 4.0)
B28: 2.0 * atan(sqrt((1.0 + B6) / (1.0 - B6)) * tan(B27 / 2.0)) * 180.0° / 3.141592653589793
B29: B28 + B8
C28: 2.0 * atan(sqrt((1.0 + C6) / (1.0 - C6)) * tan(C27 / 2.0)) * 180.0° / 3.141592653589793
D28: 2.0 * atan(sqrt((1.0 + D6) / (1.0 - D6)) * tan(D27 / 2.0)) * 180.0° / 3.141592653589793
E28: 2.0 * atan(sqrt((1.0 + E6) / (1.0 - E6)) * tan(E27 / 2.0)) * 180.0° / 3.141592653589793
F28: 2.0 * atan(sqrt((1.0 + F6) / (1.0 - F6)) * tan(F27 / 2.0)) * 180.0° / 3.141592653589793
G28: 2.0 * atan(sqrt((1.0 + G6) / (1.0 - G6)) * tan(G27 / 2.0)) * 180.0° / 3.141592653589793
H28: 2.0 * atan(sqrt((1.0 + H6) / (1.0 - H6)) * tan(H27 / 2.0)) * 180.0° / 3.141592653589793
I28: 2.0 * atan(sqrt((1.0 + I6) / (1.0 - I6)) * tan(I27 / 2.0)) * 180.0° / 3.141592653589793
affInterMars = true