Eléments orbitaux

Grand axe

Excentricité

Inclinaison

Noeud ascendant

Périastre

Anomalie vraie

Ω

ω

M

V = 0 rad (slider from 0 to 180*°, step 1*°)
a = 10 (slider from 5 to 15, step 1)
e = 0.01 (slider from 0.01 to 0.95, step 0.01)
k: (2 * e) - e^(3) / 4
a_V: V + (k * sin(V))
b: (a * sqrt(1 - (e * e)))
c: (a * e)
F = (0, 0, 0)
VV = (0, 4, 0)
U = (4, 0, 3.552713678800501E-15)
ω = 0 rad (slider from 0*° to 6.283185307179586, step 1*°)
Ω = 1.7453292519943295 rad (slider from 0*° to 6.283185307179586, step 1*°)
i = 0.017453292519943295 rad (slider from 0 to 45*°, step 0.5*°)
h = Circle(F, 8, xOyPlane)
M1: {{(cos(i) * cos(Ω)), (-sin(Ω)), (cos(Ω) * sin(i))}, {(cos(i) * sin(Ω)), cos(Ω), (sin(i) * sin(Ω))}, {(-sin(i)), 0, cos(i)}}
U' = ApplyMatrix(M1, U)
VV' = ApplyMatrix(M1, VV)
noe = Line(F, VV')
N, N' = Intersect(noe, h)
orb = Plane(F, U', VV')
per = OrthogonalLine(F, orb)
U'' = Rotate(U', ω, per)
VV'' = Rotate(VV', ω, per)
aps = Line(F, VV'')
bb = Line(F, U'')
s_C = Sphere(F, c)
C, C' = Intersect(aps, s_C)
F' = Mirror(F, C)
s_B = Sphere(F, b)
B', B = Intersect(bb, s_B)
d = Ellipse(F, F', B)
bb' = Segment(B, B')
Y = (0, 14, 0)
u = Vector(F, Y)
Ω' = Angle(Y, F, N)
P, A = Intersect(aps, d)
P' = Rotate(P, a_V, per)
f = Ray(F, P')
S = Intersect(f, d)
ω' = Angle(N, F, P)
V' = Angle(P, F, S)
g = Segment(F, S)