r = 2 (slider from 0 to 3, step 0.1)
A = Point(xAxis)
f = OrthogonalLine(A, xAxis)
M = Point(f)
B = Point(Circle[M, r])
c = Circle(M, B)
P = Point(c)
a = OrthogonalLine(M, f)
D = Intersect(a, yAxis)
O = Intersect(xAxis, yAxis)
b = Segment(M, B)
α = Angle(D, M, P)
d = Segment(M, P)
e = OrthogonalLine(P, a)
C = Intersect(a, e)
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)
β_1: 2 α
H = Intersect(c, a, 2)
m = Segment(M, H)
Q = Rotate(H, 2 α, M)
β = Angle(H, M, Q)
n = Segment(Q, M)
p = OrthogonalLine(Q, i)
I = Intersect(i, p)
q = OrthogonalLine(Q, a)
J = Intersect(a, q)
s = Segment(Q, J)