0

0

1

1

x

y

x

y

A = Point(xAxis)
A' = Mirror(A, yAxis)
a = OrthogonalLine(A, xAxis)
b = OrthogonalLine(A', xAxis)
B = Point(a)
c = Ray(A, B)
A'_1 = Mirror(A, yAxis)
C = Point(c)
C' = Mirror(C, yAxis)
f(x) = 1 / x
D = Point(a)
in: y(D) / y(C)
out: f(in)
E = Point({x(A'), (out * y(C'))})
d = Segment(D, E)
F = Point({in, out})
e = OrthogonalLine(F, xAxis)
g = OrthogonalLine(F, yAxis)
G = Intersect(e, xAxis)
h = Segment(F, G)
H = Intersect(g, yAxis)
i = Segment(F, H)
k: x^(2) / (x(A))^(2) + y^(2) / (y(C))^(2) = 1