A: (5, 5, 10)
B: (5, 5, 0)
C: (5, 0, 0)
D: (0, 5, 0)
E: (0, 0, 0)
q1, f, g, h, i = Polygon(E, D, B, C)
I = Point(zAxis)
a, F, G, H, faceDEIF, faceBDFG, faceBCHG, faceCEIH, faceFGHI, edgeEI, edgeDF, edgeBG, edgeCH, edgeFI, edgeFG, edgeGH, edgeHI = Prism(q1, I)
j = Segment(E, G)
J = Point(zAxis)
K = Point(yAxis)
L = Point(yAxis)
b = Plane(J, K, L)
M = Point(xAxis)
N = Point(xAxis)
c = Plane(M, N, J)
P = PointIn(c)
d = Plane(P, xOyPlane)
e = Circle((0, 0, z(P)), P, space)
Q: (0, 0, z(P))
k = Circle(zAxis, Q)
p = Circle(zAxis, P)
O, R = Intersect(b, p)
l = OrthogonalLine(O, b)
m = OrthogonalLine(P, c)
S = Intersect(m, l)
n = Segment(Q, S)
T = Intersect(n, p, 1)
o = Sphere(S, P)