B: (9, 0, 0)
D: (0, 9, 0)
A = Intersect(yAxis, zAxis)
C_{1} = Point(Circle[A, Distance[D, A], Segment[D, A]])
a, C, H_{1}, E, F, G, faceACFD, faceDFJG, faceADGH, faceACIH, faceCFJI, faceGHIJ, arêteDF, arêteAD, arêteAC, arêteCF, arêteDG, arêteFJ, arêteGJ, arêteGH, arêteAH, arêteHI, arêteCI, arêteIJ = Cube(D, A, C_{1})
P: (1, 0, 0)
R: (9, 3, 9)
Q: (0, 0, 3)
O = Midpoint(E, C)
p = Plane(R, Q, P)
f = OrthogonalLine(O, p)
H = Intersect(f, p)