r = 2 (slider from 0 to 3, step 0.1)
A = Point(xAxis)
f = OrthogonalLine(A, xAxis)
M = Point(f)
a = OrthogonalLine(M, f)
D = Intersect(a, yAxis)
B = Point(Circle[M, r])
c = Circle(M, B)
P = Point(c)
e = OrthogonalLine(P, a)
C = Intersect(a, e)
O = Intersect(xAxis, yAxis)
b = Segment(M, B)
d = Segment(M, P)
α = Angle(D, M, P)
g = Segment(C, P)
E = Point(Circle[D, α])
h = Segment(D, E)
i = OrthogonalLine(E, a)
j = OrthogonalLine(P, i)
F = Intersect(i, j)
k = Segment(E, F)
l = Segment(A, M)