a(ε)

a

x_0

x_0-δ_2

x_0+δ_1

x0 = 2
f1(x) = 1 / (x - x0)^(2)
p1 = Invert(f1)
a = 2.6 (slider from 1 to 4, step 0.1)
x1: p1(a)
p2(x) = (-sqrt(1 / x)) + 2
x2: p2(a)
F = Point({x0, 0})
A = Point({0, a})
b = OrthogonalLine(A, yAxis)
B = Intersect(f1, b, (1.37982632872054, 2.6000000139738715))
c = OrthogonalLine(B, b)
C = Point({x2, 0})
d = Segment(B, C)
D = Intersect(f1, b, (1.37982632872054, 2.6000000139738715))
e = OrthogonalLine(D, xAxis)
E = Point({x1, 0})
G = Point(yAxis)
h = Segment(A, G)
f2(x) = If[x < x0, 1 / (x - x0)^(2)]
H = Intersect(f1, b, (2.62017367127946, 2.6000000139738737))
g = Segment(H, E)