e = 3 (slider from -5 to 5, step 1)
b = 2 (slider from -5 to 5, step 1)
c: x^(2) / (e^(2) + b^(2)) + y^(2) / b^(2) = 1
d: x^(2) / (e^(2) - b^(2)) - y^(2) / b^(2) = 1
F_1: ((-e), 0)
F_2 = Mirror(F_1, yAxis)
A, B, C, S = Intersect(c, d)
xs: sqrt(e^(4) - b^(4)) / e
ys: b^(2) / e
D, E = Intersect(c, yAxis)
F, G = Intersect(c, xAxis)
f = OrthogonalLine(F, xAxis)
g = OrthogonalLine(G, xAxis)
h = OrthogonalLine(E, yAxis)
i = OrthogonalLine(D, yAxis)
H, I = Intersect(d, xAxis)
j = OrthogonalLine(I, xAxis)
k = OrthogonalLine(H, xAxis)
J = Intersect(f, h)
K = Intersect(k, h)
L = Intersect(j, h)
M = Intersect(h, g)
N = Intersect(g, i)
O = Intersect(i, j)
P = Intersect(i, k)
Q = Intersect(f, i)
Vieleck1, p, o, l, k_1 = Polygon(P, O, L, K)
Vieleck2, q, n, m, j_1 = Polygon(Q, N, M, J)
r = Line(P, L)
s = Line(O, K)
t: y = (sqrt(e^(4) - b^(4)) / (e^(2) - b^(2)) * (x - xs)) + ys
a: y = (((-sqrt(e^(4) - b^(4)))) / (e^(2) + b^(2)) * (x - xs)) + ys