O = Point(xAxis)
U = (3, 0)
u = Vector(O, U)
D: {{2, 0}, {0, 1.5}}
P: {{-2, 3}, {3, 1}}
evec_1 = Line((0, 0), (Element[P, 1, 1], Element[P, 2, 1]))
P_{-1} = Invert(P)
A: ((P * D) * P_{-1})
evec_2 = Line((0, 0), (Element[P, 1, 2], Element[P, 2, 2]))
f = Line(U, evec_1)
g = Line(U, evec_2)
Uj = Intersect(evec_2, f)
Ui = Intersect(g, evec_1)
Au: (A * u)
AU: (A * u)
i = Line(AU, evec_2)
j = Line(AU, evec_1)
AUj = Intersect(evec_2, j)
AUi = Intersect(i, evec_1)
uj = Segment(O, Uj)
auj = Segment(O, AUj)
ui = Segment(O, Ui)
aui = Segment(O, AUi)
dj: auj / uj
di: aui / ui
A1: 1
A2: 2
A3 = 3
A4 = 4
A5 = 5
A6 = 6
A7 = 7
A8 = 8
A9 = 9
A10 = 10
A11 = 11
A12 = 12