Inverseur de Peaucellier
Author: Patrice Debart

B = (3.2, -5.12967)
A = (3.2, 4.84302)
dAB = Line(A, B)
r_1 = 3
r_2 = 2
k: sqrt(r_1^(2.0) - r_2^(2.0))
O = (0, 0)
c_1 = Circle(O, r_1 - r_2)
c_3 = Circle(O, r_1 + r_2)
c = Circle(O, r_1)
M = Point(dAB)
c' = Circle(M, r_2)
P = Intersect(c, c', 1.0)
Q = Intersect(c, c', 2.0)
dPQ = Line(P, Q)
M' = Mirror(M, dPQ)
MP = Segment(M, P)
PM' = Segment(P, M')
QM' = Segment(M', Q)
QM = Segment(Q, M)
OQ = Segment(O, Q)
OP = Segment(O, P)
OM = Segment(O, M)
lieuM' = Locus(M', M)