a = 16 (slider from 0 to 20, step 0.5)
b = 9 (slider from 0 to 20, step 0.5)
c = 5 (slider from 0 to 20, step 0.5)
r = 3.2 (slider from 0 to 10, step 0.1)
xy: x^(2) + y^(2) - (a * x + b * y)^(2) / (a^(2) + b^(2) + c^(2)) = r^(2)
d: x^(2) + y^(2) - (a * x - c * y)^(2) / (a^(2) + b^(2) + c^(2)) = r^(2)
A_1 = Intersect(xAxis, yAxis)
B = Point({0, b})
f = Segment(A_1, B)
A = Point({a, 0})
g = Segment(A_1, A)
B_1 = Point({-b, 0})
h = Segment(A_1, B_1)
C = Point({0, -c})
i = Segment(A_1, C)
G = Midpoint(B, A)
H = Midpoint(A, C)
A' = Mirror(A_1, H)
A'_1 = Mirror(A_1, G)
I = Midpoint(B_1, C)
A'_2 = Mirror(A_1, I)
poly1, a_1, b_1, a'_1, c_1 = Polygon(A_1, B, A'_1, A)
poly2, a_2, c_2, a', e_1 = Polygon(A_1, A, A', C)
poly3, a_3, d_1, a'_2, e_2 = Polygon(A_1, B_1, A'_2, C)
j = Segment(A_1, A'_1)
k: x^(2) + y^(2) - (b * x + c * y)^(2) / (a^(2) + b^(2) + c^(2)) = r^(2)
e = Segment(A_1, A'_2)
l = Segment(A_1, A')
F = Point(yAxis)
J = Point(xAxis)
m = Segment(F, J)
K = Point(m)
n = Ray(A_1, K)
L = Intersect(xy, n)
mp1 = Locus(L, K)
M = Point(yAxis)
p = Segment(M, J)
N = Point(p)
q = Ray(A_1, N)
O = Intersect(d, q, 1)
mp2 = Locus(O, N)
P = Point(h)
s = Segment(M, P)
Q = Point(s)
t = Ray(A_1, Q)
R = Intersect(k, t, 1)
mp3 = Locus(R, Q)