variables orbites
Author: aaa

Date origine

Date

1 ua

Mars

a

long. 0

γ

a en ua

eccentricity

Earth

lg périhelion

Year

Day

min

h

year périh

Month

month péri

day péri

h péri

min péri

julian days

annual period in Earth days

B8 = 1.796589572125403 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
B21 = 20
B20 = 15
B19 = 4
B18 = 1
B17 = 2009
jjp_T: 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
B4 = 1.49598034E8
B11 = 5
B12 = 22
B13 = 0
B14 = 0
B10 = 2011
jj0: 367.0 * B10 - floor(7.0 * (B10 + floor((B11 + 9.0) / 12.0)) / 4.0) - floor(3.0 * floor(floor(B10 + floor((B11 - 9.0) / 7.0)) / 100.0 + 1.0) / 4.0) + floor(275.0 * B11 / 9.0) + B12 + 1721029.0 - 0.5 + B13 / 24.0 + B14 / 24.0 / 60.0
t = 3 (slider from -300 to 6000, step 1.0E-4)
jj: jj0 + t
B1 = 1.4959787E8
B7 = 365.256
B16 = 5.388704065532492 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
C4 = 2.27939184E8
C8 = 5.865353484252144 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
C7 = 686.98
C17 = 2009
C19 = 21
C18 = 4
C20 = 9
C21 = 46
jjp_M: 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
C16 = 0.6614797865058508 rad (slider from 0 to 6.283185307179586, step 0.017453292519943295)
alp: floor((floor(jj + 0.5) - 1867216.25) / 36524.25)
a: floor(jj + 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)
hr: (jj + 0.5 - floor(jj + 0.5)) * 24.0
ht: floor(hr)
mt: floor((hr - ht) * 60.0)
jjs = Mod(floor(jj - 0.5) + 0.5 + 1.5, 7.0)
C13: ht
C14: mt
e: floor((b - d) / 30.6001)
jr: b - d - floor(30.6001 * e) + Mod[jj + 0.5, 1.0]
mr: e - If[e < 13.5, 1.0, 13.0]
yr: c - If[mr > 2.5, 4716.0, 4715.0]
C12: floor(jr)
C11: mr
C10: yr
B6 = 0.017
C6 = 0.093
C1 = 0
B5: B4 / $B$1
CB51 = 0
C5: C4 / $B$1
B22: jjp_T
C22: jjp_M
B24: jj
AM_T: (B24 - B22) * 360.0 / B7
u0_T: 3.141592653589793 / 180.0 * AM_T
u1_T: u0_T + (u0_T + B6 * sin(u0_T) - u0_T) / (1.0 - B6 * cos(u0_T))
u2_T: u1_T + (u0_T + B6 * sin(u1_T) - u1_T) / (1.0 - B6 * cos(u1_T))
u3_T: u2_T + (u0_T + B6 * sin(u2_T) - u2_T) / (1.0 - B6 * cos(u2_T))
u4_T: u3_T + (u0_T + B6 * sin(u3_T) - u3_T) / (1.0 - B6 * cos(u3_T))
u_T: u4_T + (u0_T + B6 * sin(u4_T) - u4_T) / (1.0 - B6 * cos(u4_T))
v_T: 2.0 * atan(tan(u_T / 2.0) * sqrt((1.0 + B6) / (1.0 - B6))) * 180.0° / 3.141592653589793
ρ_T: B5 * (1.0 - B6^(2.0)) / (1.0 + B6 * cos(v_T))
θ_T: v_T + B8
E: (ρ_T; θ_T)
S = (0, 0)
F2_T: (2.0 * B5 * B6; B8 + 180.0°)
el_T = Ellipse(S, F2_T, B5)
u0_M: 2.0 * 3.141592653589793 * (B24 - C22) / C7
fM(x) = x + (u0_M + C6 * sin(x) - x) / (1.0 - C6 * cos(x))
u_M = Iteration(fM(x), u0_M, 4.0)
v_M: 2.0 * atan(tan(u_M / 2.0) * sqrt((1.0 + C6) / (1.0 - C6))) * 180.0° / 3.141592653589793
θ_M: v_M + C8
ρ_M: C5 * (1.0 - C6^(2.0)) / (1.0 + C6 * cos(v_M))
M: (ρ_M; θ_M)
F2_M: (2.0 * C5 * C6; C8 + 180.0°)
el_M = Ellipse(S, F2_M, C5)
texte1: "" + (If[jjs ≟ 0.0, "Sunday", If[jjs ≟ 1.0, "Monday", If[jjs ≟ 2.0, "Tuesday", If[jjs ≟ 3.0, "Wednesday", If[jjs ≟ 4.0, "Thursday", If[jjs ≟ 5.0, "Friday", "Saturday"]]]]]]) + " " + (floor(jr)) + "/" + mr + "/" + yr + "    " + C13 + "h" + C14 + "min"

Tool: ekep
Input: C5, um0, um1, um2, um3, um4
Output: um

um1: um0 + (um0 + C5 * sin(um0) - um0) / (1.0 - C5 * cos(um0))
um2: um0 + (um1 + C5 * sin(um1) - um1) / (1.0 - C5 * cos(um1))
um3: um0 + (um2 + C5 * sin(um2) - um2) / (1.0 - C5 * cos(um2))
um: um0 + (um4 + C5 * sin(um4) - um4) / (1.0 - C5 * cos(um4))