Edit-2
Manipulate[
Module[{F, o, P, AB},
F = Rationalize[{0.4, 0.5}, 0];
o = {0, 0};
P[t_] = AngleVector@t;
AB[t_] :=
AngleVector@s /.
NSolve[{(AngleVector@s - F) . (AngleVector@s -
F) == (AngleVector@s - P@t) . (AngleVector@s - P@t),
0 <= s <= 2 π}, s];
Show[Table[
Graphics[{{Hue[Rescale[t, {0, 2 π}]],
Line[AB@t]}}], {t, .1, τ, .1}],
Graphics[{Circle[], Point[{o, F}], Line[{F, P@τ}],
Line[AB@τ]}], PlotRange -> 1.2]], {τ, 0, 2 π}]
Edit-1
Use the same setting of the questionor.
Clear[line, plot, n, lines];
line[t_] := Module[{O, P, A, B, F}, O = {0, 0};
P = AngleVector[t];
F = {0.4, 0.5};
{A, B} =
RegionIntersection[Circle[], PerpendicularBisector[{P, F}]][[1]];
Graphics[{Hue[Rescale[t, {0, 2 π}]],
Line@{A, B}}, PlotRange -> 1.2]];
plot[t_] := Module[{O, P, A, B, F}, O = {0, 0};
P = AngleVector[t];
F = {0.4, 0.5};
{A, B} =
RegionIntersection[Circle[], PerpendicularBisector[{P, F}]][[1]];
Graphics[{Circle[], {Line[{P, F}],
Line@{A, B}}, {Point[{O, P, F}]}}, PlotRange -> 1.2]]
n = 40;
lines = Table[line[t], {t, Subdivide[2 π, n]}];
Manipulate[Show[lines[[1 ;; j]], plot[j/(2 π)]], {j, 0, n, 1}]
Edit-0
Need to be improve since the RegionIntersection
too slow in my cases.
F = Rationalize[{0.4, -0.5}, 0];
o = {0, 0};
P[t_] = AngleVector[t];
e[t_] = λ*P@t /.
Simplify[
Solve[{EuclideanDistance[λ*P@t, F] == 1 - λ,
0 <= λ <= 1, 0 <= t <= 2 π}, λ, Reals],
0 <= t <= 2 π][[1, 1]];
c[t_] = Midpoint[{P@t, F}];
plot[t_] =
Graphics[{Circle[], Point[{o, F}], Point[{e@t, c@t}], Red,
Point[P@t], {Green,
Line[{{e@t, o}, {e@t, P@t}, {e@t, F}}]}, {Dotted,
Line[{F, P[t]}]}}];
n = 40;
lines = Table[
Graphics[{{ColorData["Rainbow"][Rescale[t, {0, 2 π}]],
RegionIntersection[Disk[], InfiniteLine[{c[t], e[t]}]]}}], {t,
Subdivide[2 π, n]}];
Manipulate[
Show[plot[(j - 1)/(2 π)], lines[[1 ;; j]]], {j, 1, n, 1}]