D = (0, 3)
f(x) = 1 + 1 / x
B = Point(xAxis)
A: (x(D), 0)
C: (x(A) - (x(B) - x(A)), 0)
c = Circle(D, A)
i = OrthogonalLine(A, xAxis)
E = Intersect(c, i, 2)
g = Segment(B, E)
h = Segment(C, E)
L_1 = Sequence(Intersect[c, Segment[(f(k), 0), E], 2], k, x(B), x(C), abs((x(B) - x(C)) / 20))
L_2 = Sequence(Segment[(k, 0), (k, f(k))], k, x(B), x(C), abs(x(B) - x(C)) / 20)
L_3 = Sequence((k, f(k)), k, x(B), x(C), abs(x(B) - x(C)) / 20)
L_4 = Sequence(Segment[L_1(k), L_3(k)], k, 1, Length[L_1])