S = (0, 0)
F = (1.72, -0.66)
a = Line(F, S)
b = OrthogonalLine(S, a)
F' = Mirror(F, b)
l = OrthogonalLine(F', a)
c = Parabola(F, l)
P = Point(c)
l_1 = Segment(F, P)
d = OrthogonalLine(P, l)
L = Intersect(l, d)
l_2 = Segment(L, P)
e, f = AngularBisector(l_2, l_1)