Peaucellier–Lipkin Linkage [gnuplot]

YouTube   Formulation Find the coordinates of all joints \begin{eqnarray*} OA=OP=r\\ PB=PD=a\\ AB=BC=CD=DA...