Definition eq_distance A B C D := Py A B A = Py C D C.Definition eq_product A B C D P Q R S := Py A B A * Py C D C = Py P Q P * Py R S R.Definition tangent O1 A O2 B :=
Py O1 O2 O1 * Py O1 O2 O1 / (2+2) + 
Py O1 A O1   * Py O1 A O1 / (2+2) +  
Py O1 B O1   * Py O1 B O1 / (2+2) -
Py O1 O2 O1 * Py O1 A O1 / 2 -
Py O1 O2 O1 * Py O1 B O1 / 2 -
 Py O1 A O1 * Py O1 B O1 / 2.Definition m_ratio Y U V r := 1+r<>0 /\ on_parallel_d Y U U V (r/(1+r)).Definition inversion P Q O A := Py O Q O <> 0 /\ on_line_d P O Q (Py O A O / Py O Q O).Definition eq_angle A B C D E F := 
  S A B C * Py D E F = S D E F  * Py A B C.Definition co_circle A B C D :=
  S C A D * Py C B D = S C B D * Py C A D.

Lemma check_co_circle : forall A B C D, co_circle A B C D <-> eq_angle C A D C B D.