r = 1.4 (slider from 1 to 5, step 0.01)
O: (0, 0, r)
e = Sphere(O, r)
A = (-3.3506826035236568, -0.3964474906598632, 5)
O' = Dilate(O, abs(A - O), A)
f, f_1, f_2 = Cone(O', A, (r * abs(A - O)))
g = IntersectPath(xOyPlane, f)
B, C = Focus(g)
h = Ray(A, O)
a = OrthogonalPlane(C, Line[B, C])
D = Intersect(h, a)
b = Sphere(D, abs(C - D))