四棱锥外接球

A = Point(yAxis)
B = Point(xAxis)
C = Intersect(xAxis, zAxis)
D = Point(zAxis)
a, faceABC, faceACD, faceABD, faceBCD, edgeAC, edgeAB, edgeBC, edgeCD, edgeAD, edgeBD = Pyramid(C, A, B, D)
b = Plane(edgeAD, edgeBD)
f = OrthogonalLine(C, b)
g = Line(A, edgeBC)
h = Line(B, edgeAC)
F = Intersect(g, h)
i = Line(D, edgeBC)
j = Line(D, edgeAC)
k = OrthogonalLine(A, j, space)
G = Intersect(k, j)
l = Line(G, i)
m = OrthogonalLine(B, i, space)
H = Intersect(i, m)
n = OrthogonalLine(H, l, space)
I = Intersect(l, n)
J = Intersect(n, m)
K = Intersect(m, l)
L = Intersect(i, j)
M = Intersect(k, h)
q1, g_1, i_1, j_1, d = Polygon(G, I, J, D)
q2, f_1, a_1, c, b_1 = Polygon(F, A, C, B)
q3, j_2, b_2, c_1, d_1 = Polygon(J, B, C, D)
q4, i_2, j_3, b_3, f_2 = Polygon(I, J, B, F)
p = Segment(I, C)
N = Midpoint(I, C)
q5, d_2, g_2, a_2, c_2 = Polygon(D, G, A, C)
q6, g_3, i_3, f_3, a_3 = Polygon(G, I, F, A)
e = Sphere(N, C)