Author: Daniel Mentrard

R = 7.5 (slider from 0.1 to 10, step 0.1)
d = 3.1 (slider from -5 to 5, step 0.1)
S = (6.889305968514359, 0.41565990876470266)
Zoom = 1 (slider from 0.1 to 5, step 0.1)
View3d = (-4.389828957325238, 1.1039535754990089)
ya: y(View3d)
costheta: cos(ya)
sintheta: sin(ya)
xa: x(View3d)
cosphi: cos(xa)
xy: -(Zoom) * cosphi * sintheta
sinphi: sin(xa)
B_1 = (-3.8125503981287236, -6.511682801593476)
Z: (x(B_1), y(B_1) + Zoom * costheta)
w = Vector(B_1, Z)
X: (x(B_1) + Zoom * sinphi, y(B_1) + xy)
u = Vector(B_1, X)
g = 1 (slider from -5 to 5, step 0.1)
k = 1 (slider from -5 to 5, step 0.1)
a = OrthogonalLine(View3d, yAxis)
b = OrthogonalLine(View3d, xAxis)
c: cos(ya)
Y: (x(B_1) + Zoom * cosphi, y(B_1) + Zoom * sinphi * sintheta)
v = Vector(B_1, Y)
liste2 = Sequence((x(S), y(S)) + (R * cos(i) + cos(i) + cos(d)) * u + (R * sin(i) + sin(i) + cos(d)) * w, i, 0.0, 6.28319, 0.1)
e = 42.23 (slider from 0 to 62.83, step 0.01)
f = Conic(Element[liste2, 1.0], Element[liste2, 5.0], Element[liste2, 10.0], Element[liste2, 15.0], Element[liste2, 20.0])
liste1 = Sequence((x(S), y(S)) + (R * cos(i / 10.0) + cos(i / 10.0) + cos(d)) * u + (R * sin(i / 10.0) + sin(i / 10.0) + cos(d)) * w, i, e, e, 0.1)
r = 0.751 (slider from 0.001 to 5, step 0.001)
liste3 = Sequence((x(Element[liste1, 1.0]), y(Element[liste1, 1.0])) + (r / 10.0 * cos(i) + cos(i)) * u + (r / 10.0 * sin(i) + sin(i)) * w, i, 0.0, 6.28319, 0.1)
h = Conic(Element[liste3, 1.0], Element[liste3, 5.0], Element[liste3, 10.0], Element[liste3, 15.0], Element[liste3, 20.0])
lune = 5.56 (slider from 0 to 6.24, step 0.01)
n = 16.8 (slider from 1 to 100, step 0.1)
liste4 = Sequence((x(Element[liste1, 1.0]), y(Element[liste1, 1.0])) + (r / 10.0 * cos(n * i) + cos(n * i)) * u + (r / 10.0 * sin(n * i) + sin(n * i)) * w, i, lune, lune, 0.1)
Speed = 1.3 (slider from 0 to 5, step 0.1)
p = Circle(S, Zoom)