Hauteur = 10 (slider from 0 to 10, step 0.1)
Rayon = 5 (slider from 0 to 5, step 0.1)
f(x) = sqrt(Rayon² - x²) / 4 - 8
A = Point(yAxis)
g(x) = -f(x) - 16
B: (Rayon, -8)
C: (-Rayon, -8)
D = Point(g)
E = Point(g)
F = Point(f)
G = Point(f)
c = Conic(D, E, F, G, C)
d = Segment(A, B)
e = Circle(A, Hauteur)
H = Intersect(e, yAxis, 2)
h = Segment(A, H)
u = Vector(A, H)
k = Translate(c, u)
I = Translate(C, u)
J = Translate(B, u)
i = Segment(C, I)
j = Segment(B, J)
K = (31.075000000000003, 11.274999999999995)
p = Circle(K, Rayon)
L = Point(p)
l = Segment(K, L)
m = OrthogonalLine(L, l)
n: 2 π Rayon
q = Circle(L, n / 2)
M = Intersect(q, m, 2)
N = Intersect(q, m, 1)
r = OrthogonalLine(M, m)
s = OrthogonalLine(N, m)
t = Circle(M, Hauteur)
c_1 = Circle(N, Hauteur)
O = Intersect(t, r, 2)
P = Intersect(c_1, s, 2)
Q, m_1, n_1, p_1, o = Polygon(M, N, P, O)
R = Midpoint(O, M)
a_1 = OrthogonalLine(R, o)
d_1 = Mirror(p, a_1)