A = Point(yAxis)
B = Intersect(xAxis, yAxis)
C = Point(Circle[B, Distance[A, B], Segment[A, B]])
a, D, E, F, G, H, faceABCD, faceADHE, faceABFE, faceBCGF, faceCDHG, faceEFGH, arêteAD, arêteAB, arêteBC, arêteCD, arêteAE, arêteDH, arêteEH, arêteEF, arêteBF, arêteFG, arêteCG, arêteGH = Cube(A, B, C)
P = Midpoint(B, F)
Q = Midpoint(E, H)
d = Segment(A, Q)
M = Midpoint(G, H)
N = Midpoint(D, H)
I = Midpoint(A, B)