$d$

$O$

$B$

2c=δ±\sqrt{\delta^2-4b^2}

$b$

$\delta$

$F'_1$

$F_1$

$F'_2$

$F_2$

A = (-1.34, 5.02)
B_1 = (-1.46, -2.16)
d = Line(A, B_1)
O = (4.04, -0.34)
b_1 = Line(O, d)
B = Point(b_1)
a = OrthogonalLine(O, d)
b = Segment(B, O)
c = OrthogonalLine(B, d)
C = Intersect(c, d)
δ = Segment(B, C)
c1: (δ + sqrt(δ^(2) - 4 * b^(2))) / 2
e = Circle(O, c1)
F'_1, F_1 = Intersect(e, a)
f = Ellipse(F_1, F'_1, B)
c2: (δ - sqrt(δ^(2) - 4 * b^(2))) / 2
g = Circle(O, c2)
F'_2, F_2 = Intersect(g, a)
h = Ellipse(F_2, F'_2, B)
k = Circle(O, δ / 2)
D, E = Intersect(k, b_1)