A = (-1.98, 2.48)
B = (-0.478788880540947, 3.1047716686018703)
f = Line(A, B)
C = (-1.3160131138583442, 0.975299638002869)
D = (0.8476232497780185, 1.0480269107301416)
g = Line(C, D)
P = (-2.1705585684037985, 1.8662087289119587)
h = Line(P, C)
i = Line(A, P)
F = Intersect(g, i)
E = Intersect(f, h)
O = Intersect(f, g)
eq1: (((x(O) - x(D))^(2) + (y(O) - y(D))^(2)) * ((x - x(C))^(2) + (y - y(C))^(2))) = (((x(O) - x(C))^(2) + (y(O) - y(C))^(2)) * ((x - x(D))^(2) + (y - y(D))^(2)))
eq2: (((x(O) - x(B))^(2) + (y(O) - y(B))^(2)) * ((x - x(A))^(2) + (y - y(A))^(2))) = (((x(O) - x(A))^(2) + (y(O) - y(A))^(2)) * ((x - x(B))^(2) + (y - y(B))^(2)))
F' = Mirror(F, eq1)
E' = Mirror(E, eq2)
j = Line(E', D)
k = Line(B, F')
Q = Intersect(j, k)