Ich schreib mal die Vektoren als Zeilen statt Spalten:
Vektor PX  ist (x1 - p1  ;   x2  -  p2 )  und QX  ( x1 - q1 ; x2 - q2 )
senkrecht: Skalarprodukt = 0
 (x1 - p1)*( x1  -  q1 )    +   ( x2  -  p2 )*(X2 - q2 ) = 0
x1^2  -  x1q1  - p1x1 + p1q1   + x2^2 - p2x2  - q2x2  +  p2q2 = 0 
x1^2  + x1(-q1  - p1)+ p1q1   + x2^2 + x2( -p2  - q2) + p2q2 = 0 
x1^2  + x1(-q1  - p1)  +(1/4)(q1+p1)^2 - (1/4)(q1+p1)^2  + p1q1  
                 + x2^2 + x2( -p2  - q2) +(1/4)(p2+q2)^2 -(1/4)(p2+q2)^2   +  p2q2 = 0 
(x1 -0,5(p1+q1))^2 - (1/4)(q1+p1)^2  + p1q1 + (x2-0,5(p2+q2))^2  -(1/4)(p2+q2)^2   +  p2q2 = 0 
(x1 -0,5(p1+q1))^2 +(x2-0,5(p2+q2))^2  =     (1/4)(q1+p1)^2  -p1q1  +   (1/4)(p2+q2)^2   - p2q2 
(x1 -0,5(p1+q1))^2 +(x2-0,5(p2+q2))^2  =  (1/4)* ((p1-q1)^2 +(p2-q2)^2)
Das gibt einen Kreis um den Mittelpunkt von PQ mit Radius (1/2)LängePQ
Das ist der Satz des Thales !