A = (-0.528461804082998, 0.5624850961342412)
B = (9.24, -3.22)
S_1 = (4.24, 2.62)
D = (6.426420527575388, 1.2691308945380333)
k_1 = Circle(S_1, D)
p, l = Tangent(A, k_1)
g = Line(S_1, D)
S_2 = Point(g)
k_2 = Circle(S_2, D)
h = OrthogonalLine(D, g)
i, j = AngularBisector(h, p)
S_3 = Intersect(g, i)
k_3 = Circle(S_3, D)
a = 0 (slider from 0 to 9, step 1)
C = Point(g)
c = Circle(C, D)
E = Point(i)
f = OrthogonalLine(E, p)
F = Intersect(p, f)
d = Circle(E, F)
G = Intersect(k_1, p, 1)
k = OrthogonalLine(G, p)
e = Ellipse(S_1, S_2, D)
H = Point(k)
I = Point(e)
m = Line(I, S_1)
J, K = Intersect(k_1, m)
q = Circle(I, J)
r = Circle(H, G)
S_4 = Intersect(e, k, 1)
k_4 = Circle(S_4, G)