eq1: x^(2) / 16 + y^(2) / 25 = 1
α = 1.4311699866353502 rad (slider from 0° to 360°, step 0.017453292519943295)
A: (4, 0)
O: (0, 0)
f = Segment(O, A)
D: (6; α)
g = Segment(O, D)
M = Intersect(eq1, g)
E: (6; (2 * α))
h = Segment(O, E)
N = Intersect(eq1, h)
i = Segment(N, M)
j = OrthogonalLine(O, i)
B = Intersect(j, i)
c = Circle(M, B)