P = O

D = P_1

F=P_3

O = Intersect(xAxis, yAxis)
D = (3.5426881220300985, 5.222300450211408)
P_2 = (5.041709074458816, 3.3185954871990995)
F = (3.886836694557285, 1.790087925564719)
P_4 = (0.9782372845234504, 3.011813684383035)
P_5 = (0.8298215712885249, 4.728219127373029)
c = Conic(D, P_2, F, P_4, P_5)
ox: y = 0
oy: x = 0
od = Line(O, D)
of = Line(O, F)
df = Line(D, F)
cd = Line(D, ox)
ef = Line(F, ox)
w = Polar(O, c)
A = Intersect(od, w)
C' = Mirror(D, Circle[Midpoint[O, A], O])
B_{1} = Intersect(w, of)
E' = Mirror(F, Circle[Midpoint[O, B_{1}], O])
c'e' = Line(C', E')
I = Intersect(df, c'e')
fc' = Line(C', F)
de' = Line(D, E')
J = Intersect(fc', de')
ij = Line(I, J)
U = Intersect(ox, w)
nn = Diameter(ox, c)
N_1 = Intersect(cd, nn)
N_2 = Intersect(nn, ef)
C = Mirror(D, N_1)
E = Mirror(F, N_2)
M = Intersect(nn, ox)
my = OrthogonalLine(M, ox)
c_{OU} = Circle(Midpoint[O, U], O)
b', b = Intersect(c_{OU}, my)
Mb = Segment(M, b)
Mb' = Segment(M, b')
tm, tm' = Tangent(M, c_{OU})
T = Intersect(c_{OU}, tm)
c_m = Circle(M, T)
B, B' = Intersect(c_m, ox)
MT = Segment(M, T)