Sam Loyd

B = Intersect(xAxis, yAxis)
D: (2, 2.5)
A: (0, 2.5)
C = Point(xAxis)
ox: y = 0
top: y = 5
med: y = 2.5
A' = Rotate(A, 90°, B)
c = Semicircle(A', C)
F = Intersect(c, yAxis, 1)
CF = Segment(F, C)
v_{FA} = Vector(F, A)
C' = Translate(C, v_{FA})
ex = Line(C', ox)
E = Intersect(ex, CF)
EC' = Segment(E, C')
AB = Segment(A, B)
CD = Segment(D, C)
M = Midpoint(A, D)
AM = Segment(A, M)
MD = Segment(M, D)
W = Point(AM)
WW: W
u = Vector(A, M)
I_0 = Translate(W, Vector[(2.5 * u)])
Ix = Ray(I_0, u)
I = Point(Ix)
c_I = Circle(I, W)
II: I
Q = Intersect(c_I, top, 1)
P = Intersect(c_I, ox, 1)
d = CircleArc(I, Q, P)
P' = Translate(P, Vector[(-u)])
P'' = Translate(P, Vector[(-2 * u)])
P''' = Translate(P, Vector[(-3 * u)])
l1 = Sequence(Translate[d, Vector[((k - 5) * u)]], k, 1, 8)
I_L: W + (sqrt(26) / 2 * u)
W_L: A + (0.8 * u)