milieu

A = Intersect(xAxis, yAxis)
c = Circle(A, 1)
B = Point(c)
n = 4 (slider from 3 to 100, step 1)
a: 360 / n
l: (2 * sin(π / n))
d = Circle(B, l)
C = Intersect(c, d, 1)
b, D, E, F, G, H, I, faceBCDE, faceBEIF, faceBCGF, faceCDHG, faceDEIH, faceFGHI, arêteBE, arêteBC, arêteCD, arêteDE, arêteBF, arêteEI, arêteFI, arêteFG, arêteCG, arêteGH, arêteDH, arêteHI = Cube(B, C, xOyPlane)
J = Midpoint(arêteFI)
K = Midpoint(arêteFG)
L = Midpoint(arêteGH)
M = Midpoint(arêteHI)
N = Midpoint(arêteDH)
O = Midpoint(arêteEI)
P = Midpoint(arêteBE)
Q = Midpoint(arêteBF)
R = Midpoint(arêteBC)
S = Midpoint(arêteCG)
T = Midpoint(arêteDE)
U = Midpoint(arêteCD)
e = true
V = (-0.27402698770496414, 0.8557992405381025, 0)
W = Intersect(arêteGH, arêteFG)
poly1, i, t, c_1, q = Polygon(I, T, C, Q)