A = (3.78, -3.88)
B = (15.02, -3.16)
f = Line(A, B)
C = Point(f)
g = OrthogonalLine(C, f)
S_1 = Point(g)
c = Circle(S_1, C)
C_{2} = Point(c)
h = Line(S_1, C_{2})
S_2 = Point(h)
d = Circle(S_2, C_{2})
a = 0 (slider from 0 to 7, step 1)
D = (13.373972856786084, -2.7425072985080496)
ω = Circle(C, D)
E = Intersect(c, ω, 2)
F = Intersect(c, ω, 1)
i = Line(E, F)
d' = Mirror(d, ω)
G = Intersect(d', i, 1)
j = OrthogonalLine(G, i)
H = Center(d')
I = Intersect(f, j)
J = Midpoint(I, G)
k = OrthogonalLine(J, j)
e = Circle(G, Segment[J, G])
K, L = Intersect(e, j)
p = Circle(H, L)
M, N = Intersect(p, k)
q = Circle(J, G)
r = Circle(N, Radius[q])
s = Circle(M, Radius[q])
l, m = Tangent(i, d')
r' = Mirror(r, ω)
q' = Mirror(q, ω)
s' = Mirror(s, ω)
m' = Mirror(m, ω)
S_6 = Center(r')
S_3 = Center(m')
S_5 = Center(q')
S_4 = Center(s')