s = 0.5
a = 1
b = 4
c: (x - (2 * a))^(2) - y^(2) = s^(2)
d: x^(2) / a^(2) + y^(2) / b^(2) = 1
Q = Point(d)
e = OrthogonalLine(Q, yAxis)
O = (0, 0)
f = Line(O, Q)
E = Intersect(c, f, 1)
g = OrthogonalLine(E, xAxis)
P = Intersect(g, e)
Ortslinie1 = Locus(P, Q)
Vieleck1, e_1, q, c_1 = Polygon(Q, P, E)