O = (-1.5, 3.2333333333333334)
A = (4.066666666666666, 3.2)
a = Line(O, A)
b = OrthogonalLine(O, a)
B = Rotate(A, 270°, O)
c = Circle(O, B, A)
P = Point(c)
d = OrthogonalLine(P, b)
Q = Intersect(d, b)
e = Circle(Q, Distance[O, P])
M, M' = Intersect(e, d)