Largeur de l'angle

Demi Grand Axe

Temps

L'étoile

Angle

La planette

Anomalie Vraie

Start

Stop

Montrer le cercle équivalent

Montrer la loi des aires

largeur = 5 (slider from 0 to 5, step 0.1)
a = 8.1 (slider from 0 to 20, step 0.1)
T: sqrt(a^(3)) * 365.25
excentricité = 0.73 (slider from 0 to 2, step 0.01)
f_1(x) = excentricité * sin(x)
t = 5965 (slider from 0 to T, step 1)
t': t - largeur / 100 * T
M': 2 * π * t' / T
g_2: y = x - M'
Z' = Intersect(f_1, g_2, (3.73113080250078, -0.405863350693768))
M: 2 * π * t / T
f_2: y = x - M
Z = Intersect(f_1, f_2, (3.9322552148336065, -0.5188982036780476))
Sun = (0, 0)
b: sqrt(1 - excentricité^(2)) * a
ea: excentricité * a
c: ea
O: (-c, 0)
F: (-2 * c, 0)
Orbit = Ellipse(Sun, F, a)
OC = Circle(O, a)
Aphélie, Périhelie = Intersect(Orbit, xAxis)
E_1: x(Z)
P = Rotate(Périhelie, E_1, O)
E = Angle(Périhelie, O, P)
d = Segment(O, P)
ll = OrthogonalLine(P, xAxis)
X = Intersect(ll, xAxis)
h = Segment(X, P)
Planet = Intersect(Orbit, h)
i = Segment(Sun, Planet)
θ = Angle(Périhelie, Sun, Planet)
cercle = false
E_2: x(Z')
P' = Rotate(Périhelie, E_2, O)
α = Angle(Périhelie, O, P')
ll' = OrthogonalLine(P', xAxis)
X' = Intersect(ll', xAxis)
h' = Segment(P', X')
A = Intersect(Orbit, h', 1)
bestrekenperk, planet, sun, a_1 = Polygon(Sun, A, Planet)
ellipse = true
Vitesse = 0.84 (slider from 0 to 2, step 0.01)