M = Point(xAxis)
c = Circle(M, 1)
P = Point(c)
a = Segment(M, P)
b = OrthogonalLine(P, xAxis)
Q = Intersect(b, xAxis)
d = Segment(Q, P)
α = Angle(xAxis, a)
O = Intersect(xAxis, yAxis)
B = Point(Circle[O, α])
e = Segment(O, B)
f = OrthogonalLine(B, xAxis)
g = OrthogonalLine(P, f)
C = Intersect(f, g)
h = Segment(B, C)