Inverseur de Peaucellier
Author: Patrice Debart

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_2 = Circle(O, r_1 + r_2)
c = Circle(O, r_1)
M = (2, 2)
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)
dOM = Ray(O, M)
MP = Segment(M, P)
PM' = Segment(P, M')
QM' = Segment(M', Q)
QM = Segment(Q, M)
Γ = Circle(O, k)