A = Intersect(xAxis, zAxis)
B = Point(xAxis)
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)
I = Point(arêteAB)
J = Point(arêteCG)
f = Line(I, J)
K = Point(arêteAE)
L = PointIn(faceCDHG)
g = Line(K, L)
M = Point(f)
N = Point(f)
h = Segment(I, J)
i = Segment(I, N)
j = Segment(M, J)
O = Point(g)
P_{1} = Point(g)
Q = Point(arêteCG)
p = Plane(E, H, I)
k = Line(L, K)
P = Intersect(k, p)