Paramètres d'une orbite

Plan écliptique

Plan Orbite

Noeuds

Orbite 1

Orbite 2

Orbite 3

Orbite

Attention, le centre du repère n'est pas le centre de l'ellipse, 
mais l'un des foyer où est placé le corps principal (Soleil, 
Terre... suivant le système étudié).

     direction
point vernal (γ)

Normale au
plan orbital

a = 2.9 (slider from 0.1 to 10, step 0.1)
e = 0.28 (slider from 0 to 1, step 0.01)
i = 0 (slider from 0 to 180, step 0.1)
Ω = 0 (slider from 0 to 360, step 0.1)
ω = 0 (slider from 0 to 360, step 0.1)
c: (a * e)
F': ((-2 * c), 0)
F: (0, 0)
orbite1 = Ellipse(F, F', a)
orbite4 = Rotate(Rotate[Rotate[Ellipse[F, F', a], ω°, zAxis], ((-i))°, yAxis], Ω°, zAxis)
fi = If(i > 0, true, false)
fΩ = If(Ω ≟ 0, false, true)
fω = If(ω ≟ 0, false, true)
Pecl = Plane(xAxis, yAxis)
fplecl = false
fplanorb = false
vorb = Rotate(Vector[Rotate[Rotate[Vector[(0, 0, 0), (0, 0, 1 + a / 2)], ω°, zAxis], ((-i))°, yAxis]], Ω°, zAxis)
Porb = OrthogonalPlane(F, vorb)
lignd = IntersectPath(Pecl, Porb)
N_1, N_2 = Intersect(lignd, orbite4)
fnoeuds = false
forb1 = false
forb2 = false
forb3 = false
forbe = false
F'' = Rotate(Rotate[Rotate[F', ω°, zAxis], ((-i))°, yAxis], Ω°, zAxis)
orbe = Ellipse(F'', F, a, Porb)
sgaps = Rotate(Rotate[Rotate[Segment[(a - c, 0), ((-(a + c)), 0)], ω°, zAxis], ((-i))°, yAxis], Ω°, zAxis)
Périhélie = Rotate(Rotate[Rotate[(a - c, 0), ω°, zAxis], ((-i))°, yAxis], Ω°, zAxis)
Aphélie = Rotate(Rotate[Rotate[((-(a + c)), 0), ω°, zAxis], ((-i))°, yAxis], Ω°, zAxis)
b: sqrt(a^(2) - c^(2))
orbite2 = Rotate(Ellipse[F, F', a], ((-i))°, yAxis)
orbite3b = Rotate(orbite2, Ω°, zAxis)
orbite3 = Rotate(Rotate[Ellipse[F, F', a], ((-i))°, yAxis], Ω°, zAxis)
orbite4b = Rotate(orbite3b, ω°, Line[F, vorb])
orbite4c = Rotate(Rotate[Rotate[Ellipse[F, F', a], ((-i))°, yAxis], Ω°, zAxis], ω°, Line[F, vorb])
C: ((-c), 0)
texte2: "Demi-grand axe   " + a + "

excentricité           " + e + "

Inclinaison             " + i + "°

Long. noeud asc. " + Ω + "°

Arg. périhélie        " + ω + "°"