X_0

Y_0

X

Y

Q-D věta

A_{1} = Point(xAxis)
V = Point(zAxis)
a, d, b = Cone(A_{1}, V, 1)
D_{2} = Point(zAxis)
e, f, g = Cone(D_{2}, V, 1)
h = Line(A_{1}, D_{2})
E = (1.4342215073760907, -0.8174607476893083, 0)
i = Line(E, yAxis)
F = Point(zAxis)
p = Plane(F, i)
c = IntersectPath(p, a)
G, H, I, J = Corner(c)
k = IntersectPath(p, e)
K, L, M_{1}, N = Corner(k)
q: y = 0
j = IntersectPath(q, p)
r = IntersectPath(q, a)
O, P, Q, R = Corner(r)
s = IntersectPath(q, e)
S, T, U, V_{1} = Corner(s)
l = Line(V, P)
m = Line(V, O)
n, t = AngularBisector(l, j)
S_1 = Intersect(h, t)
f_1 = OrthogonalLine(S_1, p)
F_1 = Intersect(j, f_1)
o = Sphere(S_1, F_1)
S_2 = Intersect(h, n)
g_1 = OrthogonalLine(S_2, p)
F_2 = Intersect(j, g_1)
u = Sphere(S_2, F_2)
M = Point(c)
h_1 = Line(V, M)
A = Intersect(l, j)
B = Intersect(m, j)
C, W = Intersect(m, o)
Z, B_1 = Intersect(m, u)
p_1 = OrthogonalPlane(C, h)
c_1 = IntersectPath(p_1, o)
q_1 = OrthogonalPlane(Z, h)
d_1 = IntersectPath(q_1, u)
C_1 = Intersect(h_1, c_1)
D_1 = Intersect(h_1, d_1)
i_1 = IntersectPath(p, p_1)
D = Intersect(i_1, j)
v = false