A = Intersect(zAxis, xOyPlane)
B = (4, 0, 1.7763568394002505E-15)
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(a)
J = Point(a)
K = PointIn(faceABCD)
b = Plane(I, J, K)
poly1, L, M, N, O, P, Q, c, d, e, f, g, h = IntersectPath(b, a)