O = (0.6, 2.9)
A = (4.1, 2.6)
c = Circle(O, A)
a = Ray(O, A)
Y = Midpoint(O, A)
B = Mirror(Y, A)
b = OrthogonalLine(B, a)
P = Point(c)
d = Line(P, A)
Q = Intersect(d, b)
M = Midpoint(P, Q)