A = (7.276260056213695, 9.408062128177479)
B = (7.319335721974814, 4.556040004006374)
d = OrthogonalLine(A, Line[A, B])
c = Circle(A, 3 Length[Vector[A, B]])
Z, D = Intersect(c, d)
C = Translate(B, Vector[Vector[A, D]])
poly1, AB, BC, CD, AD = Polygon(A, B, C, D)
M = Point(AB)
N = Intersect(Circle[B, Segment[A, M]], Segment[B, C])
P = Intersect(Circle[C, Segment[A, M]], Segment[C, D])
Q = Intersect(Circle[D, Segment[A, M]], Segment[D, A])
a = 2 (slider from -5 to 5, step 0.1)
α = 2 (slider from -5 to 5, step 0.0010)
β = 4 (slider from -5 to 15, step 1.0E-4)
f(x) = a (x - α)^(2) + β
poly2, m, n, p, q = Polygon(M, N, P, Q)
W: (Length[Vector[A, M] 2 / Length[Vector[A, B]]], poly2 * 12 / poly1)
loc1 = Locus(W, M)