(* shows that we need A<>B in elim_py_on_line_d_middle Lemma essai : forall A P Q Y : Point, forall r:F, r<> 0 -> on_line_d Y P Q r -> False. *)

Theorem elim_py_on_line_middle :
forall A B P Q Y : Point,
on_line Y P Q ->
Py A Y B = Py A P B + P ** Y / Q ** P * (Py A P B - Py A Q B + 2* Py Q P Q) -
P ** Y / Q ** P * (1 - P ** Y / Q ** P) * Py Q P Q.