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