A = Intersect(zAxis, xAxis)
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êteEF)
J = Point(arêteEH)
K = Point(arêteBF)
p = Plane(J, K, I)