α = 0 rad (slider from 0*° to 6.283185307179586, step 2*°)
E_1 = (-13.2424198724697, 4.470084464094601)
O = (-13.107488695473327, 4.784923877086137)
u, v, w = Repère1(O, E_1)
E' = Point3D(u, v, w, {2 * cos(α) + 2 * sin(α), 2 * sin(α) - 2 * cos(α), 2})
A = Point3D(u, v, w, {2 * cos(α) + 2 * sin(α), 2 * sin(α) - 2 * cos(α), -2})
E'' = Mirror(E', A)
B = Point3D(u, v, w, {2 * cos(α) - 2 * sin(α), 2 * sin(α) + 2 * cos(α), -2})
Z: A + (B + E' - 2 * A) / sqrt(2)
W_1: 2 * A + 4 * w - B
B' = Mirror(B, A)
c = Conic(E'', B, Z, E', B')
poly1, a, e', w_1, b' = Polygon(A, E', W_1, B')
W = PointIn(poly1)
b = Line(A, W)
E_b = Intersect(c, b, 1)
cos1 = If(Defined[E_b], x(E_b - A) / x(A - B), 0)
sin1 = If(1 - cos1^(2) > 0, sqrt(1 - cos1^(2)), 0)
E: A + cos1 * (A - B) + 4 * sin1 * w
c' = Translate(c, E - A)
H': E - A + B
I: H' + 4 * w
G': E + A - B
J: G' + 4 * w
poly2, g, h, i, j = Polygon(G', H', I, J)
K = PointIn(poly2)
d = Line(E, K)
F_b = Intersect(c', d, 1)
cos2 = If(Defined[F_b], x(F_b - E) / x(A - B), -1)
sin2 = If(1 - cos2^(2) > 0, sqrt(1 - cos2^(2)), 0)
C = Point3D(u, v, w, {-2 * cos(α) - 2 * sin(α), -2 * sin(α) + 2 * cos(α), -2})
D = Point3D(u, v, w, {-2 * cos(α) + 2 * sin(α), -2 * sin(α) - 2 * cos(α), -2})
F: E + cos2 * (A - B) + 4 * sin2 * w
G: F + D - A
H: E + D - A
e = Segment(A, B)
f = Segment(B, C)
k = Segment(C, D)
l = Segment(D, A)
m = Segment(A, E)
n = Segment(E, H)
p = Segment(E, F)
q = Segment(H, G)
r = Segment(G, F)
s = Segment(H, D)
L: A + cos1 * (A - D) + 4 * sin1 * w
M: B + cos1 * (A - D) + 4 * sin1 * w
N: B + cos1 * (B - A) + 4 * sin1 * w
P: C + cos1 * (B - A) + 4 * sin1 * w
Q: C + cos1 * (D - A) + 4 * sin1 * w
R: D + cos1 * (D - A) + 4 * sin1 * w
t = Segment(A, L)
a_1 = Segment(L, M)
b_1 = Segment(M, B)
c_1 = Segment(B, N)
d_1 = Segment(N, P)
e_1 = Segment(P, C)
f_1 = Segment(C, Q)
g_1 = Segment(Q, R)
h_1 = Segment(R, D)
poly6, a_3, b_4, m_1, l_1 = Polygon(A, B, M, L)
poly3, c_2, b_2, n_1, p_1 = Polygon(C, B, N, P)
poly4, b_3, a_2, d_2, c_3 = Polygon(B, A, D, C)
poly5, c_4, q_1, r_1, d_3 = Polygon(C, Q, R, D)
poly7, d_4, h_2, e_2, a_4 = Polygon(D, H, E, A)
poly8, b_5, a_5, l_2, m_2 = Polygon(B, A, L, M)
poly9, c_5, d_5, r_2, q_2 = Polygon(C, D, R, Q)
poly10, e_3, h_3, g_2, f_2 = Polygon(E, H, G, F)
poly11, f_3, e_4, h_4, g_3 = Polygon(F, E, H, G)
poly12, d_6, c_6, b_6, a_6 = Polygon(D, C, B, A)
poly13, g_4, h_5, e_5, f_4 = Polygon(G, H, E, F)
poly14, e_6, a_7, d_7, h_6 = Polygon(E, A, D, H)
poly15, p_2, n_2, b_7, c_7 = Polygon(P, N, B, C)
poly16, g_5, h_7, e_7, f_5 = Polygon(G, H, E, F)

Tool: Repère1
Origine, E_1
Input: O, A
Output: u, v, w

E: O - (1, 0)
E' = Rotate(E, 5*°, O)
poly1, b, d, e, f, g, h, i, j, k, l, m, n, p, q, r, s, t, a_1, b_1, c_1, d_1, e_1, f_1, g_1, h_1, i_1, j_1, k_1, l_1, m_1, n_1, p_1, q_1, r_1, s_1, t_1, a_2, b_2, c_2, d_2, e_2, f_2, g_2, h_2, i_2, j_2, k_2, l_2, m_2, n_2, p_2, q_2, r_2, s_2, t_2, a_3, b_3, c_3, d_3, e_3, f_3, g_3, h_3, i_3, j_3, k_3, l_3, m_3, n_3, p_3, q_3, r_3, F, G, H, I, J, K, L, M, N, P, Q, R, S, T, U, V, W, Z, A_1, B_1, C_1, D_1, E_1, F_1, G_1, H_1, I_1, J_1, K_1, L_1, M_1, N_1, O_1, P_1, Q_1, R_1, S_1, T_1, U_1, V_1, W_1, Z_1, A_2, B_2, C_2, D_2, E_2, F_2, G_2, H_2, I_2, J_2, K_2, L_2, M_2, N_2, O_2, P_2, Q_2, R_2, S_2, T_2, U_2, V_2, W_2, Z_2, A_3, B_3, C_3, D_3 = Polygon(E, E', 72)
c = Circle(O, 1)
A = PointIn(poly1)
a: x = x(A)
B = Intersect(c, a, 2)
B' = Rotate(B, 90*°, O)
C: (x(B'), y(O) - x(A - O) * y(A - O) / x(B' - O))
D: (x(O), y(O) + sqrt(1 - y(A - O)^(2) - y(C - O)^(2)))
u = Vector(O, A)
v = Vector(O, C)
w = Vector(O, D)

Tool: Point3D
u, v, w, liste des coordonnées
Input: u, v, w, liste1
Output: A

liste1: {-2, -2, -2}
A: (0, 0) + Element[liste1, 1] * u + Element[liste1, 2] * v + Element[liste1, 3] * w

Tool: ADroite1 (command blablabvlal)
Input: A, B, C
Output: a

a: (x(A) - x(B)) * y(C) - (y(A) - y(B)) * x(C) > (x(A) - x(B)) * y(A) - (y(A) - y(B)) * x(A)

Tool: ADroite
Input: A, B, C
Output: a

a: (x(A) - x(B)) * y(C) - (y(A) - y(B)) * x(C) ≥ (x(A) - x(B)) * y(A) - (y(A) - y(B)) * x(A)