s

r

s

F = Point(yAxis)
F' = Mirror(F, xAxis)
d = OrthogonalLine(F', yAxis)
B = Point(d)
g = OrthogonalLine(B, d)
J = Point(g)
C = Point(g)
s = Ray(C, J)
O = Intersect(xAxis, yAxis)
A = Point(d)
f = OrthogonalLine(A, d)
h = OrthogonalLine(C, g)
D = Point(h)
E = Point(s)
t1, j, i, c = Polygon(E, C, D)
e = Slope(c)
text1: "eccentricity = " + e + ""
G = Point(f)
k = Line(G, c)
H = Intersect(d, k)
t2, g_1, a, h_1 = Polygon(A, H, G)
p = Circle(F, h_1)
q = Circle(A, H)
I = Intersect(q, f, 1)
l = OrthogonalLine(I, f)
P = Intersect(p, l, 2)
K = Intersect(p, l, 1)
loc1 = Locus(P, G)
loc2 = Locus(K, G)
m = OrthogonalLine(P, l)
L = Intersect(d, m)
n = Segment(L, P)
r = Segment(P, F)