O = Intersect(a, b)
P = (1.338, 1.758)
d = OrthogonalLine(P, a)
H = Intersect(a, d)
I = Point(a)
c = Line(P, I)
A_1 = Midpoint(O, I)
e = Circle(A_1, O)
M = Intersect(e, d, 2)
f = Circle(O, M)
g = Segment(M, O)
h = Segment(M, I)
i = OrthogonalLine(P, d)
P_1 = Intersect(i, g)
k = Circle(O, P_1)
B = Intersect(k, b, 1)
A = Intersect(f, b, 1)
j = Segment(P, P_1)
l = OrthogonalLine(B, b)
C = Intersect(f, l, 1)
m = OrthogonalLine(C, a)
F = Intersect(a, m)
F' = Mirror(F, O)
p = Ellipse(F, F', B)