P = (0, 0, 1)
Q = (1, 2, 2)
g = Line(Q, P)
A = Intersect(g, a)
u = OrthogonalVector(a)
B: A - u * 4
b = Plane(B, a)
C = Intersect(g, b)
v = Vector(C, A)
D: A + v
c = Sphere(C, 4)
d = Sphere(D, 4)
e = OrthogonalLine(C, a)
E, F = Intersect(e, c)
f = OrthogonalLine(D, a)
G, H = Intersect(f, d)