S_1 = (5.72, 1.44)
B = (7.452806139534895, 5.013655906976741)
k_1 = Circle(S_1, B)
S_2 = (10.54, -1)
D = (8.26, -0.16)
k_2 = Circle(S_2, D)
C = Intersect(k_1, k_2, 2)
f = Tangent(C, k_2)
A = (8.033937127194273, 2.600998785732822)
ω = Circle(C, A)
a = 0 (slider from 0 to 6, step 1)
C_{1} = Intersect(k_2, ω, 2)
F = Intersect(k_2, ω, 1)
g = Line(C_{1}, F)
G = Intersect(k_1, ω, 2)
H = Intersect(k_1, ω, 1)
h = Line(G, H)
i, j = AngularBisector(h, g)
k, l = AngularBisector(f, h)
E = Intersect(i, k)
I = Intersect(j, l)
m = OrthogonalLine(E, f)
J = Intersect(f, m)
c = Circle(E, J)
d = Circle(I, Radius[c])
k_3 = Mirror(c, ω)
k_4 = Mirror(d, ω)
K, L = Intersect(k_3, d)
S_3 = Center(k_3)
S_4 = Center(k_4)