h = 6.5 (slider from -100 to 100, step 0.1)
A = (4.599999999999998, -1.7000000000000006, 0.2511275246717144)
B = (1.4243977078620738, 1.5280718980479193, 1.5822129375466725)
f = Line(A, B)
C = (1.909801560873829, 0.7047423868339227, 0.0025937146569230407)
c = Circle(f, C)
b = OrthogonalPlane(C, f)
D = Intersect(f, b)
d = Sphere(D, C)
P = (-0.7333011374179179, 2.9970496992644087, 0.6625941770104262)
g = Segment(P, D)
r = Segment(D, C)
d' = Dilate(d, h / (g^(2) - r^(2)), P)
D' = Dilate(D, h / (g^(2) - r^(2)), P)
i = OrthogonalLine(P, b)
E_1 = Intersect(i, b)
j = Segment(P, E_1)
α = Angle(D, E_1, P)
D'' = Dilate(D, h / g^(2), P)
F = Midpoint(P, D'')
a = OrthogonalPlane(F, g)
E = Intersect(i, a)
e = Sphere(E, P)
k = IntersectConic(e, d')
H' = Point(k)
l = Line(P, H')
H = Intersect(l, b)